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PART  I:  PLAN  OF  TEST 


by 

D.A.  Hill 


ABSTRACT 

Previous  determinations  of  drag  coefficients  for  both  steady 
and  unsteady  flow  are  reviewed  for  all  reported  experiments  at  critical 
and  supercritical  Reynolds  numbers.  The  unsolved  problems  remaining 
from  earlier  DRES  work  were  found  to  be:  end  effects,  dust  drag  loading, 
the  discrepancies  between  steady  and  unsteady  flow  drag  forces,  and  the 
need  for  a further  data  analysis.  A new  analysis  is  also  presented  for 
the  drag  coefficient  as  a function  of  flow  Mach  number  and  cylinder 
length-diameter  ratio. 

A plan  of  test  is  presented  for  drag  force  measurements  on 
cylinders  during  Event  DICE  THROW,  a 628-ton  AN/FO  explosion  held  in 
October  1976.  The  objectives  of  this  experiment  were  to  provide  drag 
force  support  data  for  the  lattice  mast  and  all  antennas  in  Event  DICE 
THROW,  and  to  extend  the  table  of  data  for  drag  coefficients  to  diameters 
up  to  18  inches.  The  free-flight  measurement  technique,  using  both 
cameras  and  velocity  transducers  for  data  sources,  was  employed  on  seven 
cylinders  at  20,  10  and  7 psi  overpressure  levels.  The  expected  motion 
of  each  cylinder  was  calculated,  based  on  the  predicted  value  for  the 
drag  coefficient.  Finally,  suggestions  for  the  data  analysis  and  for 
further  research  in  this  area  are  discussed.  ^ 

(U) 
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Part  I:  Plan  of  Test 


by 

D.A.  Hill 


1 . INTRODUCTION 

Event  DICE  THROW  was  a large  blast  trial  conducted  on  October  6, 
1976  at  the  White  Sands  Missile  Range,  New  Mexico,  U.S.A.  The  explosive 
charge  was  628  tons  of  an  ammonium  nitrate-fuel  oil  mixture  (AN/FO)  at 
ground  level,  considered  similar  in  explosive  energy  yield  to  previous 
500-ton  TNT  explosions  (PRAIRIE  FLAT,  DIAL  PACK  and  MIXED  COMPANY).  DRES 
was  requested  by  the  Canadian  Forces  (Appendix  A)  to  design  and  carry 
out  an  experiment  on  further  measurements  of  unsteady-flow  drag  forces  on 
cylinders.  The  plan  of  test  for  this  experiment  is  described  separately 
as  Part  I of  this  report  because  the  author  was  relocated  at  Defence 
Rerearch  Establishment  Ottawa  prior  to  the  DICE  THROW  event.  The  execution 
and  analysis  of  the  experiment  were  carried  out  by  other  DRES  workers. 

These  aspects  of  the  experiment  are  described  by  Dr.  A.W.  Gibb  in  Part  II 
of  the  report. 
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2 . REVIEW  OF  UNSTEADY  DRAG  FORCES  ON  CYLINDERS 
2.1  General  Introduction 

When  a fluid  flows  in  a given  direction  towards  and  around 
an  object,  the  force  on  the  object  in  the  flow  direction  is  defined 
as  the  "drag  force"  and  any  force  in  a perpendicular  direction  is  de- 
fined as  "lift".  The  classic  text  in  the  area  of  drag  is  Hoerner's 
book  (1958),  "Fluid-Dynamic  Drag",  which  summarizes  virtually  all  of 
the  work  prior  to  that  date.  The  drag  forces  are  further  divided  into 
skin  friction  drag,  due  to  tangential  stresses  in  the  fluid,  and  form 
or  pressure  drag,  due  to  normal  stresses.  From  dimensional  analysis 
(Owczarek,  1968)  the  drag  force  can  be  shown  to  be  given  by  Cp  (%>U2)A 
where  A is  the  cross-sectional  area  of  the  object  presented  perpendi- 
cularly to  the  flow,  U is  the  fluid  speed,  p the  fluid  density  and  C^ 
the  dimensionless  drag  coefficient.  It  can  also  be  shown  (Owczarek, 
1968)  that  for  a set  of  objects  of  the  same  shape  and  orientation  but 
arbitrary  size,  CQ  = C0(M,R,y)  where  M and  R are  the  Mach  and  Reynolds 
numbers  for  the  flow  and  y is  the  specific  heat  ratio  of  the  gas.  Note 
that  surface  roughness  on  the  object  effectively  changes  the  "shape" 
and  therefore  the  CQ  value.  The  form  drag  force  can  be  computed  by 
knowing  the  pressure  on  each  point  on  the  surface  of  the  object  and 
summing  up  the  component  forces.  This  technique  has  been  used  to  cal- 
culate the  drag  force  on  a cylinder  by  Bishop  and  Rowe  (1967).  As 
the  flow  speed  increases  so  that  R becomes  >RC(RC>  critical  Reynolds 
number,  = (3-5)  x 105)  the  drag  coefficient  changes  from  1.2  to  0.3 
(Owczarek,  1968).  This  can  be  explained  by  the  change  in  the  pres- 
sure distribution  around  the  cylinder  (Gowen  arid  Perkins,  1953,  figure 
5)  as  the  flow  pattern  becomes  turbulent.  For  discussions  of  the 
relationship  between  the  laminar  and  turbulent  boundary  layers  and  the 
von  Karman  vortex  street,  see  Hoerner  (1958),  Owczarek  (1968)  or 
Martin  et  al.  (1965). 

The  main  difference  in  drag  force  between  the  flow  behind 
a shock  front  and  that  in  a wind  tunnel  is  in  the  early  stage  of  the 
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flow  as  the  shock  front  diffracts  around  the  object.  Shadowgraphs 
and  drawings  showing  this  process  are  given  by  Martin  et  al.  (1965  and 
1967)  and  by  Mellsen  and  Naylor  (1969).  The  effect  of  diffraction  on 
the  drag  force  during  the  first  few  milliseconds  of  flow  is  sketched  by 
Wolkovitch  (1968,  figure  4).  The  net  impulse  imparted  to  a cylinder  by 
this  diffraction  phase  has  been  fitted  empirically  in  figure  3.1  of  Long 
et  al.  (1975),  a result  used  to  estimate  the  motion  of  the  cylinders  for 
DICE  THROW. 

2.2  Previous  PRES  Work 

Experiments  similar  to  the  present  one  have  been  carried  out  in 
previous  large  scale  blast  trials.  Those  in  which  the  free-flight  measure- 
ment method1  was  used  include  PRAIRIE  FLAT  (Mellsen,  1969a),  DIAL  PACK 
(Mellsen,  1971)  and  MIXED  COMPANY  (Mellsen,  1974).  The  results  for  the 
circular  and  square  cross-section  cylinders  tested  in  those  experiments 
are  presented  in  Table  1.  In  each  case  the  flow  started  at  supercritical 
Reynolds  number,  R>Rc  (Rc  = (3-5)  x 10s)  and  usually  progressed  to  subcriti- 
cal  values.  In  each  case  average  drag  coefficients  were  reported,  even 
though  large  changes  in  CQ  are  known  to  occur  in  steady  flow  over  this 
Reynolds  number  range.  In  addition,  for  overpressure  levels  above  12  psi 
the  compressibility  effects  due  to  M>Mci  = 0.48  are  obscured  by  the  use  of 
a single  average  drag  coefficient.  The  use  of  average  drag  coefficients 
to  cover  important  changes  in  M and  R values  is  less  desirable  than  a more 
detailed  analysis. 

2.3  Review  of  Unsolved  Problems 

The  values  of  the  drag  coefficient  are  known  to  be  a function 
of  the  four  variables:  M,  R,  L/D  and  steady  vs  unsteady  flow.  For  the 

l 

In  the  free-flight  method,  the  cylinder  is  permitted  to  accelerate 
through  the  air  in  response  to  the  blast  wave.  Acceleration 
(hence  drag  force)  is  derived  as  the  slope  of  the  measured 
Yelocity-time  history  and  curvature  of  the  measured  displacement- 
time history  of  the  cylinder. 
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regions  of  those  variables  of  interest  there  exist  no  complete  or  com- 
prehensive sets  of  data.  The  available  limited  data  will  be  analysed 
and  presented  below  to  indicate  some  of  the  scientific  problems  associated 
with  the  present  state  of  knowledge  of  drag  forces  on  cylinders  for  the 
flow  conditions  relevant  to  DICE  THROW. 

The  measured  Cp  values  from  finite  length  cylinders  in  DIAL 
PACK  (Mel  1 sen,  1971)  have  been  examined  as  a function  of  Mach  number 
(M)  and  compared  to  other  values  reported  for  infinite  length  cylinders 
in  Table  2.  Although  the  data  are  rough,  one  can  estimate  the  aver- 
age value  of  the  ratio 

Cn(M,  finite  L/D,  unsteady  flow) 

n(M)  = -2 

Cp(M,°°  L/D,  steady  flow) 

The  average  values  of  n1  = 0.65  ± 0.18  for  M>Mci  and  n ' = 1.9  + 0.35  for 

M<M  will  be  used  later.  The  importance  of  M , to  the  Cn  value  is 
ci  ci  D 

shown  in  Figure  1,  where  CfJ (M ) is  plotted  for  different  D/L  and  R 

values.  The  dramatic  change  in  the  value  of  Cn  at  M=M  is  due  to  a 

D ci 

change  in  the  flow  pattern  around  the  cylinder  resulting  from  the  local 
flow  speed  becoming  supersonic  at  the  top  and  bottom  of  the  cylinder 
cross-section. 

Figure  2 presents  the  same  data  replotted  as  Cp(D/L)  for 

different  (M,  R)  values.  It  can  be  seen  that  for  M<M„,  and  R<R  the 

ci  c 

classical  results  show  that  L/D  ratio  in  the  range  of  5-20  gives  a Cp 
figure  appreciably  different  from  Cp ( L/D  = °>) . 

The  above  results  are  summarized  in  Table  3.  Different 
regions  of  Mach  and  Reynolds  numbers  are  marked  off  by  the  critical 
values  of  M = 0.48,  Mc?  = 0.9  and  Rc  = (3-5)  x 105.  Only  the  regions 
labelled  I - 1 1 1 are  of  interest  to  the  present  experiment  but  the  re- 
sults for  the  other  regions  are  presented  for  completeness  and  for  com- 
parison purposes. 

Using  the  above  analysis  as  a reference  point  it  was  con- 
cluded that  there  remain  several  important  unsolved  problems  from 
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previous  DRES  work.  These  are: 

2.3.1  End  effects.  The  effect  of  using  finite  length  cylinders  ; 

has  not  been  taken  into  account.  The  previous  use  of  side  plates  near 

to  and  parallel  to  the  cylinder  end  is  of  no  help  as  the  end  plate  j 

must  make  complete  contact  with  the  cylinder  to  avoid  the  end  effects. 

The  severity  of  the  effect  is  estimated  by  the  values  of  n" 
calculated  in  Table  3 (n"  = 1 implies  end  effects  are  negligible)  where 
n"  = 1.5  (Region  I)  and  n"  = 0.65  (Region  III).  These  rough  figures 
serve  to  indicate  that  end  effects  must  be  investigated  more  carefully. 

2.3.2  Unsteady  vs  steady  flow.  In  previous  reports  (Naylor  and  Mel  1 sen, 

1973;  Mellsen,  1969b)  it  was  concluded  that  the  steady-flow  values  for  Cp 
agreed  reasonably  well  with  the  unsteady-flow  results.  The  ratio  4>  of 
Table  3 indicates  the  difference  in  the  two  values  for  finite  length 
cylinders  (<)>  = 1 implies  that  Cp  (steady)  = Cp  (unsteady)).  The  values  <j> 

= 1.2  ± 0.3  (Region  III)  and  $ = 1.9  ± 0.5  (Region  I)  are  estimated  from 
the  data  available.  It  is  evident  that  more  measurements  with  reduced 
end  effects  and  further  analysis  are  required  for  this  problem,  particu- 
larly for  Region  I.  j 

2.3.3  Dust.  The  presence  of  dust  in  the  flow  behind  the  blast  wave, 
especially  after  about  50  ms  of  flow,  has  usually  been  noted  on  previous 
camera  data.  This  causes  an  unknown  amount  of  additional  drag  loading, 

due  to  a larger  effective  density  for  the  "fluid"  flowing  past  the  cylinder. 

Since  dust  would  not  be  present  for  masts  and  antennas  at  sea,  the  dust 
drag  loading  should  be  measured  and  subtracted.  An  attempt  would  be  made 
to  condition  the  soil  in  front  of  the  cylinders  for  reducing  dust.  Note 
that  n'  = 1.9  ± 0.4  for  Region  I could  be  an  indication  of  dust  loading 
as  well  as  the  problem  of  end  effects  and  of  steady  vs^  unsteady  flow. 

2.3.4  Complete  data  analysis.  No  realistic  error  estimates  have  been 
put  on  the  Cp  values  at  each  point  in  the  flow,  nor  on  the  averaged  values. 

Without  these,  subsequent  blast  loading  predictions  lose  precision.  The 
differences  in  the  measured  camera  and  transducer  data  (up  to  30%  either 
way)  have  not  been  explained,  nor  has  it  been  shown  whether  or  not  they 
are  within  the  measurement  error  for  each  technique. 
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3.  SELECTION  OF  PROJECT  OBJECTIVES 

After  consideration  of  the  needs  of  the  Canadian  Forces, 
and  of  the  state  of  our  basic  knowledge  of  drag  forces  on  cylinders 
from  blast  waves,  the  following  goals  were  selected  for  the  DICE  THROW 
drag  project.  The  abbreviated  titles  such  as  "support”  and  "table" 
are  intended  for  use  in  Tables  4 and  5. 

(i)  Support 

Provide  drag  force  data  in  support  of  experiments  on  the 
lattice  mast,  the  VHF/UHF  Polemast  antenna  and  the  whip  antenna,  i.e., 
test  cylinders  with  diameters  similar  to  the  antenna  diameters. 

(ii)  Table 

Extend  the  current  table  of  data  of  unsteady  flow  drag  coef- 
ficients to  larger  cylinder  diameters,  for  effectively  infinite  length 
cylinders.  This  is  for  both  future  potential  applications  and  scien- 
tific interest,  and  applies  both  to  the  diffraction  phase  loading  and 
to  the  pressure  drag  forces. 

(iii)  End  Effects 

Measure  the  differences  between  finite  length  and  effectively 
infinite  length  cylinders  both  for  this  project  and  for  correlation  to, 
and  possible  correction  of,  previous  DRES  work. 

(iv)  Dust 

Estimate  the  extra  drag  force  due  to  dust  for  the  purpose  of 
subtracting  it  from  the  total  measured  drag  force. 

( v)  Unsteady  Flow  versus  Steady  Flow 

Obtain  more  data  from  all  of  the  cylinders  for  comparison 
of  the  steady  flow  and  unsteady  flow  values  of  the  drag  coefficient  at 
the  Mach  and  Reynolds  numbers  used. 

(vi)  Diffraction 

Collect  data  for  a contribution  to  the  scientific  study  of 


UNCLASSIFIED 


UNCLASSIFIED  /9 

diffraction  for  effectively  infinite  length  cylinders  in  Regions  I and 
III  of  Table  3. 
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4.  THE  EXPERIMENTAL  DESIGN 

4.1  Selection  of  Cylinders 

From  considerations  of  both  goals  and  practical  limitations, 
the  methods  and  measurements  described  below  were  selected. 

The  free-flight  method  was  used  again  because  it  had  proven 
successful  on  previous  trials. 

Since  the  blast  trial  is  a one-shot  experiment  with  a mini- 
mum of  equipment  proofing  beforehand,  as  few  changes  as  possible  were 
made  to  the  method.  Only  those  changes  necessary  for  correcting  pre- 
vious problems  or  for  meeting  a new  cylinder  diameter  requirement  were 
carried  out. 

The  problem  of  end  effects  has  been  discussed.  One  of  the 
solutions  considered  was  to  add  a hemisphere  to  each  end  of  a finite 
length  cylinder.  Since  the  CQ  value  for  a hemisphere  is  well  known, 
the  idea  would  be  to  subtract  off  the  "known"  end  effect  and  deduce 
the  drag  force  per  unit  length  on  the  effectively  infinite  length 
cylinder  in  the  middle.  However,  the  drag  coefficient  of  the  sum  of 
the  parts  is  unlikely  to  equal  the  sum  of  the  drag  coefficients  of  the 
parts,  because  the  flow  patterns  of  a hemisphere  and  a cylinder  do  not 
join  smoothly,  i.e.,  the  cylinder-flow  vortices  have  no  counterpart  in 
the  flow  pattern  around  the  hemisphere. 

An  alternate  solution  to  the  problem  of  end  effects  is  to 
attach  thin  perpendicular  end  plates  to  the  finite  length  cylinder. 

The  truncated  flow  patterns  here  would  be  closer  to  the  ideal  effec- 
tively-infinite  length,  especially  during  the  diffraction  phase  of 
the  loading.  This  design  also  more  closely  approximates  a finite- 
length  cylinder  joined  at  each  end  to  other  perpendicular  members, 
with  or  without  gusset  plates.  By  changing  the  end  plate  thickness 
(causing  negligible  additional  drag  forces)  one  can  adjust  the  total 
cylinder  weight  and  hence  the  cylinder  range  of  motion  to  permit 
optimum  measurement  with  tne  velocity  transducers. 
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One  disadvantage  of  the  end  plates  is  that  in  a sidewise 
blast-force  anomaly  such  as  occurred  at  the  25  psi  level  in  MIXED  COM- 
PANY (Mellsen,  1974)  the  force  strikes  the  end  plates  at  a glancing  angle 
causing  diagonal,  possibly  oscillatory,  motion  of  the  whole  cylinder. 

There  could  also  be  a problem  if  a turbulent  boundary  layer  flow  is  set 
up  over  the  end  plate,  causing  excessive  end  plate  drag.  The  worst  case 
is  estimated  below  to  be  a 10%  effect. 

From  examination  of  shadowgraphs  of  the  vortex  shedding  process, 
(Naylor  and  Mellsen,  1973,  figures  10B  and  C)  an  end  plate  extension  of 
1.5  diameters  at  the  end  of  the  cylinder  appears  desirable.  However,  to 
limit  the  problems  of  end  plate  drag  and  vibration  an  extension  of  only 
1 diameter  was  chosen. 

The  final  selection  of  cylinder  diameters,  end  plate  diameters 
and  overpressure  levels  is  shown  in  Table  4<  The  same  length-diameter 
ratio  (L/D)  was  chosen  for  cylinder  numbers  1,  3,  4 and  5 so  that  the 
effects  of  end  plates  and  Mach  number  regions  for  fixed  L/D  could  be 
studied.  The  value  of  L/D  = 5 was  chosen  to  keep  the  length  of  the 
largest  cylinders  at  a reasonable  value,  7.5  ft.  The  known  end  effects 
are  not  appreciably  changed  when  the  L/D  ratio  is  changed  from  5 to  10 
(Hoerner,  1958,  Section  3,  figure  28).  Table  5 lists  the  ways  in  which 
each  cylinder  contributed  to  the  goals  of  the  project. 

By  restricting  cylinder  diameters  to  three  values,  more  overlap 
in  the  measurements  was  possible.  The  9.5  in.  diameter  was  included  to 
support  the  VHF/UHF  antenna  project,  as  there  were  no  DRES  drag  measure- 
ments of  a similar  diameter  at  10  psi. 


For  the  smaller  diameter,  a value  of  3.5  in.  was  chosen  as  being 
equivalent  to  the  larger  diameter  vertical  members  to  be  used  in  the  test 
lattice  mast,  and  as  a reasonable  value  to  use  for  correlation  with  the 
whip  antenna,  which  is  tapered  from  2.5  in.  to  5.5  in. 

For  drag  data  in  support  of  the  structures  projects,  the  drag 
information  must  cover  the  first  quarter  period  of  the  fundamental  mode  of 
oscillation  of  the  related  structure.  From  previous  DRES  work  and  new 
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estimates,  this  period  was  estimated  to  be  about  70  ms  or  less  for  the 
VHF/UHF  antenna,  < 50  ms  for  the  lattice  mast  (25  - 30  ms  previously 
(Long  and  Laidlaw,  1973))  and  < 100  ms  for  the  18  fn.  cylinder  (86  ms 
for  the  DIAL  PACK  TACAN  mast  (Coffey,  1971)). 

For  completing  the  table  of  data,  the  basic  scaling  para- 
meters for  fluid  dynamic  drag  are  known  to  be  M and  R.  But  for  a series 
of  near-identical  equivalent  charge-weight  explosions  such  as  we  have 
in  the  500-ton  series,  then  the  psi  level  and  diameter  are  an  equivalent 
pair  of  variables.  These  relationships  are  shown  in  Figure  3 where  our 
selection  of  cylinders  is  also  marked. 

This  figure  shows  the  regions  below  and  above  the  critical 
M value  of  0.48  and  the  critical  R value  of  5 x 105.  Each  straight 
line  through  the  origin  represents  the  range  of  M and  R values  traversed 
by  one  cylinder  after  the  explosion.  The  straight  lines  through  the 
origin  make  the  simplifying  assumption  that  M R,  that  is,  they  assume 
that  the  speed  of  sound,  viscosity  and  density  changes  involved  in  the 
shock  wave  approximately  cancel  out.  This  is  an  over-simplification  but 
is  sufficient  for  the  purposes  of  this  drawing. 

Each  straight  line  through  the  origin  represents  a different 
diameter  and  the  latter  are  recorded  on  the  drawing  in  inches.  A1 1 the 
other  ranges  of  M and  R covered  in  literature  available  to  us  are  also 
recorded  for  R > 5 x 105  or  R > 1 x 105  and  M > 0.48.  The  legend  ex- 
plains the  meaning  of  the  curves,  and  the  reference  letters  A-J  are  the 
same  references  as  listed  in  Table  3. 

The  paths  to  be  taken  by  our  seven  cylinders,  as  well  as  the 
starting  position  for  each  cylinder,  are  also  marked. 

4.2  Selection  of  Instrumentation 

A1 1 cylinders  were  instrumented  with  a velocity  transducer  at 
each  end,  and  were  photographed  with  a high-speed  camera  (framing  rate 
of  order  1000  frames/second).  The  transducers  on  each  end  were  cali- 
brated and  recorded  separately. 
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The  overlap  between  cameras  and  transducers  was  felt  to  be 
necessary  because  experience  has  shown  that  either  recording  device  may 
fail  when  subjected  to  the  violent  air  blast  loading  from  the  explosion. 
The  position  where  each  cylinder  first  strikes  the  ground  was  marked 
and  measured  as  a further  fail-safe  measurement.  This  one  measurement 
gives  the  total  impulse  received  by  the  cylinder  from  the  blast  if  negli- 
gible lift  is  assumed. 

A second  reason  for  the  dual  recording  system  is  that  each 
technique  gives  secondary  information  unique  to  itself.  That  is: 

(i)  velocity  transducers 

- record  from  each  cylinder  end  for  angle  of  yaw, 

- measure  degree  of  cylinder  vibration,  and 
frequency  of  vibration. 

(ii)  cameras 

- record  if  a cylinder  end  falls  off  the  mount  due  to 
ground  shock,  as  it  did  for  one  cylinder  about  100 
ms  before  time  of  arrival  for  MIXED  COMPANY, 

- indicate  through  image  obscuration  the  presence  of 
dust  in  the  blast  wave, 

- observe  any  flying  debris  hitting  the  cylinders, 

- permit  vertical  and  horizontal  motions  and  rotational 
motion  about  cylinder  axes  to  be  measured, 

- record  time  of  arrival  of  shock  (red  ribbon  flutter). 

The  vertical  component  of  motion  could  be  measured  in  order  to 
check  for  possible  lift  due  to  non-horizontal  blast  flow  direction  and/or 
a choking  effect  due  to  the  flow  between  a large  cylinder  and  the  ground. 

4.3  Predicted  Motion  of  Cylinders 

As  each  cylinder  was  a model  for  a particular  structural  member 
(whether  represented  in  DICE  THROW  or  not),  it  was  necessary  that  the 
duration  of  recorded  motion  for  each  cylinder  be  not  less  than  the  first 
quarter  period  of  the  fundamental  mode  of  oscillation  of  the  particular 
structure.  Each  cylinder  was  designed  with  a mass  large  enough  to  cause 
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the  predicted  motion  to  fall  within  the  desired  distance  of  travel 
over  the  desired  time  interval. 

The  Cp  values  chosen  for  these  calculations  and  the  detailed 
calculations  themselves  are  given  in  Appendix  B. 

The  internal  structure  of  each  cylinder  was  designed  to  be 
more  rigid  against  transverse  vibrations  than  previous  designs,  and 
to  have  sufficiently  strong  walls  to  withstand  collapse  when  hit  with 
the  full  initial  force  of  the  shock  front  (see  below). 

The  estimates  are  conservative  by  20 % in  the  Cp  values  used 
and  by  20%  in  the  desired  recording  time.  This  40%  safety  factor  was 
felt  to  be  sufficient  to  protect  against  likely  errors  in  Cp  values  and 
against  the  possibility  of  an  unexpectedly  strong  shock  front  from  the 
628-ton  AN/FO  charge. 

4.4  Mechanical  Design  of  Cylinders 

The  overall  arrangement  of  cylinders  and  supporting  equipment 
is  illustrated  in  Figure  4,  the  assembly  photo.  Most  of  the  equipment 
is  similar  to  that  described  in  the  DIAL  PACK  report  (Mellsen,  1971) 
with  a few  changes  as  described  below.  The  outside  diagonal  support 
plates  were  designed  to  minimize  the  possibility  of  side-plate  vibrations 
and  yet  keep  cross  bars  out  of  the  centre  flow  region  between  the 
cylinder  and  the  ground.  In  fact,  there  will  be  a certain  amount  of 
flow  choking  in  that  region.  The  worst  case  is  for  the  18  in.  diameter 
cylinders  centred  6 ft.  off  the  ground  which  leave  a separation  gap  of 
only  63  ins.  or  3.5  diameters  between  the  cylinder  and  the  ground. 

Assuming  a flat  solid  ground  surface  and  a vertical  original  shock  front, 
the  method  of  flow  images  can  be  used  along  with  previous  DRES  "soli- 
dity-effect" measurements  to  estimate  the  size  of  the  ground  interference 
effect  on  the  Cp  value.  The  flow  pattern  between  the  idealized  flat 
solid  ground  and  the  cylinder  3.5  diameters  above  it  is  regarded  as 
identical  to  the  flow  pattern  between  two  isolated  cylinders  of  7 diameter 
separation.  Approximate  DRES  data  (Long  and  Laidlaw,  1973,  Table  2.1) 
for  a related  experiment  indicate  that  the  effect  for  this  separation  is 
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likely  not  more  than  10%.  The  flow  will,  however,  produce  a small  lift 
force  on  the  cylinder  which  will  be  detected  through  the  camera  measure- 
ments of  the  vertical  motion  of  the  cylinder.  Another  reason  for  the 
5-6  ft.  height  of  the  cylinders  above  ground. is  to  reduce  the  dust  drag 
loading. 

The  camera  marker  plate  served  the  dual  purpose  of  providing 
an  upper  scale  marked  off  in  inches  to  record  the  cylinder  position,  and 
a lower  scale  to  provide  black/white  contrast  for  possible  use  in  esti- 
mating relative  dust  density  as  a function  of  time.  The  camera  marker 
plate  and  the  support  plates  were  9 in.  from  the  cylinder  so  as  not  to 
impede  the  flow  of  air  around  the  end  of  the  cylinder.  Because  of  this, 
the  camera  lens  was  adjusted  to  have  both  the  scale  and  the  cylinder  in 
focus  at  the  same  time,  and  a constant  geometrical  scale  factor  was 
applied  to  the  position  measurements  to  relate  them  to  true  cylinder 
motion. 

The  cylinders  were  connected  to  the  velocity  transducers  by  a 
system  modified  from  MIXED  COMPANY.  The  connecting  rod  was  longer  to  re- 
duce the  error  in  the  horizontal  speed  measurement  due  to  the  approximately 
2 in.  drop  in  height  under  the  force  of  gravity  which  occurs  in  the  first 
100  ms  of  flight.  The  flexible  aluminum  segment  of  the  connector  reduced 
the  chance  of  broken  magnets  resulting  from  yaw  in  the  cylinder  motion. 

For  the  same  reason,  the  gimbal  mounts  were  farther  from  their  support 
plates  than  in  previous  trials. 

The  internal  structure  of  one  of  the  cylinders  is  depicted  in 
Figure  5.  The  internal  support  discs  and  the  1 in.  rod  through  the  centre 
were  to  strengthen  the  structure  against  the  transverse  vibrations 
observed  to  interfere  with  the  velocity  measurements  on  all  previous 
similar  trials.  They  also  provide  strength  against  the  possible  collapse 
of  the  front  cylinder  wall  when  hit  by  the  full  force  of  the  initial 
shock  front.  Pertinent  design  calculations  are  contained  in  Appendix  C. 

As  some  transverse  oscillation  of  these  cylinders  was  expected,  the  funda- 
mental vibration  frequency  of  each  cylinder  with  its  end  plates  was 
measured  prior  to  the  trial.  There  is  some  evidence  that  in  certain  pre- 
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vious  blast  trials  the  observed  oscillations  grow  and  diminish  in  a 
way  that  may  be  analogous  to  beats  between  the  fundamental  oscillations 
and  the  forcing  function  of  the  frequency  of  vortex  shedding  which  is 
related  to  the  Strouhal  number. 

4. 5 The  Associated  Dust  Density  Experiment 

The  drag  loading  from  dust  has  been  estimated  to  be  as  high 
as  20%  of  the  total  force  by  Naylor  and  Mellsen  (1973)  but,  since  it 
has  never  been  measured  accurately,  it  may  be  larger.  Consequently, 
the  ground  for  100  ft  in  front  of  the  cylinders  was  treated  to 
minimize  the  presence  of  dust,  and  one  or  two  simple  passive  dust 
collectors  were  installed  at  each  overpressure  level.  It  was 
hoped  that  they  would  provide  an  estimate  of  the  total  mass  of  dust 
in  the  shock  wave  and  the  particle  size  distribution. 
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5.  DATA  ANALYSIS 

5.1  Error  Analysis 

The  data  analyses  of  both  velocity  transducer  and  camera  data 
were  similar  to  those  employed  for  data  from  the  MIXED  COMPANY  Event 
(Mellsen,  1974).  However,  in  previous  analyses,  there  has  never  been  a 
discussion  of  the  errors  in  the  Cp  values  calculated  for  a given  cylinder 
at  each  point  in  time.  In  the  DICE  THROW  analysis,  the  effect  of  noise 
and  fluctuations  in  the  raw  data  was  reflected  in  a quoted  error  in  the 
Cp  values  for  each  method.  This  permitted  a realistic  appraisal  of  the 
consistency  of  results  from  camera  and  transducer  data. 
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6.  SUGGESTIONS  FOR  FURTHER  RESEARCH 

Many  factors  are  known  to  affect  the  C^  values,  but  not 
all  have  yet  been  investigated.  In  most  cases,  it  would  take  too 
many  experiments  to  map  out  a set  of  data  of  limited  importance.  Con- 
sequently, the  following  effects  were  not  examined  in  DICE  THROW: 

(a)  Turbulence 

The  degree  of  turbulence  behind  the  shock  front  is 
known  to  be  important,  especially  for  multiply-reflected 
waves. 

(b)  Surface  Roughness 

The  effect  of  varying  degrees  of  cylinder  surface  rough- 
ness has  been  shown  to  be  important  in  wind  tunnel  work. 
Despite  the  fact  that  real  masts  will  have  rough  surfaces  due 
to  ice  covering  and  general  weathering,  this  is  considered  to 
be  too  complex  a problem  for  the  current  program.  The  masts 
and  antennas  tested  were  aerodynamical ly  "smooth";  for  that 
reason,  smooth  cylinders  were  used  in  the  aerodynamic  drag 
experiment. 

( c ) Shading  and  Sol idity 

"Shading"  and  "solidity"  are  two  interference  effects  on 
the  mast  structure.  "Shading"  refers  to  mutual  interference 
effects  of  two  or  more  members  in  a plane  parallel  to  the  flow 
direction;  "solidity"  is  a similar  effect  for  members  in  a 
plane  perpendicular  to  the  flow  direction.  There  exist  some 
simple  DRES  results  of  limited  generality  but  it  is  recognized 
that  more  data  are  needed. 
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PART  II:  TEST  RESULTS  AND  ANALYSIS 


by 

A.W.M.  Gibb 


ABSTRACT 

Results  are  presented  from  an  experiment  to  measure  aero- 
dynamic drag  on  circular  cylinders  under  unsteady  flow  conditions  in 
a long  duration  free-field  blast  wave  during  Event  DICE  THROW.  These 
results  provided  drag  loading  information  required  for  analysis  of  the 
structural  response  tests  on  Canadian  Navy  masts  and  antennae  also 
conducted  during  Event  DICE  THROW.  Seven  cylinders,  distributed  at 
nominal  20,  10  and  7 psi  peak  overpressure  locations  and  spanning 
three  different  diameters  (3.5,  9.5  and  18  inches),  were  studied.  The 
18-inch  diameter  cylinder  at  20  psi  with  48-inch  diameter  end  plates 
was  partially  destroyed  by  a sidewise  blast  pressure  anomaly  travel- 
ling from  east  to  west.  No  useful  data  were  obtained  for  this  cylinder, 
but  the  remaining  six  cylinders  yielded  valid  data.  A free-flight  method, 
developed  in  earlier  trials  (PRAIRIE  FLAT,  DIAL  PACK,  MIXED  COMPANY),  was 
employed  to  measure  time-dependent  drag  pressures.  For  every  cylinder, 
one  velocity  transducer  was  attached  to  each  end  of  the  central  shaft 
to  record  cylinder  velocity  vs  time,  while  a high-speed  camera  recorded 
displacement  vs  time.  Cylinder  acceleration,  and  hence  drag  pressure, 
was  obtained  from  the  slope  and  curvature,  respectively,  of  these  curves. 
Generally  good  agreement  was  obtained  between  results  derived  from  camera 
and  transducer  data.  Dynamic  pressure  (needed  to  extract  drag  coeffi- 
cients) was  calculated,  assuming  a Friedlander-type  overpressure  decay, 
from  ground-level  gauge  measurements  of  overpressure-time  histories  at 
the  20,  10  and  7 psi  peak  overpressure  locations.  Some  cylinders  were 
fitted  with  extended  end  plates  to  reduce  end  effects.  Comparison  of  re- 
sults for  cylinders  with  and  without  extended  end  plates  indicated  the 
presence  of  substantial  end  effects  at  critical  and  supercritical  Reynolds 
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numbers.  Dust  samples  were  collected  at  each  cylinder  location  on 
vertical  aluminum  channels  filled  with  grease.  These  samples,  com- 
bined with  camera  records,  suggest  that  dust  loading  was  insignifi- 
cant at  the  initial  cylinder  positions  5 or  6 feet  above  ground. 
Measured  drag  coefficients  for  Mach  number  <0.4  were  in  agreement  with 
steady-state  values  for  Reynolds  numbers  in  the  range  (4  - 30)  x 1 05 , 
but  were  lower  than  steady-state  values  in  the  range  (30  - 40)  x 105. 

(U) 
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PART  II:  TEST  RESULTS  AND  ANALYSIS 


A.W.M.  Gibb 


I.  INTRODUCTION 

Event  DICE  THROW,  a free-field  blast  trial  employing  a 628-ton 
Ammonium  Nitrate/Fuel  Oil  explosive  charge,  took  place  on  October  6,  1976 
at  White  Sands  Missile  Range,  New  Mexico.  One  of  the  Canadian  projects, 
the  measurement  of  aerodynamic  drag  on  right  circular  cylinders,  using 
the  free-flight  method,  is  the  subject  of  this  report.  The  project  was 
undertaken  at  the  request  of  Director  Maritime  Facilities  and  Resources 
to  provide  blast-loading  information  for  the  lattice  mast  (Laidlaw,  1977), 
pole  mast  (Price  and  Coffey,  1977)  and  whip  antennae  (Price,  1977)  which 
underwent  structural  response  testing  during  this  trial. 

This  project  was  a continuation  of  research  begun  in  Operation 
PRAIRIE  FLAT  (Mellsen,  1969a)  and  continued  in  Events  DIAL  PACK  (Mellsen, 
1971  and  Naylor,  1973)  and  MIXED  COMPANY  (Mellsen,  1974).  The  choice  of 
cylinder  diameter,  peak  overpressure  location,  and  cylinder  design  for 
Event  DICE  THROW  are  discussed  in  detail  in  the  Plan  of  Test  which  com- 
prises the  first  part  of  this  report.  A brief  summary  only  will  be  pre- 
sented here  for  the  sake  of  continuity. 

Seven  cylinders  of  circular  cross-section  were  employed.  Their 
basic  properties  are  sutmarized  in  Table  7.  Two  of  the  diameters  em- 
ployed, 3.5  in.  and  9.5  in.,  were  chosen  because  they  correspond  closely 
to  the  diameters  of  the  main  structural  members  of  the  related  structures 
(3.5  in.  - whip  antenna  and  lattice  mast;  9.5  in.  - pole  mast).  The 
third  diameter,  18  in.,  was  included  to  support  future  mast  designs.  The 
cylinders  were  located  at  the  same  peak  overpressure  levels  as  their 
related  structures  (3.5  in.  diameter  at  10  psi,  9.5  in.  diameter  at 

UNCLASSIFIED 


/ 22 


UNCLASSIFIED 


7 psi).  An  additional  3.5  in.  diameter  cylinder  was  located  at  20  psi 
peak  overpressure.  The  major  unresolved  problem  chosen  for  study  in 
this  test  was  the  influence  of  end  effects  for  finite-length  cylinders 
on  the  measured  drag  coefficient.  With  this  goal  in  mind,  the  cylinders 
of  a given  diameter  were  grouped  in  pairs.  In  each  pair,  one  cylinder 
had  end  plates  with  the  same  diameter  as  the  cylinder  diameter;  the 
second  cylinder  had  end  plates  with  a diameter  which  was  3 times  the 
cylinder  diameter.  The  purpose  of  the  extended  end  plates  was  to  eli- 
minate end  effects  by  cutting  off  the  air  flow  around  the  ends  of  the 
cylinder. 

The  methods  of  data  recording  were  the  same  as  those  deve- 
loped and  used  in  previous  trials  employing  the  free-flight  method. 
Velocity  transducers  were  used  on  all  test  cylinders  to  record  cylinder 
velocity  as  a function  of  time.  In  addition,  a high-speed  camera, 
operating  at  approximately  1000  frames/second,  was  stationed  at  each 
cylinder  location  to  record  cylinder  displacement  as  a function  of 
time.  The  slope  of  the  velocity- time  curve  and  the  curvature  of  the 
displacement-time  curve  provided  independent  measurements  of  cylinder 
acceleration,  and  hence  drag  force,  as  a function  of  time.  The  camera 
records  also  provided  secondary  information  on  possible  complicating 
factors  such  as  cylinder  rotation  and  the  presence  of  solids  (both 
fine  dust  or  massive  particles)  in  the  blast  wave. 

All  of  the  measurements  reported  herein  refer  to  the  drag 
phase  of  loading  on  the  cylinder.  No  measurements  of  loading  during 
the  initial  shock  diffraction  phase  are  reported. 
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2.  APPARATUS 

2.1  Test  Cylinders  and  Mounts 

A typical  test  set-up  is  shown  in  Figure  4. 

Figure  6 indicates  the  relative  position  of  each  cylinder  with 
respect  to  Ground  Zero. 

At  each  location,  support  for  the  cylinders  was  provided  by  two 
vertical  rectangular  plates  made  of  0.25  inch  steel,  with  their  bottom 
edges  fastened  firmly  to  a concrete  base.  Additional  rigidity  for  these 
plates  was  provided  by  triangular  support  in  the  form  of  two  one  inch 
diameter  steel  bars  welded  to  the  outside  of  the  support  plates  at  an 
angle  of  30  degrees  approximately  5 feet  above  ground  level.  The  lower 
ends  of  these  bars  were  set  into  the  concrete  base.  This  arrangement 
provided  a more  rigid  support  than  the  arrangement  used  in  previous  trials 
(horizontal  parallel  bars  connecting  the  inner  surfaces  of  the  support 
plates). 

The  design  of  the  test  cylinders  is  discussed  extensively  in  the 
Plan  of  Test  and  the  cylinder  properties  are  summarized  in  Table  9.  The 
construction  of  a typical  cylinder  is  illustrated  in  Figure  5.  Briefly, 
each  was  a right  circular  cylinder  with  a solid  centre  shaft  of  0.75  or 
1.0  inch  diameter  which  extended  14  inches  beyond  each  end  of  the  cylinder. 
Flats  were  cut  in  the  shaft  nine  inches  from  each  end  of  the  cylinder,  and 
the  cylinder  was  suspended  between  the  support  plates  with  the  flats  resting 
on  the  tops  of  the  support  plates.  The  purpose  of  the  flats  was  to  pre- 
vent the  cylinders  from  rolling  off  the  supports  under  the  influence  of 
small  gusts  of  wind  prior  to  the  shot.  The  coefficient  of  sliding  friction 
between  support  plate  and  cylinder  shaft  was  minimized  by  application  of 
silicone  grease  to  the  top  of  the  support  plate. 

There  is  one  additional  feature  of  the  cylinder  design  which  was 
added  at  an  advanced  date  and  was  not  discussed  in  the  Plan  of  Test.  In 
an  attempt  to  dampen  shaft  vibrations  by  providing  better  mechanical  coup- 
ling between  the  central  shaft  and  the  shell  of  each  cylinder,  0-ring  grooves 
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were  cut  in  the  outer  edges  of  the  interior  support  discs  (Figure  5) 
and  3/16  or  1/4  inch  diameter  0-rings  were  inserted.  This  addition 
helped  to  fill  in  any  small  gaps  between  the  inner  surface  of  the 
cylindrical  shell  and  the  outer  edges  of  the  interior  support  discs. 

As  well,  an  additional  plywood  support  disc  was  placed  4 inches  in 
from  each  end  of  the  cylinder.  Wood  was  used  because  of  its  low  den- 
sity; the  intention  was  to  provide  extra  support  for  the  shaft  close 
to  the  cylinder  end  without  seriously  affecting  uniformity  of  mass 
distribution  along  the  cylinder  length.  The  additional  support  close 
to  the  cylinder  end  was  provided  in  an  attempt  to  increase  the  vib- 
ration frequency  of  the  end  of  the  cylinder  shaft;  this  was  felt  to  be 
necessary  because  calculations  had  shown  that  the  fundamental  vibration 
frequency  of  the  solid  centre  shaft  could  be  low  enough  to  complicate 
analysis  of  the  velocity  transducer  data. 

2.2  Velocity  Transducers 

2.2.1  General . The  velocity  transducers  for  measuring  cylinder  velo- 
city directly  as  a function  of  time  consisted  of  seven  pairs  of  Hew- 
lett-Packard Sanborn  7LV9  transducers.  On  a given  cylinder,  two  trans- 
ducers were  used,  one  coupled  by  a mechanical  linkage  to  each  end  of 
the  cylinder  shaft.  The  transducer  signals  were  recorded  separately, 
on  a tape  recorder  with  nominal  DC  to  4KHz  recording  bandwidth.  This 
provided  two  independent  measurements  of  velocity  for  each  cylinder. 

A close-up  view  prior  to  the  shot  showing  the  transducer  coil  in  its 
gimbal  mount,  and  the  mechanical  linkage  which  couples  the  magnet  in- 
side the  coil  to  the  end  of  the  cyl inder  shaft,  is  presented  in  Figure 
7. 

2.2.2  Principle  of  operation.  The  transducers  consist  of  two  toroidal 
hollow-core  wire-wound  coils  in  series,  along  the  conmon  core  of  which 
travels  a permanent  magnet.  The  magnet  is  mechanically  coupled  to  the 
object  whose  velocity  is  to  be  measured.  When  the  object  moves,  pulling 
the  magnet  through  the  coils,  the  flux  lines  of  the  magnet  cut  the 
coils  and  induce  a voltage  in  them  which  is  proportional  to  the  magnet 
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velocity.  The  two-coil  design  gives  rise  to  two  regions  of  equal  length 
over  which  the  transducer  velocity  response  is  uniform,  separated  by  a 
short  region  of  highly  non-uniform  response.  Over  the  second  working 
length,  however,  the  sensitivity  is  reduced  by  a factor  of  two  and  the 
polarity  reversed  with  respect  to  the  first  working  length.  A typical 
constant-velocity  trace  would  have  the  following  shape: 


The  useful  working  lengths  are  denoted  by  the  symbol  L.  This 
diagram  explains  why  the  velocity  traces  from  the  cylinders  shown  in 
Figures  15  to  25  are  in  two  segments  with  a gap  between:  only  data  col- 
lected in  the  working  lengths,  L,  are  usable. 

The  Sanborn  7LV9  transducers  used  in  this  trial  had  two  nominal 
working  lengths,  each  .9  inches.  The  overall  recording  length  of  approxi- 
mately 20  inches  was  sufficient  to  permit  between  80  and  200  milliseconds 
of  cylinder  motion  to  be  recorded. 

2.2.3.  Velocity  calibration.  The  electro-mechanical  calibration  method, 
which  employed  a Kistler  standard  accelerometer  and  shake  table,  has  been 
described  in  a previous  report  (Mellsen,  1971,  Appendix  A).  The  cali- 
bration error  was  estimated  to  be  ± 3%.  In  addition,  careful  inspection 
of  transducer  traces  indicated  possible  variations  of  ± 2%  in  the  uni- 
formity of  response  over  any  working  length.  An  overall  uncertainty  of 
± b%  in  the  velocity  calibration  was  assumed  when  analyzing  the  data. 
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2. 3 Cameras 

Photosonic  framing  cameras,  operating  at  nominal  speeds  of 
1000  frames/second,  were  set  on  camera  posts  at  each  cylinder  location 
to  record  displacement  of  the  cylinders  as  a function  of  time.  A camera 
set  on  its  mounting  post  can  be  seen  in  Figure  8.  Relevant  information 
on  the  cameras  and  their  locations  is  given  in  Table  10. 

The  cameras  provided  internal  timing  marks  which  were  pro- 
jected onto  the  film  at  10  millisecond  intervals  to  permit  the  framing 
rate  to  be  established  and  the  constancy  of  the  framing  rate  over  the 
recording  interval  to  be  checked.  The  timing  mark  generators  functioned 
on  all  cameras.  To  indicate  the  time  of  arrival  of  the  shock  front  in 
the  film  frame,  a red  ribbon  was  glued  to  the  back  edge  of  the  support 
plate.  Horizontal  distance  calibration  was  provided  by  a photomarker 
plate  with  a 12-inch  scale  marked  off  in  inches.  Both  of  these  aids  can 
be  seen  clearly  in  Figure  7. 

Approximately  one  week  before  shot  day.  Test  Command  moved 
the  shot  time  forward  from  1300  hours  to  0800  hours.  This  change  pro- 
vided potentially  serious  problems  for  the  camera  recording  system.  The 
angle  and  position  of  the  sun  at  0800  hours  were  such  that  it  shone 
almost  directly  into  the  camera  lenses.  The  cylinder  ends  to  be  photo- 
graphed were  in  shadow,  and  the  high  background  light  level  caused  ex- 
tremely poor  image  contrast.  To  improve  contrast,  reflector  panels  were 
constructed  of  sheets  of  aluminum  foil  taped  to  a 4 ft  x 4 ft  sheet  of 
plywood  backing.  One  panel  was  bolted  to  each  camera  post,  and  each 
was  oriented  to  reflect  sunlight  onto  the  cylinder  end  to  be  photo- 
graphed. One  such  panel  is  visible  in  Figure  8.  These  hastily- 
constructed  reflectors  provided  sufficient  reflected  light  to  permit 
pictures  of  acceptable  contrast  to  be  recorded  at  shot  time  by  all 
cameras.  However,  a slightly  denser  cloud  cover  at  shot  time  could 
have  ruined  the  camera  experiment  entirely. 


Such  a significant  change  in  shot  time  was  unexpected  on  the 


UNCLASSIFIED  /27 

basis  of  DRES  experience  on  previous  trials.  It  should  be  recognized, 
however,  that  a decision  to  change  the  shot  time  is  possible  on  any 
future  trials.  In  view  of  the  DICE  THROW  experience,  the  following  sug- 
gestions are  made  to  minimize  the  risk  to  the  camera  records  in  future 
trials: 

(1)  The  most  serious  problem  would  arise  near  sunrise  or  sunset. 
Before  choosing  which  end  of  the  cylinder  is  to  be  photographed,  informa- 
tion should  be  obtained  on  the  predicted  position  of  the  sun  at  sunrise 
and  sunset  on  shot  day.  All  else  being  equal,  the  cylinder  end  to  be 
photographed  should  be  chosen  according  to  the  camera  orientation  which 
promises  least  interference  from  a rising  or  setting  sun. 

(2)  In  the  case  that  back-lighting  by  the  sun  or  low  light 
intensity  prove  to  be  a problem,  a contingency  plan  to  provide  sources 

of  artificial  light  at  each  camera  location  should  be  seriously  considered. 
Such  a light  source  must  be  capable  of  providing  light  of  sustained  in- 
tensity for  a period  of  approximately  200  milliseconds  from  the  time  of 
shock  arrival  at  any  location.  Candidates  for  consideration  would  be 
(i)  flares,  (ii)  flash  lamps,  or  (iii)  spotlights,  possibly  used  in  con- 
junction with  parabolic  or  planar  reflectors.  Integrity  of  the  light 
source  under  blast  loading  is  an  important  factor  to  be  considered.  In 
view  of  the  fact  that  the  test  area  is  cleared  of  personnel  approximately 
two  hours  before  shot  time,  another  factor  for  consideration  is  the 
ability  to  activate  the  light  source  by  remote  control  close  to  shot  time. 

2.4  Dust  Collectors 

Since  it  was  known  that  dust  entrained  in  the  blast  wave  could 
significantly  alter  the  measured  drag  pressure,  it  was  felt  to  be  impor- 
tant to  obtain  some  indication  of  the  contribution  of  dust  loading.  A 
series  of  simple  passive  dust  collectors  consisting  of  7-foot  high  verti- 
cal aluminum  channels  filled  with  grease  were  located  at  strategic  points 
on  the  layout  (Figure  14).  The  results  of  this  experiment  are  the  subject 
of  a separate  report  (Naylor,  1977). 

The  ground  surrounding  the  Canadian  projects  was  treated  with  a 
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sprayed-on  oil-based  coating  approximately  1/8  inch  thick.  Camera  re- 
cords and  dust  collectors  confirmed  that  the  coating  was  highly  effective 
in  suppressing  dust.  The  extent  of  the  treated  ground  can  be  clearly 
seen  in  Figure  9. 

Production  of  a relatively  dust-free  blast  wave  was  desirable 
in  order  to  relate  experimental  conditions  as  closely  as  possible  to 
real-life  conditions  encountered  by  Navy  materiel. 
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3.  DATA  ANALYSIS 
3.1  General 

The  goal  of  the  analysis  was  to  obtain  the  aerodynamic  drag 
coefficient  as  a function  of  time,  CD(t).  CD  is  defined  by  the  equation: 

PD(t)  = CD(t)  q(t)  (1)  j 

where  PD  is  drag  pressure 

CD  is  drag  coefficient 
q is  dynamic  pressure 

and  the  time-dependence  of  each  quantity  is  noted  explicitly. 

A power  series  in  time  was  fitted  to  the  velocity-time  data, 
v(t),  and  displacement-time  data,  x(t),  as  described  in  Sections  3.2  and 
3.3.  The  best-fit  polynomials  to  the  velocity-time  and  displacement-time 
data  are  recorded  in  Tables  13  and  14  respectively. 

Drag  pressure  is  related  to  the  slope  of  v(t)  and  curvature  of 
x(t),  through  the  relations: 


dv(t) 

= a(t), 

= a ( t ) 

(2) 

dt 

dt2 

and 

pD(t) 

-Sr-(t) 

(3) 

where  x is  cylinder  displacement 
v is  cylinder  velocity 
a is  cylinder  acceleration 
m is  cylinder  mass 
A is  frontal  area  of  cylinder 

t is  elasped  time  after  arrival  of  shock  front  at  cylinder. 

For  the  purposes  of  this  experiment,  dynamic  pressure,  q(t),  was 


replaced  in  Equation  (1)  by  the  closely-related  quantity  impact  pressure, 
q j ( t ) . The  derivation  of  qj(t)  from  measured  free-field  overpressure- time 
histories,  and  the  reason  for  replacing  q(t)  with  q j ( t ) , are  elaborated 
in  Sections  3.4  and  3.5. 
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3.2  Velocity  Transducer  Data 

3.2.1  Conversion  from  analog  to  digital  velocity-time  signal.  The 
analog  signals  recorded  on  magnetic  tape  during  the  trial  were  digi- 
tized after  the  trial  at  a digitizing  rate  of  16  KHz  using  an  analog- 
to-digital  converter.  Up  to  100  msec  of  smooth  baseline  (unaffected 
by  ground  shock)  was  retained  on  either  side  of  the  signal.  A straight 
line  was  fitted  to  the  baseline  data.  After  subtracting  the  baseline, 
and  accurately  determining  the  crossover  point  at  which  the  transducer 
signal  changed  polarity,  the  second  half  of  the  transducer  signal  was 
multiplied  by  -2  (see  Section  2.2.2).  Using  a digitized  1 volt  reference 
signal  of  approximately  1 second  duration  recorded  just  before  the 
transducer  signal,  the  calibration  values  in  ft/sec/volt  (which  had  been 
measured  in  the  laboratory  at  DRES  prior  to  the  trial)  were  applied  to 
the  digitized  signal  to  produce  a final  velocity-time  signal.  The 
starting  point  of  the  signal,  determined  by  inspection  of  the  digitized 
data,  was  taken  as  zero  time  in  the  subsequent  analysis.  The  uncertainty 
in  this  zero  position  is  0.25  msec.  The  first  3 msec  of  data  were 
omitted  to  exclude  from  the  fit  the  initial  diffraction  phase  which  can 
last  for  as  long  as  1.5  msec. 

After  establishing  by  visual  inspection  that  digital  smoothing 
of  the  transducer  signals  would  not  suppress  any  significant  features 
in  the  data,  a smoothing  was  performed  by  averaging  each  consecutive 
interval  of  8 points.  At  the  time  that  the  curve  fitting  was  performed, 
the  interval  between  data  points  was  0.5  msec. 

A method  of  checking  the  velocity  calibration  for  each  trans- 
ducer is  described  in  Appendix  D and  the  results  of  this  investigation 
are  recorded  in  Table  11. 

The  characteristic  transducer  response  time  (approximately 
2 msec)  was  not  fast  enough  to  follow  the  abrupt  change  in  velocity 
occurring  during  the  initial  diffraction  phase  of  shock  loading  on  the 
cylinder,  which  lasts  for  about  1 millisecond.  There  was,  therefore, 
little  point  to  analyzing  veloci ty-time  data  during  the  initial  re- 
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covery  time  of  the  transducer.  For  this  reason,  only  data  from  3 msec 
onward  were  retained  for  analysis.  No  analysis  of  the  diffraction  phase 
of  shock  loading  was  attempted. 

3.2.2  Philosophy  of  curve-fitting.  The  reasons  for  choosing  a power  series 
in  time  (polynomial  regression  method)  to  fit  the  velocity-time  data  are 
thoroughly  discussed  in  Appendix  E. 

A least  squares  criterion  was  applied  to  determine  the  best  fit 
to  the  velocity  data.  All  data  points  were  weighted  equally.  The  data 
were  fitted  by  a power  series  of  increasing  order  up  to  N=6  (7  coeffi- 
cients). It  was  necessary  to  have  some  criteria  for  determining  which 
order  of  polynomial  provided  the  best  fit  to  the  data. 

3.2.3  Factors  affecting  choice  of  best-fit  function. 

(a)  From  the  theory  of  least  squares  fitting  (Appendix  E,  Part 
II),  it  was  clear  that  the  addition  of  a higher  order  term  to  the  fitting 
function  provided  an  extra  degree  of  freedom  and  automatically  ensured 
that  a smaller  root  mean  square  deviation  was  obtained.  However,  if  the 
improvement  in  the  RMS  deviation  was  small,  the  uncertainties  in  the  fit- 
ting coefficients  themselves  actually  increased  because  the  effect  of 
adding  an  extra  degree  of  freedom  more  than  offset  the  small  improvement 
in  the  RMS  deviation.  In  this  case,  since  acceleration  was  a function 

of  these  coefficients,  the  uncertainty  in  acceleration  (and  hence  drag 
pressure)  increased  as  the  order  of  the  fit  increased.  To  minimize  this 
effect,  the  lowest  order  fit  which  satisfied  the  conditions  set  forth  in 
items  (b)  and  (c)  below  was  chosen  as  the  best  fit. 

(b)  The  drag  pressure-time  curve  was  derived  directly  from  the 
first  derivative  of  the  curve  fitted  to  the  velocity-time  data.  As  the 
order  of  the  fitting  polynomial  was  increased,  in  general  the  shape  of 
the  derived  drag  pressure-time  curve  changed.  However,  if  the  scatter  in 
the  velocity-time  data  was  not  too  large,  a point  was  quickly  reached  for 
which  addition  of  the  next  higher  order  term  in  the  fitting  function  did 
not  significantly  change  the  shape  of  the  derived  drag  pressure-time 
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curve  (i.e.,  fitting  procedure  had  produced  a stable  first  derivative). 
One  could  then  claim  with  some  confidence  that  adding  the  next  higher 
order  term  produced  no  significant  new  information  on  the  shape  of  the 
drag  pressure-time  curve.  In  other  words,  if  a fit  of  order  N and  a 
fit  of  order  N+l  produced  essentially  the  same  shape  for  P^( t ) , the  fit 
of  order  N was  chosen  as  the  one  which  contained  all  of  the  available 
information  on  the  shape  of  PD(t). 

(c)  While  the  overall  shape  of  the  drag  pressure-time  curve 
cannot  be  known  with  certainty,  two  characteristics  of  the  curve  are 
known,  based  on  the  shape  of  the  dynamic  pressure  curve:  (1)  initially, 
the  slope  is  negative;  (2)  at  large  times,  when  the  dynamic  pressure 
approaches  zero  asymptotically,  the  drag  pressure  approaches  zero  asymp- 
totically. 

When  one  performs  a least  squares  fit  to  the  velocity  data 
with  a power  series,  however,  there  are  no  special  constraints  to 
guarantee  that  the  drag  pressure  derived  from  the  best-fit  polynomial 
will  have  these  two  characteristics.  When  the  scatter  in  the  velocity 
data  is  small,  the  data  themselves  provide  sufficient  constraint  to 
ensure  that  the  curve  for  drag  pressure  derived  from  the  fit  will  have 
these  characteristics.  However,  if  the  scatter  in  the  velocity  data  is 
sufficiently  large,  the  derived  drag  pressure  may  fail  to  exhibit  either 
one  or  both  of  these  characteristics.  When  this  occurs,  it  is  necessary 
to  exercise  judgment  to  rule  out  such  a fit  on  the  grounds  that  it  is 
not  physically  reasonable.  This  is  particularly  true  for  cases  in  which 
the  derived  drag  pressure  drops  to  zero  and  becomes  negative  before  the 
end  of  the  fitting  region  is  reached. 

Since  one  is  relying  entirely  upon  the  data  themselves  to 
constrain  the  fitting  function  to  have  the  proper  behaviour  at  late 
times,  it  would  seem  to  be  important  for  the  transducer  time  record  to 
extend  up  to  the  point  where  the  slope  of  the  velocity-time  curve  is 
approximately  zero  (i.e.,  where  the  dynamic  pressure  is  negligible  com- 
pared to  its  initial  value).  While  this  was  the  case  for  all  of  the 
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transducer  data  collected  in  Event  DICE  THROW  (except  for  a couple  of 
truncated  signals  caused  by  broken  magnets,  etc.),  it  was  not  always  the 
case  in  previous  trials,  where  transducers  with  shorter  working  lengths 
were  used. 

In  summary,  before  a fitted  function  was  accepted  as  an  accurate 
description  of  the  variation  of  drag  pressure  with  time,  three  conditions 
had  to  be  satisfied: 

(1)  reasonable  limits  on  uncertainty  (low  order  of  polynomial), 

(2)  stable  first  derivative,  and 

(3)  correct  physical  behaviour  at  early  times  (when  dynamic 
pressure  is  large)  and  at  later  times  (when  drag  pressure 
is  decreasing  to  zero  asymptotically). 

If  the  fitting  function  failed  to  meet  these  criteria  for  a 
particular  data  set,  then  all  high  order  fits  were  rejected  as  unsuitable 
and  a linear  fit  to  drag  pressure  was  chosen  as  the  best-fit  function. 
Since  the  dynamic  pressure  curve  is  non-linear,  it  is  clear  (Equation  1) 
that  a linear  fit  cannot,  in  general,  represent  the  detailed  shape  of 
PD(t).  It  is,  however,  useful  in  showing  the  trend  of  PQ(t)  and  can  be 
considered  as  a coarse  averaging  function.  For  data  sets  which  fit  Pp(t) 
with  a linear  function,  it  should  be  noted  that  the  detailed  shape  of 
the  derived  CD(t)  curve  is  not  significant.  Only  values  averaged 
over  25  msec  intervals  were  accepted  as  meaningful  in  these  cases. 

3.2.4  Effect  of  non-random  fluctuations  in  velocity-time  data.  As  was 
the  case  in  all  previous  trials,  non-random  fluctuations  were  evident  in 
all  transducer  signals.  These  could  be  subdivided  into  two  categories: 

(1)  pure  (damped)  sinusoidal  oscillations;  (2)  irregular  fluctuations. 

(1)  Pure  Sinusoidal  Oscillations,  Large  single-frequency 
oscillations  were  observed  in  the  velocity-time  spectra  from  transducers 
attached  to  cylinders  3,  4 and  5.  By  inspection  of  the  corresponding 
camera  records,  it  was  established  that  they  were  oscillations  of  the 
solid  centre  shaft  of  the  cylinder  to  which  the  transducers  were  attached. 
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The  fact  that  the  amplitudes  of  oscillation  were  larger  in  general 
than  those  observed  in  velocity-time  traces  from  previous  trials  may 
be  attributed  to  two  factors:  (1)  the  centre  rod  in  this  trial  was 
continuous  and  of  order  5 to  8 feet  in  length  (cf.  the  short  - approxi- 
mately 1 foot  long  - end  rods  of  previous  trials);  (2)  the  absolute 
length  of  rod  protruding  beyond  the  cylinder  end  was  14  inches  for  DICE 
THROW,  compared  to  only  6 inches  for  previous  trials.  Each  of  the 
above  factors  (1)  and  (2)  would  tend  to  favour  a higher  initial  dis- 
placement for  the  end  of  the  centre  rod  relative  to  the  cylinder  axis. 

An  attempt  was  made  to  remove  the  prominent  single-frequency 
component  by  performing  a Fourier  transform  of  the  velocity-time  spec- 
trum, subtracting  the  unwanted  component,  and  reconstructing  the  filtered 
velocity-time  function  by  performing  an  inverse  Fourier  transform.  Due 
to  the  short  length  of  the  v(t)  spectrum,  the  presence  of  gaps  at  the 
beginning  and  in  the  middle  of  the  spectrum,  and  the  fact  that  the 
oscillation  was  damped,  it  proved  impossible  to  apply  a sufficiently 
precise  frequency  filter  which  would  remove  the  oscillatory  component 
without  simultaneously  distorting  the  shape  of  the  main  velocity-time 
signal . 

The  next  approach  employed  was  an  attempt  to  fit  the  velocity- 
time spectrum  with  a function  of  the  form: 

V ( t ) = V] ( t)  + v 2 ( t ) (4) 

where 

Vj(t)  = ax  + a2t  + a3t2;  v2(t)  = a4e~a5tsin(2ira6t  + a7).(5) 

This  function,  which  includes  an  explicit  sinusoidal  term,  contained 
seven  fitting  parameters  (aj-  - - - a7).  A least  squares  best-fit 
criterion  was  adopted,  and  the  best-fit  function  was  found  by  a para- 
meter search  method.  Because  the  computer  time  consumed  increases 
rapidly  with  the  number  of  parameters  using  this  method,  it  proved 
necessary  to  limit  the  polynomial  part  of  the  fitting  function  to 
second  order  in  order  to  keep  the  overall  number  of  parameters  manage- 
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able.  While  it  was  recognized  that  a second  order  polynomial  was  not 
of  sufficiently  high  order  to  describe  adequately  the  v(t)  trace,  the 
results  of  this  fit  were  still  of  interest  because  the  model  function 
attempted  to  take  into  account  the  damped  sinusoidal  oscillation.  It 
was  useful  to  compare  the  best-fit  coefficients  a^  a2,  a3  obtained 
using  this  function  with  the  best-fit  coefficients  obtained  by  linear 
least  squares  fitting  procedure  using  a second  order  power  series  only 
(described  in  Appendix  E).  The  results  of  the  two  procedures  are  sum- 
marized in  Table  12.  They  indicate  that,  at  least  for  a second  order 
polynomial  fit  to  velocity,  the  two  methods  give  similar  answers  for 
the  polynomial  describing  the  velocity-time  curve.  In  fitting  the  velo- 
city-time records  with  a polynomial  of  order  higher  than  2,  the  implicit 
assumption  is  made  that  any  oscillations  about  the  mean  velocity  are 
sufficiently  rapid  that  they  do  not  seriously  affect  the  choice  of  best- 
fit  polynomial.  The  results  in  Table  12  lend  some  credence  to  this 
assumption. 

(2)  Irregular  Fluctations.  For  cylinders  2,  6 and  7,  irregular 
fluctuations  were  superimposed  on  the  sinusoidal  oscillations.  The  fluc- 
tuations were  most  likely  caused  by  friction  between  the  moving  magnet 
and  the  surrounding  coil  housing.  This  type  of  fluctuation  is  the  most 
problematic  because  it  is  neither  completely  random  nor  periodic.  It  is 
impossible  to  assess  accurately  the  uncertainty  in  the  fitted  function 
attributable  to  the  irregular  nature  of  these  fluctuations.  No  attempt 
was  made  to  do  so,  and  in  the  error  analysis  it  was  assumed  that  the 
fluctuations  were  random.  Because  the  fluctuations  are  rapid  compared  to 
the  slope  of  the  v(t)  curve,  this  is  not  a bad  assumption.  There  is  no 
a priori  means  of  knowing  if  this  is  the  case,  however;  it  must  be  assumed. 

3.3  Camera  Data 

3.3.1  Use  of  film  reader.  Developed  films  from  the  high-speed  cameras 
were  analyzed  with  the  aid  of  a precision  film  reader.  Timing  marks  pro- 
jected onto  the  film  at  10  millisecond  intervals  were  used  to  establish 
the  framing  rate.  A horizontal  distance  scale  in  each  film  frame  was  pro- 
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vided  by  a photomarker  plate  attached  to  the  support  plate  nearest 
the  camera  and  marked  off  over  a 12  inch  interval  in  alternate  black- 
and-white  1-inch  wide  bands.  The  zero  of  coordinates  was  defined  for 
each  film  frame  to  be  the  junction  of  the  photomarker  plate  with  the 
vertical  back  edge  of  the  support  plate. 

For  those  cameras  with  50  mm  focal  length  lenses  (Table  10), 
non-linearity  across  the  field  of  view  could  be  neglected.  For  those 
cameras  with  13  mm  lenses,  a correction  had  to  be  applied.  It  would 
be  an  advantage  in  future  trials  to  use  cameras  with  focal  lengths 
longer  than  13  mm  to  minimize  the  importance  of  this  correction. 

The  measuring  position  on  the  cylinder  was  defined  by  the 
junction  of  alternate  black  and  white  sectors  painted  onto  the  end 
plates.  Because  the  end  plate  is  9 inches  farther  from  the  camera 
than  the  photomarker  plate,  a simple  geometrical  correction  had  to  be 
applied  to  the  measured  position  coordinates.  It  was  also  necessary 
to  apply  a correction  to  account  for  motion  of  the  camera  and  mounting 
post  under  blast  loading.  This  correction  was  necessitated  by  a 
shift  in  the  optic  axis  of  the  camera  with  respect  to  the  object  being 
photographed.  Since  no  cross-hairs  were  present  on  the  camera  lens 
to  define  the  optic  axis,  a reference  point  on  the  edge  of  the  film 
frame  was  used  to  keep  track  of  the  camera  motion.  The  accuracy  of 


position  measurement  was  estimated  to  be  ± 0.04  inches  before  any  cor- 
rections were  applied. 

3.3.2  Philosophy  of  curve-fitting.  The  same  considerations  which 
governed  the  fitting  of  the  velocity-time  data  discussed  in  Section  3.2 
applied  to  the  fitting  of  the  displacement-time  record  from  the  high- 
speed cameras,  except  that  one  was  interested  in  the  second,  rather 
than  first,  derivative,  and  the  record  was  continuous.  In  addition, 
the  ability  to  observe  the  cylinder  end,  rather  than  the  end  of  the 
cylinder  shaft,  meant  that  the  oscillations  of  the  cylinder  shaft,  so 
prominent  in  the  velocity  transducer  data,  were  absent  in  the  camera 
data. 
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3.4  Free-Field  Overpressure  Measurements 

Side-on  pressure  gauges  mounted  at  ground  level  were  used  to 
record  overpressure-time  histories  at  strategic  points  on  the  Canadian 
layout.  These  measurements  are  the  subject  of  a related  Suf field  Tech- 
nical Paper  (Winfield,  1977).  Four  side-on  pressure  gauges  were  located 
in  the  vicinity  of  the  20  psi  overpressure  position,  six  gauges  in  the 
vicinity  of  the  10  psi  overpressure  position,  and  four  gauges  in  the 
vicinity  of  the  7 psi  overpressure  position. 

Each  overpressure-time  curve  was  assumed  to  follow  the  empi- 
rical Fried lander  decay  formula: 

P(t)  = PqF  (6) 

* 

where 


F = 


1 - r 


-k^ 

X 


with  pQ  = peak  overpressure  (psi) 

t+  = duration  of  positive  overpressure  phase 
k = Friedlander  decay  constant  (empirically  determined). 

The  positive  duration,  t+,  was  determined  by  visual  inspection 
of  the  digitized  pressure-time  records. 

The  overpressure  impulse,  I,  defined  by: 

;♦ 

I = J p(t)  dt 


(7) 


was  obtained  by  numerical  integration  of  the  area  under  the  measured 
pressure-time  record  from  t = o to  t = t . Integration  of  Equation  6 
from  t = o to  t = t leads  to  the  equation: 


I 


~ (1  - e'k) 
k2 


(8) 


The  function  on  the  right  is  an  unique  function  of  the  decay 
constant  k only.  This  function  was  plotted  and  the  value  of  k determined 
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graphically  for  each  pressure  gauge  by  calculating  the  ratio  — y- 

using  experimental  values  of  I,  pQ  and  t+  determined  directly  from  the 
measured  pressure-time  records.  Once  the  parameters  I,  pQ,  t+  and  k 
were  determined  for  each  gauge  at  a given  nominal  peak  overpressure 
location,  a best  value  was  determined  for  each  parameter  by  averaging 
the  results  from  all  the  gauges  at  that  peak  overpressure  location. 

The  scatter  in  the  values  of  the  parameters  about  the  mean  value  was 
used  to  provide  an  estimate  of  the  uncertainty  in  each  parameter.  The 
results  of  this  procedure  are  summarized  in  Table  15. 

3.5  Impact  Pressure  Calculations 

The  dynamic  pressure  q and  impact  pressure  qj  were  assumed 
to  decay  as  F2  (Glasstone,  1957),  i.e., 


q(t)  = qQF2 

PlU)  = ql0F2 


(9) 

(10) 


where  the  peak  dynamic  pressure,  qQ,  is  determined  from  the  Rankine- 
Hugoniot  relations  at  the  shock  front  to  be  (Glasstone,  1957): 

n 2 


5 ro 
2 (P0  + 7pa5 


(11) 


and  the  peak  impact  pressure  is  determined  to  be  (Kinney,  1962): 


*Io 


qn  + Mo 

P.) 


(12) 


where 

Po  = peak  overpressure 
p,  = ambient  pressure. 

a 

It  has  been  the  practice  in  recent  years  at  our  Establish- 
ment to  define  drag  coefficient  in  terms  of  impact  pressure,  rather 
than  dynamic  pressure  (see  Equation  1),  because  drag  force  for  com- 
pressible fluids  is  directly  related  to  impact,  rather  than  dynamic, 
pressure.  This  practice  has  been  continued  in  this  part  of  the  pre- 
sent report.  The  ratios  of  impact  pressure  to  dynamic  pressure  at  the 
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20.1,  9.7  and  6.7  psi  peak  overpressure  locations  were  1.103,  1.039, 

1.022,  respectively,  based  on  Equations  9,  10,  11  and  12. 

In  calculating  impact  pressure  for  a given  peak  overpressure 
location,  the  average  values  of  the  parameters  from  all  of  the  pressure 
gauges  at  that  overpressure  location,  as  listed  in  Table  15,  were  used. 

To  obtain  an  estimate  of  the  uncertainty  in  impact  pressure  at  that  loca- 
tion, impact  pressure-time  curves  were  calculated  using  the  values  of  pQ, 
t+  and  k determined  for  each  gauge.  At  10  msec  intervals,  the  RMS  deviation 
of  these  curves  from  the  "average"  curve  was  calculated.  This  RMS  deviation, 
expressed  as  a percentage,  is  plotted  as  a function  of  time  for  each  peak 
overpressure  location  in  Figure  32. 

The  impact  pressure  curves  used  in  the  data  analysis  are  plotted 
in  Figure  31 . 

3.6  Calculation  of  Mach  and  Reynolds  Numbers 

Free  stream  Mach  and  Reynolds  numbers  were  calculated  using 
standard  definitions.  Fluid  velocity  was  assumed  to  decay  as: 


where 


u = u F 
o 


F is  defined  in  Equation  6, 


and  uQ  was  derived  in  terms  of  pQ  and  pa  from  the  Rankine-Hugoniot 
relations  across  the  shock  front. 

The  temperature  of  the  flow  behind  the  shock  front  was  approxi- 
mated by  the  isentropic  relation.  The  kinematic  velocity  was  described 
by  a power  series  in  temperature,  where  the  coefficients  of  the  series  were 
obtained  by  fitting  a power  series  to  values  of  kinematic  viscosity  for  air 
at  specific  temperatures. 

Details  of  these  calculations  are  given  in  a previous  report 
( Mel  1 sen,  1974). 
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4.  RESULTS 

4. 1 General 

A qualitative  summary  of  results  is  presented  in  Table  27. 
Cylinder  1 failed  to  undergo  free  flight  and  suffered  severe  damage  due 
to  a blast  anomaly  which  moved  up  the  east  side  of  the  Canadian  sector. 
The  anomaly  is  discussed  in  Section  4.10.  Post-shot  views  of  the 
cylinder  and  the  east  and  west  support  plates  are  presented  in  Figures 
10,  11  and  12  respectively.  Figure  13  presents  a post-shot  view  of 
cylinder  5,  showing  the  lateral  displacement  caused  by  a lateral  pres- 
sure component  associated  with  the  blast  anomaly.  Figure  14  is  a 
post-shot  view  of  one  of  the  7-foot  high  greasy  stake  dust  collectors. 

The  height  of  the  dust  cloud  is  indicated  by  the  change  in  reflectivity 
of  the  greased  surface  at  a height  of  3 to  3.5  feet. 

Figures  15  through  25  present  results  of  fits  to  velocity 
transducer  data  and  camera  data  for  cylinders  2 through  7 (except  for 
cylinder  2 where  only  transducer  data  are  available).  Figures  26 
through  30  and  Tables  16  through  26  present  derived  drag  pressures  and 
drag  coefficients  for  cylinders  3 through  7 from  both  camera  and  trans- 
ducer data.  A qualitative  summary  of  results  for  each  cylinder,  with 
qualifying  remarks,  is  presented  in  Table  27.  A summary  of  best  values 
for  the  measured  drag  coefficients  is  presented  in  Table  28. 

4.2  Velocity  Transducer  Data 

The  data  from  east  and  west  ends  on  a given  cylinder  were 
analyzed  separately.  The  best-fit  curve  is  drawn  through  the  velocity- 
time data  as  a solid  line.  The  dotted  lines  in  the  figures  represent 
± one  root  mean  square  (RMS)  deviation  in  the  scatter  of  data  about  the 
best-fit  curve.  It  was  assumed  for  purposes  of  error  analysis,  and  in 
the  absence  of  better  information,  that  the  RMS  deviation  had  a constant 
value  at  all  points  on  the  curve.  The  acceleration  curve  is  the  first 
derivative  of  the  best-fit  velocity  function.  The  dotted  lines  on  the 
acceleration-,  drag  pressure-,  and  drag  coefficient-time  curves,  however. 
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represent  ± three  standard  deviations  (99%  confidence  interval). 

Power  series  with  terms  up  to  the  fourth  power  in  time  (5  para- 
meters) were  fitted  to  the  velocity  data  for  cylinders  3,  4,  6 and  7. 
Velocity  data  for  cylinders  2 and  5 were  fitted  by  power  series  with 
terms  up  to  the  second  power  in  time  (3  parameters). 

4.3  High-Speed  Camera  Data 

Cameras  recorded  the  motion  of  the  west  cylinder  ends  only.  In 
the  figures,  the  best-fit  power  series  curve  is  drawn  as  a solid  line 
through  the  displacement-time  data.  Dotted  lines  representing  ± one  RMS 
deviation  are  also  drawn,  but  are  not  evident  on  most  drawings  because  the 
deviation  is  so  small.  On  the  velocity-,  acceleration-,  drag  pressure-, 
and  drag  coefficient-time  curves,  however,  the  dotted  lines  represent  ± 
three  RMS  deviations  about  the  best-fit  curve  (99%  confidence  interval). 
Displacement  data  for  cylinders  3,  4,  5 and  6 were  fit  by  a power  series 
with  terms  up  to  the  fifth  power  in  time  (6  parameters).  Displacement  data 
for  cylinder  7 were  fit  by  a power  series  with  terms  up  to  the  third  power 
in  time  (4  parameters). 

4.4  Drag  Coefficient  vs  Reynolds  Number  - Before  Corrections 

(Figure  36) 

For  a given  record,  the  drag  pressure-time  curve  was  divided  by 
the  appropriate  impact  pressure-time  curve  (Figure  31)  to  obtain  drag  coef- 
ficient as  a function  of  time,  Cp(t).  The  CQ(t)  curves  are  plotted  in 
Figures  26  through  30  for  each  cylinder.  Dotted  lines  represent  ± three 
standard  deviations  of  uncertainty. 

In  Figure  36,  the  resultant  Cp(t)  curves  are  plotted  as  a function 
of  Reynolds  number  for  all  cylinders  for  which  the  free  stream  Mach  number 
is  less  than  the  critical  value  Mc  = 0.48.  For  this  Mach  number  range,  the 
dependence  of  CD  upon  Mach  number  is  slight. 

In  Figure  36,  data  which  required  a linear  fit  to  drag  pressure 
have  not  been  included  because  a linear  fit  was  felt  to  provide  no  de- 
tailed information  on  the  shape  of  the  drag  pressure-time  curve  (see  Sec- 
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tion  3.2).  The  Cp(t)  values  are 
three  reasons: 


presented  as  bands  of  uncertainty  for 


(1)  To  permit  a visual  comparison  of  the  relative  accuracies 
of  the  velocity  transducer  and  high-speed  camera  techniques. 

(2)  To  emphasize  that,  for  a given  data  record,  the  uncertainty 
in  the  derived  drag  coefficients  is  not  constant  across  the  record.  The 
uncertainty  is  least  near  the  middle  of  each  record. 

(3)  To  show  the  measure  of  agreement  between  drag  coefficients 
obtained  using  the  velocity  transducer  and  those  obtained  using  the  high- 
speed camera  technique. 

In  this  figure,  and  in  the  plots  of  Cp(t)  in  Figure  26  through 
30,  no  uncertainty  in  the  calculated  impact  pressure  has  been  included. 
This  has  been  done  so  that  the  ratio  of  drag  coefficients  with  and  without 
extended  end  plates  could  be  formed  directly  to  assess  the  importance  of 
end  effects.  In  such  a ratio  impact  pressure  cancels  out,  so  the  un- 
certainties in  Figures  26  through  30  and  Figure  36  are  the  appropriate 
ones  to  use  for  assessing  end  effects. 

Included  for  completeness  in  Figure  36  is  a solid  curve  re- 
presenting drag  coefficients  measured  in  a wind  tunnel  under  steady- 
state  flow  conditions  (Delaney  and  Sorensen,  1953;  Gowen  and  Perkins, 
1953).  The  extension  of  the  steady-state  results  to  higher  Reynolds 
numbers  (Roshko,  1961)  is  represented  by  the  dotted  portion  of  the  curve. 


I ; 


4.5  End  Effects  (Figure  37) 

The  flow  of  air  over  the  ends  of  finite-length  cylinders  can 
produce  a measured  drag  coefficient  which  is  different  from  the  value 
that  would  be  measured  for  an  infinitely-long  cylinder.  Since  the  main 
structural  members  of  the  Canadian  Navy  masts  and  antennae  tested  in 
DICE  THROW  either  had  relatively  large  length/diameter  (L/D)  ratios  or 
were  attached  at  either  end  to  other  members,  drag  coefficients  for 
infinite-length  cylinders  were  the  appropriate  input  to  structural  res- 
ponse calculations  for  these  structures. 
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In  the  present  experiment,  thin  extended  end  plates  were  attached 
to  the  ends  of  some  cylinders  (Figures  4 and  5)  to  prevent  air  flow  around 
the  ends  of  the  cylinders,  thereby  eliminating  end  effects. 

There  was,  of  course,  a contribution  to  the  overall  drag  on  the 
cylinder  due  to  drag  on  the  end  plates  themselves.  However,  because  the 
end  plates  had  bevelled  knife  edges,  and  the  air  flow  was  expected  to  be 
nearly  parallel  to  the  faces  of  the  end  plates,  the  main  contribution  to 
end  plate  drag  was  skin  friction  drag,  for  which  the  maximum  drag  coef- 
ficient, according  to  Hoerner  (1958),  is  0.008.  Using  this  value  for  Cp, 
the  fractional  contribution  by  the  end  plates  to  the  overall  measured 
cylinder  drag  pressure  was  assumed  to  be  given  by: 


^CD^end  plate 


^D^end  plate  + ^CD^cyl inder 


where 


''cylinder 
^D^end  plate 

(r  ) 

' D cylinder 


^end  plate  t^ie  exPosec^  surface  area  of  the  two 

end  plates 

A_.-i.j_,! is  the  frontal  area  of  the  cylinder 

= 0.008 

= measured  value  from  experiment. 

These  calculated  contributions  from  the  end  plates  to  the  measured 
drag  coefficient  (approximately  10%  for  cylinder  7,  7%  for  cylinder  5,  and 
less  than  1%  for  cylinder  2)  were  then  subtracted  from  the  measured  coef- 
ficients to  produce  a set  of  corrected  coefficients  appropriate  to  in- 
finite-length cylinders.  It  is  these  corrected  coefficients  which  are 
plotted  in  Figure  38. 


Since  the  experiment  included  pairs  of  identical  cylinders  with 
and  without  extended  end  plates,  at  the  same  peak  overpressure  locations,  it 
was  possible  to  measure  end  effects  directly  by  forming  the  ratio 

Cn  (with  extended  end  plates,  corrected  for  end  plate  drag) 

U • 

Cp  (without  extended  end  plates) 
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These  ratios  were  formed,  using  drag  coefficients  averaged 
over  25  msec  intervals,  for  cylinders  6 and  7 (velocity  transducer 
data)  and  cylinders  4 and  5 (camera  data).  The  results  are  plotted 
in  Figure  34  as  a function  of  elapsed  time  after  shock  arrival. 

For  cylinders  6 and  7,  with  an  L/D  ratio  of  17,  the  first 
two  points  in  Figure  37  covering  50  msec  of  motion  correspond  to  Rey- 
nolds numbers  (9-5)  x 10s,  i.e.,  just  above  the  critical  Reynolds  number 
range  (5-3)  x 105.  For  these  points,  the  CQ  ratio  lies  below  1.  A 
slight  tendency  for  the  ratio  to  decrease  with  time  after  shock  arrival 
is  noted.  All  points  are  consistent  with  the  weighted  average  value  of 
0.78. 

For  cylinders  4 and  5,  which  have  an  L/D  ratio  of  5,  all  of 
the  data  points  correspond  to  Reynolds  numbers  in  the  supercritical 
range  (16-9)  x 105.  Initial  values  of  the  Cp  ratio  are  substantially 
greater  than  1 and  are  not  consistent  with  unity  within  error.  More- 
over, the  ratio  increases  markedly  for  later  times.  The  average  value 
for  the  first  50  msec  of  motion  is  1.43  while  for  the  second  50  msec 
of  motion  it  is  1.95.  The  average  value  for  the  first  100  msec  of 
motion  is  1.58,  but  not  all  of  the  data  points  agree  with  this  value 
within  error. 

Due  to  the  failure  of  cylinder  1,  the  18- inch  diameter 
cylinder  with  extended  end  plates,  it  was  not  possible  to  measure 
end  effects  directly  using  cylinders  1 and  3.  In  view  of  the  large 
end  effects  observed  for  cylinders  4 and  5 in  the  supercritical 
Reynolds  number  range  with  an  L/D  ratio  of  5,  it  was  felt  that  sub- 
stantial end  effects  could  also  be  expected  for  the  18-inch  diameter 
cylinder  3,  whose  motion  spanned  a somewhat  higher  Reynolds  number  range 
and  which  also  had  an  L/D  ratio  of  5.  Since  no  direct  information 
was  available  for  cylinder  3,  the  average  end  effect  factor  of  1.58 
measured  for  cylinders  4 and  5 was  applied  to  the  camera  data  for 
cylinder  3.  These  corrected  data  are  plotted  in  Figure  38.  It  is 
notable  in  this  Figure  that,  even  after  substantial  end  effect  cor- 
rections, the  two  data  points  at  highest  Reynolds  number  (first  50  msec 
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of  motion)  for  the  18-inch  diameter  cylinder  lie  well  below  steady-state 
values.  To  obtain  agreement  with  the  steady-state  values,  the  measured 
drag  coefficients  (before  end  effect  corrections)  would  have  to  be  mul- 
tiplied by  an  approximate  factor  of  2.  If,  instead  of  an  average  end 
effect  factor,  one  employed  the  measured  values  for  each  25  msec  interval 
recorded  in  Figure  37,  the  two  points  at  highest  Reynolds  numbers  in 
Figure  38  would  be  depressed  a further  10%,  while  the  third  point  would 
be  elevated  a further  20%  to  lie  above  the  steady-state  value. 

4.6  Dust  Loading 

Both  the  greasy  dust  collectors  described  in  Section  2.4  and  the 
high-speed  camera  records  provided  qualitative  information  on  the  amount 
of  dust  entrained  in  the  blast  wave  during  the  cylinder  motion.  Both  dust 
collectors  (see  7-foot  tall  dust  collector  in  Figure  14)  and  camera  records 
confirmed  that  a significant  dust  cloud  existed  only  to  a height  of  about 
three  feet  above  ground,  and  that  relatively  little  dust  existed  at  the 
initial  cylinder  height  5 to  6 feet  above  ground.  One  would  expect  any 
dust  loading  to  increase  the  effective  drag  force  on  the  cylinder,  thereby 
increasing  the  measured  drag  coefficient.  The  insignificance  of  dust 
loading  is  supported  by  the  fact  that  the  measured  drag  coefficients  in 
Figure  38  were  consistent  with  steady-state  values  over  most  of  the  range 
of  measurement.  The  oil -based  coating  sprayed  onto  the  ground  in  the 
Canadian  sector  proved  highly  effective  in  suppressing  dust,  as  evidenced 
by  the  relatively  low  dust  levels  in  this  trial  compared  to  previous  trials 
in  which  dust  clouds  were  observed  to  rise  to  a height  of  7 to  8 feet 
(Leech  et  al. , 1973). 

4.7  Drag  Coefficient  vs  Reynolds  Number  - After  Corrections 

(Figure  38) 

Figure  38  is  a composite  semi-log  plot  showing  measured  drag  coef- 
ficients for  "infinitely  long"  smooth  cylinders  in  unsteady  flow  conditions 
for  Reynolds  numbers  from  (3  to  40)  x 105.  For  the  3.5-inch  and  9.5-inch 
diameter  cylinders,  data  from  cylinders  with  extended  end  plates  were  used, 
after  subtracting  a correction  for  end  plate  drag.  For  the  18-inch  dia- 
meter cylinder  without  extended  end  plates,  the  average  end  effect  factor 
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of  1.58  measured  for  the  9. 5-inch  diameter  cylinder  was  applied  to  the 
results  to  convert  them  to  values  appropriate  to  a cylinder  of  infinite 
length.  The  error  bars  on  the  data  points  in  this  Figure  include  the 
uncertainties  in  impact  pressure  plotted  in  Figure  31.  As  in  Figure 
36,  the  solid  line  represents  results  from  wind-tunnel  experiments 
in  steady-state  flow  conditions  (see  Section  4.4).  The  extension 
by  other  workers  of  these  results  to  higher  Reynolds  numbers  is  re- 
presented by  the  dashed  portion  of  the  curve. 

4.8  Surface  Roughness 

All  of  the  cylinders  were  sanded  and  polished,  after  deep 
scratches  were  filled  with  body-filler  compound,  to  ensure  that  all 
surface  imperfections  were  less  than  1/1000  of  the  cylinder  diameter 
and  that  all  scratches  present  were  in  the  direction  of  air  flow  over 
the  cylinder.  Under  these  conditions,  according  to  Hoerner,  the 
cylinder  could  be  considered  aerodynamically  smooth  and  the  effect  of 
surface  imperfections  on  the  air  flow  would  be  negligible. 

4.9  Drag  Coefficients  for  Supercritical  Mach  Numbers  - Cylinder  2 

For  cylinder  2 at  20.1  psi  peak  overpressure,  the  flow  Mach 

number  fell  from  0.63  to  0.50  during  the  first  20  msec  of  motion,  and 

0.50  to  0.38  during  the  next  20  msec  of  motion.  During  these  intervals, 

the  Mach  number  was  above  M . . . = 0.48,  so  a large  drag  coefficient 

cnti  cd  i 

was  expected.  For  cylinder  2,  only  one  transducer  record  provided  useful 
data  and  these  data  contained  large  irregular  fluctuations  on  the  main 
signal.  It  was  necessary  to  accept  a linear  fit  to  drag  pressure  to 
obtain  reasonable  uncertainty  limits;  only  average  drag  coefficients 
over  20  msec  intervals  were  accepted  as  meaningful  (see  Section  3.2.2). 
The  large  scatter  in  results  from  pressure  gauges  adjacent  to  the  cylin- 
ders at  the  20  psi  peak  overpressure  location  caused  a correspondingly 
large  uncertainty  in  impact  pressure,  an  uncertainty  which  reached  100% 
after  only  60  msec  of  cylinder  motion  (Figure  32).  The  net  result  was 
that  only  average  drag  coefficients  for  0-40  msec  were  obtained,  and 
these  had  large  uncertainties  associated  with  them.  The  results  were: 
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For  0-20  msec  CDavge  = 0.76  ± 0.20. 

For  20-40  msec  CDavge  = 0.98  ± 0.21. 

4.10  Blast  Anomaly 

High-speed  camera  records  showed  that  the  large  end  plates  on 
cylinder  1 distorted  and  separated  from  the  main  body  of  the  cylinder 
shortly  after  the  cylinder  left  the  support  stand.  Available  evidence 
suggests  that  a blast  anomaly,  in  the  form  of  a surface  precursor  jet 
moving  up  the  east  side  of  the  Canadian  sector,  was  responsible  for  the 
failure  of  cylinder  1.  This  anomaly  produced  a secondary  pressure  wave 
which  moved  diagonally  from  east  to  west  across  the  layout  behind  the 
main  shock  front.  The  dust-raising  precursor  jet  could  be  clearly  seen 
on  overhead  photographs  of  the  charge  just  after  detonation.  The  evidence 
for  the  laterally-moving  pressure  wave  follows: 

(1)  Small  secondary  pressure  peaks  were  observed  on  pressure 
records  at  the  20  and  10  psi  peak  overpressure  locations.  Correlation 

of  the  time  of  arrival  of  these  secondary  pulses  with  the  gauge  positions 
indicated  that  the  pressure  wave  responsible  was  moving  diagonally  from 
east  to  west. 

(2)  All  cylinders  which  translated  laterally  did  so  from  east 
to  west  (Figure  13) . 

(3)  The  west  support  stand  for  cylinder  1 had  been  twisted  to- 
ward Ground  Zero  and  the  stand  had  been  collapsed  (Figure  11).  The  east 
support  stand  was  somewhat  distorted  but  still  upright  (Figure  12).  The 
only  explanation  consistent  with  these  and  other  minor  pieces  of  evidence 
is  that  the  cylinder  or  cylinder  end  plates  delivered  a series  of  rapid 
blows  to  each  support  plate.  The  fact  that  the  west  support  plate  col- 
lapsed first  suggested  that  it  had  received  the  first  major  blow  from  the 
cylinder.  This  conclusion  in  turn  implied  that  the  cylinder  initially  had 
to  translate  laterally  from  east  to  west.  An  east-west  force  would  have 
been  required  to  produce  this  motion. 
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5.  DISCUSSION  OF  RESULTS 

5.1  General 

The  principal  information  gained  from  the  present  experiment 
is  (1)  the  variation  of  drag  coefficient  in  the  range  of  critical  and 
supercritical  Reynolds  numbers  under  unsteady  flow  conditions  of  a 
free-field  blast  wave,  and  (2)  a measurement  of  end  effects  under  these 
same  conditions  for  cylinders  with  length/diameter  ratios  of  5 and  17. 

All  of  the  measurements  recorded  in  Figures  36,  37  and  38 
are  for  Mach  number  less  than  the  critical  value  (Mc  = 0.48).  For 
M<Mc,  Cp  is  primarily  a function  of  Reynolds  number.  For  M>Mc,  C^  de- 
pends mainly  upon  Mach  number.  As  M increases  through  Mc,  an  abrupt 
rise  in  CD  from  0.3  to  approximately  1.2  is  observed.  This  is  a result 
of  the  fact  that,  for  M = M , flow  becomes  supersonic  at  some  point  on 
the  cylinder.  The  local  shock  wave  which  forms  causes  a buildup  in 
thickness  of  the  boundary  layer  and  a rapid  movement  of  the  separation 
point  forward  on  the  cylinder  with  an  attendant  rapid  rise  in  drag 
coefficient. 

The  principal  difference  between  measurements  made  in  steady 
and  unsteady  flow  arises  from  the  fact  that  the  unsteady  flow  is  pre- 
ceded by  a shock  front  which  diffracts  over  the  cylinder,  sending  re- 
flections back  and  forth  several  times  across  the  cylinder.  It  is  possible 
that  passage  of  the  shock  front  may  "condition"  the  following  air  blast 
flow  to  produce  drag  coefficients  which  are  different  from  those  one  would 
measure  in  the  steady-flow  conditions  encountered  in  wind  tunnel  tests. 

5.2  Variation  of  Drag  Coefficient  with  Reynolds  Number  (Figure  38) 

5.2.1  Cylinder  with  3.5-inch  diameter.  The  points  for  Reynolds  number 
in  the  range  5 x 105  to  8 x 105  are  in  good  agreement  with  steady-state 
values.  The  points  between  3 x 105  and  5 x 105,  in  the  critical  range, 
fall  well  below  steady-state  values.  This  result  might  be  attributable 
to  surface  roughness,  which  tends  to  move  the  critical  Reynolds  region 
toward  lower  Reynolds  numbers  (Hoerner,  1958).  However,  as  discussed 
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in  Section  4.8,  precautions  were  taken  to  ensure  that  the  cylinder  sur- 
face was  aerodynamical ly  smooth,  so  this  explanation  is  an  unlikely  one. 
It  is  more  probable  that  the  lack  of  agreement  is  caused  by  the  simple 
fact  that  the  low  order  power  series  used  to  fit  the  data  is  not  capable 
of  responding  to  the  rapid  change  in  Cp  which  occurs  in  this  range  of 
Reynolds  number. 


5.2.2  Cylinder  with  9.5-inch  diameter.  Drag  coefficients  derived  for 
the  9.5-inch  diameter  cylinder  for  Reynolds  numbers  in  the  range  9 x 105 
to  17  x 105  are  consistent  with  the  steady-state  values  within  error,  but 
tend  to  lie  somewhat  higher  on  average. 

5.2.3  Cylinder  with  18.0-inch  diameter.  The  two  C^  values  spanning  the 
first  50  msec  of  motion  (R  = 43  x 10s  to  29  x 1 0 5 ) lie  well  below  steady- 
state  values.  The  Cp  value  for  the  50-75  msec  interval  (R  = 29  x 105  to 
24  x 105)  is  slightly  larger  than  the  steady-state  value,  but  consistent 
with  it  within  error. 

It  is  possible  that  the  discrepancy  between  steady  and  unsteady 
C0  values  observed  at  highest  Reynolds  numbers  is  due  to  an  inadequate  end 
effect  correction  over  this  range  of  Reynolds  number.  However,  to  obtain 
agreement  with  steady-state  values  for  all  three  points,  it  would  be  neces 
sary  to  apply  an  end  effect  correction  which  decreased  with  time  after 
shock  arrival.  This  is  contrary  to  the  observed  end  effect  variation  for 
the  9.5-inch  diameter  cylinder. 

If  one  accepts  the  data  as  presented  in  Figure  38,  they  suggest 
that,  in  the  early  stage  of  unsteady  flow  for  Reynolds  numbers  of  order 
40  x 105,  CQ  is  lower  than  the  steady-state  value.  As  time  progresses, 
the  drag  coefficient  increases  to  a value  somewhat  higher  than,  but  con- 
sistent with,  the  steady-state  value.  The  mechanism  responsible  for  the 
increase  in  drag  coefficient  for  R> 1 06  is  not  completely  understood,  but 
Roshko  (1961)  has  pointed  out  the  strong  similarity  in  shape  of  the  CQ 
vs  R and  1/S  vs  R curves,  where  S is  Strouhal  number  (S  = (fd)/u  where  d 
is  cylinder  diameter,  u is  free-stream  velocity,  and  f is  the  frequency 
of  vortex  shedding  at  the  rear  of  the  cylinder).  This  similarity  suggests 
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that  drag  coefficient  is  related  to  the  frequency  of  vortex  shedding. 

If  the  initial  shock  front  conditioned  the  following  flow  pattern 
in  such  a way  as  to  increase  artificially  the  frequency  of  vortex 
shedding,  it  is  likely  that  a decreased  drag  coefficient  would  result. 

One  might  then  expect  drag  coefficient  to  increase  as  quasi-steady 
flow  developed,  in  agreement  with  observations. 

5.3  Comparison  with  Results  of  Other  Workers 

A limited  number  of  drag  coefficient  measurements  in  unsteady 
flow  are  available.  These  have  been  carried  out  primarily  at  AWRE  (U.K.) 
in  shock  tubes  and  by  DRES  in  previous  free-field  blast  trials.  In 
several  instances,  drag  coefficients  well  in  excess  of  steady-state 
values  were  observed.  In  the  case  of  past  DRES  results,  some  of  this 
difference  may  be  attributable  to  dust  entrained  in  the  free-field  blast 
wave.  Dust  loading  seems  to  have  been  a more  serious  problem  in  pre- 
vious trials  than  in  DICE  THROW  (see  Section  4.6).  Some  work  is  under- 
way at  DRES  to  examine  the  problem  of  dust  loading  on  circular  cylin- 
ders using  a mathematical  model  in  order  to  provide  some  theoretical 
limits  on  the  potential  seriousness  of  the  problem  for  some  represen- 
tative field  conditions. 

In  the  case  of  at  least  one  set  of  results  from  AWRE  (Bishop 
and  Rowe,  1967),  the  high  measured  drag  coefficient  of  0.67  for  M<Mc 
may  be  attributable  to  the  fact  that,  early  in  the  flow  history,  the 
flow  Mach  number  M was  >M  . The  authors  suggest  that  the  drag  coef- 
ficient measured  for  M<Mc  may  depend  upon  "conditioning"  of  the  flow 
while  M>Mc>  This  contention  that,  in  unsteady  flow  conditions,  the 
measured  drag  coefficient  may  depend  upon  the  history  of  the  flow  is 
carried  forward  in  other  work  at  AWRE  by  Martin,  Mead  and  Uppard  (1967). 
In  this  work  the  authors  show,  from  shadowgraph  records  during  the 
shock  diffraction  phase,  that  the  boundary  layer  separation  point  has 
been  moved  well  forward  on  the  cylinder,  and  they  argue  that,  because 
it  is  unlikely  to  re-attach  downstream  during  the  subsequent  flow,  an 
abnormally  high  drag  coefficient  is  expected  (in  agreement  with  obser- 
vations). 
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The  present  data  are  notable  because  they  are  consistent  with 
steady-state  values  over  a wide  range  of  Reynolds  number.  The  data  for 
the  18-inch  diameter  cylinder  are  new.  Until  now,  no  known  unsteady 
flow  measurements  existed  at  such  high  Reynolds  number  and  low  Mach 
number  (M<Mc  at  all  times).  Dryden  and  Hill  (1930)  measured  CQ  for  a 
12-foot  diameter,  120  foot  long  smoke  stack  (L/D  = 10)  in  a natural  wind 
of  about  25-40  mph,  which  corresponds  to  Reynolds  numbers  of  30  x 105  to 
50  x 105.  These  measurements  were,  however,  for  extremely  low  Mach  num- 
bers, and  were  not  made  in  a decaying  blast  wave  which  was  preceded  by  a 
shock  front. 

5.4  End  Effects  (Figure  37) 

It  has  been  shown  from  shadowgraph  records  (Martin,  Mead  and 
Uppard,  1967)  that  it  can  take  as  long  as  10  msec  for  quasi-steady  flow 
to  develop  over  the  cylinders  after  passage  of  the  shock  wave.  It  is 
somewhat  surprising,  however,  to  find  that  the  ratio  of  Cp's  for  in- 
finite/finite length  cylinders  is  strongly  varying  as  late  as  75  msec 
into  the  motion  (Figure  37).  Before  any  conclusions  can  be  drawn,  it 
will  be  necessary  to  re-analyze  the  data  to  ensure  that  the  observed 
strong  variation  in  CQ  ratio  for  the  9.5-inch  diameter  cylinders  is  not 
simply  an  artifact  of  the  data  analysis.  The  average  value  of  the  ratio 
over  100  msec  of  motion,  1.58,  is  quite  close  to  the  value  of  1.67 
measured  by  Dryden  and  Hill  (see  Section  5.3)  for  very  low  Mach  numbers 
and  R in  the  range  30  x 105  to  50  x 105.  The  Reynolds  number  range  covered 
in  Dryden  and  Hill's  experiment  brackets  the  range  covered  by  our  18-inch 
diameter  cylinder  (R  = 24  x 105  to  43  x 105).  Their  measured  end  effect 
was  1.67.  The  fact  that  this  value  is  close  to  our  average  end  effect 
factor  of  1.58  (measured  on  the  9.5-inch  diameter  cylinder)  lends  support 
to  the  decision  to  apply  the  factor  of  1.58  to  the  results  for  the  18-inch 
diameter  cylinder. 

The  results  for  the  3.5-inch  diameter  cylinder  are  consistent 
with  a Cp  ratio  of  unity  for  Reynolds  numbers  just  above  the  critical 
region.  This  ratio  appears  to  drop  below  unity  by  as  much  as  25%  as  the 
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critical  Reynolds  region  is  entered.  The  indicated  decrease  in  Cp 
ratio  may  not  be  significant,  however,  because  it  is  likely  that  the 
low-order  power  series  fit  to  the  data  is  inadequate  for  modelling  the 
rapid  change  in  Cp  which  probably  occurs  in  this  region. 

The  number  of  end  effect  measurements  in  unsteady  flow  con- 
ditions is  limited.  A comparison  with  present  results  is  made  dif- 
ficult, if  not  impossible,  by  the  fact  that  some  of  the  previous 
results  are  contradictory,  and  in  most  cases  no  estimates  of  the  un- 
certainties in  the  quoted  values  have  been  reported.  For  these  reasons, 
no  detailed  comparison  of  end  effect  results  has  been  attempted. 

5.5  Results  for  Supercritical  Mach  Numbers  - Cylinder  2 

The  initial  Cp  value  of  0.76  ± 0.20  measured  for  cylinder  2 
is  more  consistent  with  the  steady-state  results  for  a finite-length 
cylinder,  with  L/D  ratio  of  17  (Cp  = 0.9),  than  for  an  infinite-length 
cylinder  (Cp  = 1.3,  Gowen  and  Perkins,  1953).  Examination  of  the 
velocity-time  curve  in  Figure  15  indicates  that  the  initial  acceleration 
values  obtained  from  the  fitted  curve  may  well  be  too  low.  A more 
reliable  determination  of  acceleration  is  not  possible,  however,  given 
the  large  fluctuations  which  are  present  in  the  velocity-time  data. 
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6.  CONCLUSIONS 

1.  Aerodynamic  drag  coefficients  measured  for  "infinite- 
length"  aerodynamical ly  smooth  cylinders  under  unsteady  flow  conditions 
in  a long-duration  (160-300  msec)  free-field  blast  wave  were  in  generally 
good  agreement  with  steady-state  values  for  Reynolds  numbers  in  the 
range  5 x 105  to  16  x 105  and  Mach  numbers  <0.48.  In  the  critical  Rey- 
nolds number  range  3 x 105  to  5 x 105,  the  measured  drag  coefficients 
lay  well  below  the  steady-state  values,  but  this  was  felt  to  be  due  to 
the  inability  of  the  power  series  fitting  function  to  respond  to  the  very 
rapid  changes  in  drag  coefficient  expected  in  this  region.  In  the  Rey- 
nolds number  range  30  x 105  to  40  x 105,  measured  unsteady  flow  drag  coef- 
ficients were  approximately  30%  lower  than  steady-state  values.  Further 
experiments  would  be  necessary  to  establish  whether  this  difference  is 
due  to  an  inadequate  correction  for  end  effects  or  to  a real  physical 
effect  associated  with  the  diffraction  of  the  shock  front  across  the 
cylinder. 

2.  Measurement  of  drag  coefficients  for  identical  cylinders 
with  and  without  extended  end  plates  permitted  the  direct  measurement  of 
end  effects  for  finite-length  cylinders  by  forming  the  ratio  Cg(infinite)/ 
Cp(finite).  For  the  cylinders  with  a length/diameter  (L/D)  ratio  of  5, 

an  average  ratio  of  1.6  was  observed  over  a Reynolds  number  range  of 
9 x 105  to  16  x 105  (M<0.48).  For  the  cylinder  with  L/D  of  17,  an  average 
CD  ratio  of  0.8  was  observed  for  Reynolds  numbers  in  the  range  of  3 x 105 
to  8 x 105.  The  ratio  was  observed  to  decrease  as  Reynolds  number  dropped 
from  the  supercritical  to  critical  range.  Further  data  analysis  and 
experimentation  are  required  to  confirm  the  strong  increase  in  end  effect 
ratio  with  time  after  shock  arrival  which  was  observed  for  the  cylinder 
with  an  L/D  ratio  of  5. 

3.  Greasy  stake  dust  collectors  and  high-speed  camera  records 
confirmed  that  the  oil -based  coating  sprayed  onto  the  ground  in  the  Cana- 
dian sector  proved  highly  effective  in  suppressing  dust.  It  is  probable 
that  the  dust  loading  on  the  cylinders  was  negligible  during  the  first 
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100-150  msec  of  motion  over  which  measurements  were  taken. 

4.  Cylinder  1 at  the  20.1  psi  peak  overpressure  location 
failed  due  to  the  influence  of  a ground  precursor  type  of  blast  anomaly 
which  moved  up  the  east  side  of  the  Canadian  sector  and  produced  a 
secondary  pressure  wave  travelling  diagonally  from  east  to  west  across 
the  Canadian  layout. 
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APPENDIX  A 

STATEMENT  OF  REQUIREMENT 

The  original  Canadian  Forces  request  for  further  drag  cylinder 
measurements  in  Event  DICE  THROW  came  from  the  Director  of  Maritime 
Facilities  and  Resources.  The  essence  of  the  request  is  reproduced 
below  (points  3 and  4 of  the  above  reference): 

3.  "It  is  recommended  that  drag  cylinders  with  approximately 
3 inch,  18  inch  and  24  inch  diameters  be  tested  each  at  two  blast 
overpressures.  The  blast  overpressures  will  be  in  the  order  of 

10  psi  but  will  become  finalized  through  model  mast  considerations." 

4.  "The  18  inch  and  24  inch  diameter  circular  cylinders  will 
provide  an  almost  complete  spectrum  of  drag  loadings  on  cylinders. 

In  addition,  these  cylinders  should  provide  some  guidelines  for 
future  designs  of  large  pole-type  masts.  The  3 inch  diameter  cylin- 
ders will  provide  experimental  results  for  useful  correlation  to 
actual  members  of  the  model  masts." 

After  reviewing  the  current  state  of  our  knowledge  of  drag 
forces  on  cylinders  and  the  limitations  of  the  free-flight  cylinder  method 
we  concluded  that  the  original  request  for  both  18  inch  and  24  inch  dia- 
meter cylinders  is  impractical  at  the  higher  diameter  and  is  considered 
not  desirable  in  light  of  the  following  arguments: 

(i)  The  18  inch  cylinder  stretches  the  technical  limitations 
of  our  method,  and  the  24  inch  cylinder  even  more  so.  Our  18  inch  cylin- 
ders are  7.5  feet  long,  weigh  200  pounds  and  need  to  be  centred  6 feet  off 
the  ground  to  avoid  flow  restrictions  between  the  ground  and  the  cylinder. 
The  24  inch  cylinders  would  have  to  be  10  foot  long,  350  pound  objects 
centred  8 feet  off  the  ground.  There  are  problems  fabricating  a cylinder 
of  that  size  to  be  both  strong  and  yet  as  light  in  weight  as  350  pounds. 
The  expected  motion  of  both  these  cylinders  is  beyond  the  recording  limits 
of  our  velocity  transducers,  and  is  getting  beyond  the  reasonable  camera 
range  when  dust  interference  with  the  photography  is  taken  into  account. 
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(ii)  In  present  thinking,  there  are  no  actual  mast  structures 
built  or  currently  under  consideration  which  go  beyond  20  inch  diameter 
The  largest  diameter  target  in  the  past  has  been  the  TACAN  mast,  which 
had  an  O.D.  of  20  inches  in  DIAL  PACK  and  17  inches  in  MIXED  COMPANY. 
For  the  latter  test,  the  measured  strains  were  low,  so  that  larger  dia- 
meter designs  in  the  future  are  not  likely. 

(iii)  The  24  inch  diameter  is  only  33%  larger  than  the  18  inch 
diameter,  and  for  the  region  of  flow  Mach  numbers  and  Reynolds  numbers 
that  apply  to  both  of  them,  we  have  a high  degree  of  confidence  in  ex- 
trapolating the  drag  coefficient  values.  In  fact,  the  one  wind  tunnel 
result  published  in  this  area  shows  very  little  change  in  the  CQ  value 
for  the  flow  region  of  interest. 

(iv)  From  the  review  of  previous  DRES  work,  we  feel  that  the 
problem  of  end  effects,  and  possible  corrections  to  all  previous  DRES 
drag  coefficients  of  about  30",  are  more  important  to  study  this  time 
than  the  24  inch  diameter  cylinder. 

The  experiment  chosen  instead  is  described  in  the  body  of 

the  text. 
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APPENDIX  B 

CYLINDER  MOTION  ESTIMATES 

The  motion  of  each  cylinder  has  been  estimated  in  the  manner 
described  below.  From  the  review  of  Section  2,  values  of  Cp  were  estimated 

as  follows:  for  cylinders  #5  and  #7,  with  end  plates,  at  10  psi  the  value 

of  Cp  = 0.5  was  assumed;  for  cylinders  #3,  #4  and  #6,  at  10  psi  with  no 

end  plates,  the  constant  value  of  Cp  = 0.6  was  chosen;  for  cylinders  #1  and 

#2  with  end  plates  at  20  psi  the  value  of  Cp  was  chosen  to  be  1.3  for  the 

first  20  ms  (M>Mcl),  0.7  for  the  next  20  ms  (M“Mci)  and  0.5  thereafter 

(M<Mci).  These  values  are  all  about  20%  larger  than  actually  expected, 
giving  a conservative  estimate  for  the  motion.  The  diffraction  phase  of 
the  motion  is  then  estimated  using  the  results  of  Long  et  al . (1975)  for 
the  impulse,  IQ,  received  during  diffraction.  Their  result  is  re-expressed 

as*  / p = peak  overpressure  (psi) 

I0(ps1-ms)  = 3.3(p0tT(psi-ms))>-l3  , tra„sn  t1me  (ms) 

This  was  used  as  given  for  the  3.5  in.  diameter  cylinders  but  the  estimate 

was  halved  for  the  9.5  in.  and  18  in.  diameters,  as  it  is  known  that  the 

two  points  which  lie  appreciably  below  the  IQ  curve  of  Long  et  al . (1975) 
are  from  larger  (12  in.  diameter)  cylinders.  The  motion  following  the 
diffraction  phase  was  calculated  using  the  equation: 

F = ma  = CD  A ( 1 44 ) q 

where  F is  the  drag  force  in  Ibp,  m is  the  cylinder  mass  in  slugs,  a the 
acceleration  in  ft/s2,  Cp  the  dimensionless  drag  coefficient,  A the  pro- 
jected cylinder  area  (length  x diameter)  in  ft2,  and  q the  dynamic  pressure 
in  psi.  The  constant  144  converts  the  dynamic  pressure  to  1 bp/ft 2 . The 
motion  was  calculated  using  a four-step  integration  process  over  steps 
covering  20  ms  time  intervals  using  the  q values  presented  in  Table  6, 
taken  from  published  DRES  reports.  The  simple  four-step  integration  is 
able  to  predict  the  actual  motion  of  cylinders  previously  tested  to  within 
about  20%.  The  calculated  cylinder  motions  are  given  in  Table  7. 
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Also  presented  there  is  the  estimate  of  end  plate  drag  for 
each  cylinder.  The  end  plate  drag  is  calculated  for  both  the  inner 
and  the  outer  surfaces  of  each  end  plate.  As  the  air  flow  is  essen- 
tially tangential  to  the  end  plate  surface  the  drag  force  is  fric- 
tional in  nature  and  not  "form"  drag  as  for  the  cylinder. 

Since  the  air  density  and  particle  speed  vary  in  the  same 
fashion  for  both  the  cylinder  and  the  plate,  the  ratio  of  end  plate  drag 
to  cylinder  drag  must  be  equal  to  [Cp(plate)A(plate)]/[Cp(cyl . )A(cyl . )]. 

For  the  frictional  plate  drag,  assuming  the  worst  case  of  turbulent  flow, 
the  Cp  values  from  Hoerner  (1958,  Ch.  2,  Fig.  5)  are  s0.008.  Taking 
flow  compressibility  into  account  yields  about  the  same  value  of  Cp 
(Hoerner,  1958,  Ch.  15,  Fig.  13).  The  laminar  flow  Cp  values  are  appre- 
ciably lower.  The  resulting  estimates  for  end  plate  drag,  presented 
in  Table  7,  show  that  in  the  worst  case  the  additional  drag  force  is 
only  8.9%  of  the  cylinder  drag  force.  In  view  of  the  need  for  end 
plates,  this  is  a reasonable  compromise,  and  the  additional  force  can 
be  estimated  and  subtracted. 
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APPENDIX  C 


BLAST-RESISTANT  CYLINDER  DESIGN 


Most  of  the  cylinders  consist  of  hollow  tubing  with  reinforcing 
discs  and  a welded  centre  rod.  This  design  is  a compromise  between  the 
light  weight  required  for  the  desired  flight  path  and  the  necessary 
strength  to  resist  collapse  during  the  diffraction  phase  loading.  The 
maximum  loading  pressure  on  the  one  side  of  the  cylinder  is  estimated  in 
Table  8 as  the  diffraction  phase  impulse  divided  by  two  transit  times. 

The  maximum  acceleration  and  pressure  during  the  pressure  drag  phase  is 
also  tabulated  for  comparison;  it  is  the  average  value  estimated  over  the 
first  20  ms  of  the  motion. 


As  it  is  difficult  to  estimate  the  one-sided  pressure  for  collapse 
of  a cylinder,  the  design  estimates  were  made  in  another  way.  First,  the 
uniform  external  pressure  required  for  elastic  buckling  of  the  cylinder  was 
calculated  (see  below),  assuming  a cylinder  held  circular  in  cross-section 
at  equally  spaced  intervals  of  length  i.  Then  the  design  safety  factor  was 
calculated  as  the  ratio  (uniform  pressure  for  elastic  buckl ing)/(maximum 
diffraction  pressure).  The  safety  factor  should  be  considered  in  light  of 
the  following  points:  (i)  it  does  not  take  into  account  the  extra  rigidity 
introduced  by  the  welded  centre  rod;  (ii)  it  does  not  take  into  account 
the  extra  tendency  to  buckle  due  to  the  one-sided  loading  and  to  the 
effects  of  dynamic  loading. 

The  formula  for  the  pressure  for  external  collapse  is  taken  from 
Roark  (1954,  Table  XVI,  Formula  Q,  p.  318)  as: 


P' 


0.81 


Et2 


ir 


3 t2 


r 


2 


where  p1 
E 


v 

l 


uniform  external  pressure  for  elastic  collapse  (psi), 

modulus  of  elasticity  = 107  psi,  - Values  quoted  for  A1  taken 

from  Roark,  1954 

Poisson's  ratio  = 0.36  for  aluminum, 
support  disc  spacing  (in.) 
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t = wal 1 thickness  ( in. ) 
and  r = cylinder  body  radius  (in.). 

so  that  wall  thickness  is  the  main 

variable  useful  for  altering  strength  once  the  cylinder  radius  is 
chosen.  In  fact,  the  values  of  t used  in  the  design,  and  presented  in 
Table  8,  were  derived  from  this  calculation. 


Note  that  p’  « 
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APPENDIX  D 


TEST  OF  VELOCITY  CALIBRATION 

The  velocity  transducers  were  calibrated  in  the  laboratory 
prior  to  the  field  trial  using  an  electro-mechanical  method.  Since  the 
circuit  parameters  varied  slightly  in  the  field  from  those  used  in  the 
laboratory,  the  actual  velocity  calibration  may  have  been  in  error  by 
more  than  the  5%  uncertainty  allowed  in  the  data  analysis.  A check  on 
the  accuracy  of  the  calibration  is  provided  by  the  following  arguments. 

Two  positions  of  the  magnet  within  the  transducer  coil  are 
well-defined: 

(a)  the  starting  position,  measured  at  setup  time  in  the 
field  to  an  accuracy  of  ± 1/16  inch; 

(b)  the  position  at  which  the  transducer  signal  changes 
polarity  (crossover  point)  determined  to  an  accuracy 
of  ± 1/16  inch  in  the  laboratory. 

Let  the  travel  distance  (b)  - (a)  be  denoted  by  XQ.  If  the 
time  from  the  start  of  the  transducer  signal  to  the  crossover  point  is 
denoted  by  T,  then: 

T 

J v(t)dt  = XQ  (1) 

o 

The  quantity  on  the  left  can  be  determined  entirely  from  the  velocity- 
time transducer  signal.  T is  determined  directly  to  an  accuracy  of 
± 0.25  msec  by  inspection  of  the  digitized  signal.  The  function  v(t) 
is  the  best-fit  polynomial  curve  fitted  to  the  calibrated  transducer 
signal,  using  the  calibration  value  determined  in  the  laboratory.  If 
the  calibration  value  is  correct,  the  integral  on  the  left-hand  side 
should  equal  the  independently-determined  X on  the  right-hand  side 
within  the  error  of  measurement.  In  Table  11,  the  values  of: 
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T 

f v(t)dt,  and  XQ  for  each  transducer  are  listed,  as 
”0 

well  as  the  percent  difference  expressed  as  a fraction  of  XQ. 
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PART  II  - THEORY  OF  LEAST  SQUARES  (Mathews  1965) 


Consider  a set  of  N experimental  data  points  {y^ } and  a fitting 


function  f with  the  corresponding  set  of  N values  { f ^ . The  least  squares 
fit  criterion  requires  minimization  of  the  function: 


! 

i = l 


(1) 


where  a ■ is  the  standard  uncertainty  in  the  i ' th  data  point  y^ . If  f is 
a function  of  n parameters,  aj — a , 


f = f(ar- an). 


then  the  minimum  in  (1)  is  found  by  taking  the  partial  derivative  of  (1) 
with  respect  to  each  of  the  n parameters. 


N fy,  - ft)  3^ 


i=l  o(: 


3a. 


0.  m = 1 , 2, — n . 


(2) 


( i ) Linear  Least  Squares  Fit 

The  solution  of  the  equations  (2)  can  be  cast  into  a convenient 


form  if  the  f.  are  linear  functions  of  the  parameters  a.;  that  is,  if: 

j 


n 


f . = e C • a 

i im  m 

m=  I 


(3) 


with  known  coefficients  C.m.  Conditions  (2)  can  then  be  rewritten  as: 


C, 


l im  yj  : H j 


c,  c. 


m 1 1 


. m=l , 2,  — -n  (4) 


V 


1 *■  2 

°i 


It  is  now  useful  to  define  a "data  vector"  Y and  a "measurement 
matrix"  M with  components  as  follows: 


N C. 

Y = E Y. 

m i=l  1 


(5) 
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mil 


N 

= l 

1 = 1 


C.  C. 
im  u 


°i‘ 


= M 


fcm 


(6) 


Note  that  Y depends  upon  the  experimental  results  y.  and 
the  errors  ai  whereas  M depends  only  upon  the  errors  c..  M is 
symmetric.  (4)  can  then  be  rewritten  as  the  matrix  equation: 

Y = Ma  (7) 

with  the  solution: 


a = M-i  Y . (8) 

Note  that  the  vector  a is  a unique  solution  for  a given  set 
of  data  (y^,  errors  {o^,  and  given  functional  form  (as  determined 
by  the  coefficients  {C.^}).  The  only  stipulation  was  that  the  fitting 
function  should  be  linear  in  the  coefficients  a^. 

(ii)  Uncertainties  in  Parameters 

To  find  the  uncertainties  in  the  parameters,  it  is  necessary 
to  calculate  the  expected  mean-square  deviations  <(am  - a )2>  or,  more 
generally, 

<(am  - %,)<%  - »,>>  (9) 

where  am  is  a theoretical  quantity  corresponding  to  a grand  average 
value  of  am  over  a very  large  number  of  (theoretically-performed) 
identical  experiments. 

From  (8)  and  (5),  we  have: 

(a--Tm>  ' kE<""W  <*k  - V 
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APPENDIX  E 


PHILOSOPHY  OF  CURVE-FITTING  TO  VELOCITY-TIME  DATA 


PART  I - CHOICE  OF  METHOD 


The  ultimate  goal  of  the  velocity-time  measurement  is  to  obtain 
the  aerodynamic  drag  coefficient  as  a function  of  time,  CQ(t).  A priori, 
one  cannot  claim  a knowledge  of  the  detailed  shape  of  Cp(t).  Drag  pres- 
sure is  defined  as  the  product  of  drag  coefficient  CQ  and  dynamic  pres- 
sure q.  The  dynamic  pressure  q(t)  has  a shape  that  is  well-defined:  it 
is  a monotonical ly  decreasing  function  of  time  whose  slope  is  always 
negative  but  decreases  effectively  to  zero  for  large  times. 

Drag  pressure  is  related  to  the  slope  of  the  velocity-time  curve 
through  the  relations: 

d 

j-  = q and  P0  = s a 

where  v is  cylinder  velocity 

a is  cylinder  acceleration 
m is  cyl inder  mass 
A is  cylinder  frontal  area. 

Since  drag  pressure,  Pp,is  defined  to  be  the  product  of  Cp 

and  q: 

PD(t)  = CD(t)  q(t). 

Pp(t)  can  assume  a variety  of  shapes  depending  upon  the  shape  of  Cp(t). 

The  objective  of  fitting  a curve  to  the  velocity-time  data  is  to  derive 
the  main  features  of  the  drag  pressure-time  curve,  Pp(t),  without  a de- 
tailed fore-knowledge  of  the  true  shape  of  this  function. 

To  fit  the  velocity-time  data,  a polynomial  function  consisting 
of  a power  series  in  time  was  chosen  primarily  for  three  reasons: 
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(1)  If  v(t)  = z a tn,  then  the  acceleration  is  very 

n=0  n 

simply  given  by: 

N 

a(t)  = z n an  tn~ 
n=l  n 

(2)  Because  the  power  series  is  linear  in  the  coefficients 
(a  , — , a^),  a linear  least  squares  criterion  can  be  used  to  deter- 
mine the  best-fit  function.  This  theory  has  the  advantage  that  it  pro- 
vides a unique  answer  for  the  best-fit  function  for  a given  order  of 
polynomial.  A further  advantage  lies  in  the  fact  that  the  theory  pro- 
vides a straightforward  prescription  for  the  statistical  uncertainties 
in  the  coefficients  (aQ, — , a ).  This  in  turn  permits  a simple  expres- 
sion for  the  uncertainty  in  acceleration  to  be  derived  in  terms  of  the 
uncertainties  in  these  coefficients.  The  uncertainty  in  a (t)  for  any 
time  t is  an  important  quantity  because  it  translates  directly  into  an 
uncertainty  in  C^(t).  The  uncertainty  in  C^(t)  is  a measure  of  the  con- 
fidence which  one  can  place  in  the  measured  drag  coefficient  and  is  im- 
portant for  that  reason. 

(3)  Available  evidence  on  the  expected  shape  of  Cg(t),  and 
on  the  known  shape  of  q(t),  dynamic  pressure,  suggests  that,  for  the 
Mach  and  Reynolds  number  ranges  studied  in  this  experiment,  the  varia- 
tion of  Pg(t)  is  sufficiently  smooth  to  be  well -described  by  a low  order 
power  series  in  time. 

A least  squares  criterion  was  applied  to  determine  the  best 
fit  to  the  velocity  data.  All  data  points  were  weighted  equally. 


UNCLASSIFIED 


UNCLASSIFIED 


PART  II  - THEORY  OF  LEAST  SQUARES  (Mathews  1965 


•'  Consider  a set  of  N experimental  data  points  {y..}  and  a fitting 
function  f with  the  corresponding  set  of  N values  {f ^ ) . The  least  squares 
fit  criterion  requires  minimization  of  the  function:  ■ 

; fojJ2  <») 

1»1  2o^ 

where  is  the  standard  uncertainty  in  the  i’th  data  point  y^.  If  f Is 
a function  of  n parameters,  ax — an, 

f * f(at — an), 

then  the  minimum  in  (1)  is  found  by  taking  the  partial  derivative  of  (1) 
with  respect  to  each  of  the  n parameters..  . 

N (y^  - f^|  3^ 

Z »— V = 0.  m = 1,  2,— n . (2) 

1*1  of  3aro  - j 

(i)  Linear  Least  Squares  Fit 

The  solution  of  the  equations  (2)  can  be  cast  into  a convenient 
form  if  the  f.  are  linear  functions  of  the  parameters  a^;  that  is,  if: 


f , = E C.  a _ 

1 m*l  ""  " 


■ (3) 


with  known  coefficients  C.m.  Conditions  (2)  can  then  be  rewritten  as: 
z Cim  y^n  CimCit  a.  . m=l,  2, — n (4) 

* 1 7T  1 1 a 2 

°1  °1 

It  Is  now  useful  to  define  a "data  vector"  Y and  a "measurement 
Mtrla"  M with  components  as  follows: 

I - t -1*  Y,  (5) 

* 1-1  of  1 
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(6) 


Note  that  Y depends  upon  the  experimental  results  y.  and 
the  errors  whereas  M depends  only  upon  the  errors  o^.  M Is 
symmetric.  (4)  can  then  be  rewritten  as  the  matrix  equation: 


Y » Ma  (7) 

with  the  solution: 

a * M-i  Y . (8) 


Note  that  the  vector  a is  a unique  solution  for  a given  set 
of  data  {y . > . errors  {o.},  and  given  functional  form  (as  determined 
by  the  coefficients'  (C.t)):  The  only  stipulation  was  that  the  fitting 
function  should  be  linear  in  the  coefficients  a . 

(11)  Uncertainties  in  Parameters 


To  find  the  uncertainties  in  the  parameters,  it  is  necessary 
to  calculate  the  expected  mean-square  deviations  <(am  - am)2>  or,  more 
generally, 

‘•‘m  * - *,>’  (9) 

where  ?m  is  a theoretical  quantity  corresponding  to  a grand  average 

value  of  a over  a very  large  number  of  (theoretically-performed) 
m 

identical  experiments. 

From  (8)  and  (5),  we  have: 

’ W-W  < V 7k> 


Furthermore, 


<(*j  * *j)(yi  * *1}>  “ °1*  61j 


(10) 

(11) 
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where 


i * J 

i 4 J 


if  the  individual  measurements  y-  are  assumed  to  be  statistically 
independent.  Then: 


‘K  - ' a«;>  =kir„  , (r,,mk- 

kj,p,i  o 


^ili  (M" 1 ) . -i>2.  o.6. . 
mk  2 IP  „ 2 i TJ 


= (M_1),p  MPk 

- 


The  standard  error  in  am  is  given  by: 


Aam  = <(am  - am)2>  = V (M_1)mm  '•  (13) 

In  general,  the  cross  terms  <(an)  - am)(af  “ aii)>  with  are  no^  zer0- 

This  means  that  parameters  am  are  not  statistically  independent  and  their 
errors  are  correlated. 

( i i i ) Uncertainty  in  Quantity  Wh i ch^  is  a Function  of  the  Parameters 
Consider  a ouantity  f which  is  a function  of  the  vector  a 

f = f (a)  O4' 

M af  \ 

4f(1)  ■ k:,  k r* 


( 6f ) 2 = F 


k-1  l33k 


N N 

62a.  + F.  F 
k 4-’ 


af 

af] 

aa, 

aa . 

i/ 

J 

6a.  6a, 


Then  (of)2  = F 


! i#r<4v* " " «»,«.» . 

k=l  3ak|  k i=l  j=l  r i/\  ji  ■ J 
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From  (12),  we  know  that: 


<«2ak>  • <M-l)kk 

(17) 

and 

<6Vj>  - • 

Thus 

N afi  N N ( fi 

(4) 

y <18> 

In  general,  the  parameters  am  are  not  statistically  independent,  and 
the  cross-correlation  terms  (M-1)..  are  not  zero. 

" J 

(iv)  Application  to  Fitting  of  Velocity-Time  Data 

The  functional  form  of  the  velocity-time  trace  was  chosen 
to  be  a power  series  in  time: 

v(t)  = ax  + a2t  + a3t2  + — + antn_1  . 

This  function  is  linear  in  the  parameters  (aj — an)  and  can  be  re- 
written in  the  form  (3),  for  a particular  time,  ti : 


n 


where 


v(t.)  = Z 
1 m=l 
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C . a 
lm  m 


Ch 
Ci  2 

Ci  3 
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The  theory  developed  in  (i),  (ii)  and  (iii)  above  can  then 
be  applied.  The  experimental  data  points  v.  consist  of  a set  of  velo- 
cities {v..}  measured  at  particular  times  { t ^ > . Because  the  uncertain- 
ties o .j  in  these  data  points  are  not  known  a priori,  it  was  felt  that 
the  fairest  way  to  treat  the  data  was  to  assume  that  the  uncertainty 
was  the  same  for  all  data  points.  It  can  then  be  shown  that  the  best- 
fit  parameters  are  independent  of  the  particular  value  of  o.  Once  the 
best-fit  curve  was  determined,  the  root-mean  square  deviation  of  the 
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E-7 


data  points  about  the  best-fit  curve  was  taken  as  the  best  estimate 
of  a. 

From  the  fit  to  the  velocity-time  curve,  one  is  interested 
in  the  first  derivative,  given  by: 

v1  = a2  + 2a3t  + 3a4t2+— + (n-l)antn_2  . 

The  partial  derivatives  required  in  section  (iii)  to  estimate  the  un- 
certainty in  v'  for  various  times  t are  given  by: 
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2t 

(n-1 )tn”2  . 
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DRES  DRAG  FORCE  DATA  - FREE-FLIGHT  CYLINDER  METHOD 


S- 

<D 

“O  -M 
C >>  C 
<0  ID  AJ 
O -»->  -X 
"O  0)  IS) 
(1)Q  C 
•r-  O 
i-  O 


Q +-> 
_j  q: 


u 

(D  ^ 
4->  • 

Q)  C 
E — 
«CJ  — 

o 


s- 

<u 

*o  • 
c o 


>> 

o 


UNCLASSIFIED 


o o in  o o o o 

i—  r—  r—  CM  CM  CM  CM 


c\j  c\j  co  co  cr>  cr> 

CM  cm  in  VO  r—  P- 


co  co  cn 
in  m »—  r- 


00 

O O 

«- *s 

E 

00 

o 

o 

O O 

cn 

C T» 

in  in 

V ' 

i i — 

• 

•>. 

in  ^ cm  co  o in 

i i co  o o in 

oo  o 00 

in  ^ in  in  cm  cm 

i i in  ai  oo 

co  in 

cn  co 



« — *. 

in 

r- 

l/) 

cm  m 

CM 

CM 

E 

/ — .. — . 

r^. 

in 

• 

• 

4-> 

oo 

• • 

• 

r— - 

CM 

C 

X 

in  <3- 

>■ — •* 

<D 

in  o cr>  O CM  CM  N 

1 1 CM  CM  CM 

co  o in  o 

•r- 

MO  liico  lOOl  o 

i i inNinm 

in  o m oo 

O 

a) 

ii  .... 

• • 

• • 

•t— 

> 

r— 

r—  CM 

r-  p— 

o 

<1> 

O 

l/l 

p—**» 

O 

to 

s 

E 

. — v 

in 

o 

E 

r—  CO 

00 

in 

cn 

oo  in 

00 

in 

• 

• 

ro 

o 

i i 

• • 

CO 

• 

r— 

r— 

s- 

i i 



Q 

i i 

CT>  i — r—  CO  CO  1 1 

4-> 

1 1 OlOO^ 

CO  Cn  ^ CO 

+-> 

n in  in  m m i i 

IS) 

i i inNinm 

in  o in  oo 

<U 

to 



r— 



• . 

• • 

CD 

r—  CM 

r—  f— ■ 

<a 

<D 

s 

*»■ — ^ 

> 

</> 

, — ./■ — >. 

IS) 

— . 

cn 

C 

^ > — 

t 

co 

a> 

in 

• 

• 

i i 

• • 

. 

i — 

1— 

i i 

CM 

■« — <* 

' 

i i 

sj  O ID  in  CO  co  N 

i i m^tcoN 

cm  oo  oo  in 

00  ^O  ID  U3  Ol  O 

4-> 

1 1 iruo 

o o in 

to 

l S) 

II  .... 

• • 

. • 

r"“ 

CM  CM 

o in 

cnocoooinin 

r—  r—  i—  r—  O O 

o o 

•r— 

• • 

• . 

. . 

Otc 

CM  00 

in  r-  <3-  o o 

co  ro  ai  oi  n ^ 

>- 

cr>  a> 

Q_ 

Q. 

r— 

i — t—  r—  r-  CM  CM 

CM  CM 

z 

<r 

Q. 

co  co 
o o 


O r—  r—  r—  CO  i — LT) 
ONNNCMNO 


o 

lo  r—  in  in  in  in  o 

O N CM  CM  CM  CM  Q 

r—  r—  r—  r—  r—  r-  LU 

X 


oo 

o: 


>- 

o 


QC 

< 


(_> 

a: 


oo 


co 


=4=  +-  +~ 

CO  O 00  00  00  CO 

Nininininm 
• ••••• 
m co  in  in  in  in 


r-CMCO^tPNOO 


in 

cc 


>- 

o 


on 

Si 

oo 


o o o o 
.... 
in  in  in  in 


in  in  Oi  o 


* 

* * 


assumed;  # due  to  a sideways  force  from  a blast  anomaly,  no  data  were  taken; 
assumed  by  interpolation;  ()  are  camera  data;  + shaded  cylinder 
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TABLE  NO.  2 


Cp(M,L/D)  FOR  FINITE  L/D  CALCULATED  FROM  DIAL  PACK  DATA 


AND  INFINITE  L/D  FROM  THE  OPEN  LITERATURE 


Cp  Values  from  DIAL  PACK-finite  L/D 


Average  Cn 
± S.D.  u 


0.83  ± 0.15 
0.95  ± 0.26 


0.95  ± 0.41 
0.67  ± 0.29 
0.45  ± 0.07 


Note:  (1)  For  DIAL  PACK  cylinders  4-7  and  for  M > 0.30  we  have 

R > 6.5  x 105  so  that  the  flow  conditions  are  always  with 
supercritical  Reynolds  numbers. 

(2)  The  data  are  rough.  The  purpose  of  the  table  is  to 

estimate  the  two  average  quantities  in  the  last  column; 
i.e.,  for  M > M and  M < M with  M = the  first 

Cl  Cl  ci 

critical  Mach  number  = 0.48. 
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TABLE  NO.  3 

SUMMARY  OF  L/D  RATIO  EFFECT  VERSUS  MACH  AND  REYNOLDS  NUMBERS 


From  previous  steady  and  unsteady  flow  experiments, 

Cp(L/D  = finite  ) 

n,  n'»  n"  = where 

Cp(L/D  = infinite) 

r\,  n'  and  n"  are  all  the  above  same  ratio  for  different  flow  conditions 
defined  as: 

_ CD(finite  L/D,  steady  flow)  _ _ .unstMdy.  _ .u„stea<ly. 

” * n 9 n • 

Cp( infinite  L/D,  steady  flow)  "steady"  "unsteady 

, Cn(finite  L/D,  unsteady) 

We  also  define  $ = — = . 

n Cp(finite  L/D,  steady) 

- Finite  L/D  means  L/D  in  the  range  5-20  over  which  C^  does  not  change  very 
much. 

- Different  regions  of  M and  R values  are  denoted  by  Roman  numerals  I -VI. 
Only  regions  I - 1 1 1 are  of  interest  to  us. 

- References  are  given  in  brackets  (X)  and  listed  below. 

References  for  Table  No.  3 

A - Martin  et  al.  (1967) 

B - Martin  et  al . (1965) 

C - MellsenTWl)  - (STP  382)  DIAL  PACK  results  analyzed  by  Hill 
D - Mellsen  and  Naylor  (1969)  shock  tube  work 
E - Delaney  and  Sorenson  (1953),  analyzed  by  Hill 
F - Welsh  (1953)  - some  analyzed  by  Hill 
G - Gowen  and  Perkins  (1953) 

H - Naylor  and  Mellsen  (1973)  Summary  of  PRAIRIE  FLAT  and  DIAL  PACK 
I - Long  and  Laidlaw  (1973)  from  the  Lattice  Mast  II  report  - a direct 
comparison  of  Naylor's  Cn(“)  values  with  Mellsen's  Cn(finite); 
analyzed  by  Hill  u u 

J - Hoerner  (1958) 
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(see  text) 


table  (120+ms) 
overlap 
di  ffraction 


end  effects 
overlap  (70  ms) 
diffraction 


table  (100  ms) 
end  effects 
unsteady 


support  VHF  mast(70  ms) 
table 

end  effects 
unsteady 


support  VHF  mast(70  ms) 
table 

end  effect 
unsteady 
diffraction 

support  lattice  mast 
(50  ms) 

end  effects 
unsteady 


support  lattice  mast 
and  whip  antenna 
(50+  ms) 
end  effects 
unsteady 
overlap 
diffraction 
test 
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TABLE  NO.  5 


CONTRIBUTIONS  OF  EACH  CYLINDER 

TO  THE  GOALS 

GOAL 

CYLINDER  NOS. 

1 . 

Support  present  projects 

mast 

6,  7 

VHF/UHF 

4,  5 

whip 

2,  6,  7 

2. 

Table  of  data  to  support 

M<M 

r 

3,  4,  5 

future  masts  for  diffraction 

v» 

M>M 

1 2 

phase  and  Cp(t) 

1 1^1  1 

C 

3. 

End  Effects:  finite  L/D 

M>MC 

2 

vs  infinite 

M<MC 

3,  4,  5,  6,  7 

4. 

Unsteady  vs  steady  flow 

M>M 

c 

3 vs  ref.M 

M<M„ 

c 

4,  5,  6,  7 vs  ref. I 

5. 

With  unsteady  flow,  overlap 

1 wi th  3 

M>M  with  M<M„ 
c c 

2 with  6 and  7 

6. 

Scientific  study  of 

M>M_ 

1,  2 

diffraction  for 

M<M 

5 7 

infinite  length  (R>Rc) 

1 1 N 1 1 

C 

Other 

7. 

Test  needed  diameter  of  end  plate 

2,  6,  7 

8. 

Brand  new  regions  of  measurement 

1,  3 

9. 

Some  overlap  with  previous  work 

2.  4 5 6.  7 

(DRES  and  open  literature) 

u 9 “ 9 J u » ’ 
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TABLE  NO.  6 

VALUES  OF  DYNAMIC  PRESSURE  USED  IN  CYLINDER  MOTION  ESTIMATES 


Over- 

pressure 

(psi) 


Dynamic  Pressure,  q (psi ) 

Time  Period 

(ms) 

0 - 20 

20  - 40 

40  - 60 

60  - 80 

80  - 100 

6.51 

3.24 

1.57 

0.73 

— 

2.00 

1.27 

0.77 

0.46 

0.25 

1 .2 

0.80 

0.54 

0.35 

0.22 

Reference 


Values  taken  at 
mid  point  of  time 
interval 


Taken  equal  to  20.5 
psi  results  from 
Dial  Pack  cylinder  #7 

Taken  as  average  of 
Mixed  Company  #4 
(9.1  psi ) and  Dial 
Pack  #2  (11  psi ) 

Taken  from  Mixed 
Company  #10  at 
7.0  psi 
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CALCULATED  CYLINDER  MOTION 
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CALCULATIONS  FOR  BLAST  RESISTANT  CYLINDER  DESIGN 
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t Designed  to  withstand  10  psi  overpressure;  tested  at  7 psi  peak  overpressure. 
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SUMMARY  OF  CAMERAS  USED  IN  TRIAL 


Transducer 

Designation 

c 

rT 

v(t)dt  (inches) 

) 

XQ  (inches) 
(±0.06  inches) 

Per  cent 

Difference 

2 East  Broken  magnet  - 

no  signal 

9.65 

— 

2 West 

12.7 

9.65 

31.6  ** 

3 East 

10.34 

9.65 

7.2 

3 West 

9.76 

9.65 

1.1 

4 East 

10.23 

9.86 

3.8 

4 West 

9.75 

9.68 

0.1 

5 East 

10.17 

9.80 

3.7 

5 West 

9.90 

9.84 

0.6 

6 East 

9.94 

9.60 

3.5 

6 West 

9.73 

9.60 

1 , 4 

7 East 

9.77 

9.60 

1.8 

7 West 

9.68 

9.65 

0.3 

* See  Appendix  D 

**  Poor  v(t)  data  (poor  signal/noise,  large  oscillations, 
truncate d record) 
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TABLE  NO.  12 

COMPARISON  OF  FITTING  METHODS 

METHOD  A - Mean  Square  Minimization  by  (Linear)  Theory  of  Least  Squares 
Fitting  Function:  v(t)  = at  + a2t  + a3t2,  v in  ft/sec,  t in  sec. 

METHOD  B - Mean  Square  Minimization  by  Parameter  Search 

Fitting  Function:  v(t)  = aj  + a2t  + a3t  + a4e"a5t  sin(2ira6t  + a7), 

v in  ft/sec,  t in  sec. 


DATA 

BEST-FIT  COEFFICIENTS 

A.  Linear  Least  Squares 

B.  Parameter  Search 

Cylinder  3 (East  End) 

ax  = 9.0  ± 1.0 

6.9 

a2  = 224  ± 45 

308 

a3  = -962  ±412 

-1583 

Cylinder  3 (West  End) 

a!  = 8. 00±  0.8 

6.5 

a2  = 262  ± 35 

325 

a3  =-1350  ±318 

-1817 

Cylinder  4 (East  End) 

aj  = 4.5  ± 0.3 

4.0 

a2  = 93  ±9 

108 

a3  = -342  ± 43 

-367 

Cylinder  4 (West  End) 

a!  = 4.3  ± 0.2 

4.0 

a2  =112  ±6 

119 

a3  = -420  ± 34 

-408 
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TABLE  NO.  13 


BEST-FIT  POLYNOMIALS 


FOR  VELOCITY  TRANSDUCER  DATA 


v(t)  = ai  + a2t  + a3t2  + 


— +an 


v in  ft/sec,  t in  sec. 


Transducer 

Polynomial 

Best-Fit  Coefficients 

Designation 

Order 

Cylinder  End 

Number 

Value 

Uncertainty 

(1  Std  Deviation) 

2 West 

2 

1 

8.3904x10° 

1.0016x10° 

2 

7.0125xl02 

5.4556X101 

3 

| 

-5. 3077xl03 

6.0782x102 

3 East 

4 

i 

8.6017x10° 

2.0015x10° 

1 

2 

3.0837xl02 

2.5543xl02 

3 

-5.0942xl03 

9.7309xl03 

4 

6. 7784xl04 

1.3828xl05 

5 

-3. 5039xl05 

6.5335x10s 

3 West 

4 

1 

8.2889x10° 

1.5603x10° 

2.7845x102 

1.9723xl02 

il 

-4. 1 275xl03 

7. 4484x1 03 

1 

6. 3368x1 04 

1.0497x10s 

| 

Hi 

-3.8457x10s 

4.9172x10s 

4 East 

4 

i 

i 

4.1260x10° 

0.5768x10° 

H 

1.1047X102 

0.4570xl02 

-3.2360xl02 

10.667  xlO2 

II 

-3. 3886xl03 

9. 2508xl03 

Hi 

1. 7292x1 04 

2.6663xl04 

4 West 

4 

i 

3.5314x10° 

0.3704x10° 

1.8431X102 

0.3103x102 

-2.1 854x1 03 

0. 7650xl03 

1. 5365x1 04 

0.6972xl04 

-4. 3843xl04 

2.1 055x1 04 
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TABLE  NO.  13  (Cont'd) 
BEST-FIT  POLYNOMIALS 
FOR  VELOCITY  TRANSDUCER  DATA 


v(t)  = a*  + a2t  + a3t2  + — -+a  tn_1,  v in  ft/sec,  t in  sec. 

n 


Transducer 

Designation 

Cylinder  End 

Polynomial 

Order 

Best-Fit  Coefficients 

Number 

Value 

Uncertainty 

5 East 

2 

1 

2.7445x10° 

3.8581 xlO"  1 

2 

1.0981xl02 

1.5367X101 

3 

-4. 8908x1 02 

1.2920xl02 

5 West 

2 

1 

2.3405x10° 

2. 941 5x1  O' 1 

2 

1.4965/102 

1. 3521x10° 

3 

-6. 3441x1 02 

1.31 36x1 02 

6 East 

4 

1 

1.9525x10° 

0.2717x10° 

2 

2. 5509xl02 

0. 2378xl02 

3 

-2.3499xl03 

0.6169x10° 

4 

2.0936xl04 

0. 5970xl04 

5 

-6.1 043x1 04 

1 . 921 2x1 04 

6 West 

4 

1 

2,0853x10° 

0.3312x10° 

2 

2. 8674xl02 

0. 3252xl02 

3 

-2.5961x10° 

0.9472xl03 

4 

1. 5990x1 04 

1 . 0289xl04 

5 

-4. 3328x1 O4 

3.  71 95x1 04 

7 East 

4 

1 

2.2995x10° 

0.3017x10° 

2 

2. 5942xl02 

0. 2672xl02 

3 

-2.6452x10° 

0.6995x1 0° 

4 

1.4874xl04 

0.6795xl04 

5 

-3.3221 xlO4 

2.1 900x1 O4 

7 West 

4 

1 

2.5874x10° 

0.3463x10° 

2 

2. 5647xl02 

0. 3235xl02 

3 

-2.3679x10° 

0.8945x10° 

4 

1 . 31 72xl04 

0. 9 1 93x1 04 

5 

-3.2238xl04 

3. 1 373xl04 

UNCLASSIFIED 


UNCLASSIFIED 


TABLE  NO.  14 
BEST-FIT  POLYNOMIALS 
FOR  HIGH-SPEED  CAMERA  DATA 

x(t)  = aj  + a2t  + a3t2+  +antn_1,  x in  feet,  t in  seconds 


Cylinder 

Desig- 

nation 

Polynomial 

Order 

Best-Fit  Coefficients 

Number 

Value 

Uncertainty 

3 

5 

1 

-2.1168xl0"2 

0.2699xl0~2 

2 

8.8515x100 

0.5018x10° 

3 

1.7282x102 

0.2374x102 

4 

-2.3930x103 

0.6772x103 

5 

2.3814x10^ 

0.6960xl04 

6 

-8.7640xl04 

2. 5890x1 04 

4 

5 ' 

nj 

1.1 682x1 0"** 

7 . 1 695x1 0~4 

2.6303x10° 

0.1075x10° 

1.0103x102 

0.0496x10?. 

-5.5793x102 

0.9418x102 

1.2340x103 

0.7795x103 

Kfl 

-8.0000x10° 

2.3327x103 

5 

5 

mm 

-1 . 9645xl0“3 

1 . 7078xl0"3 

1.4761x100 

0.2673x100 

He1 

8.5265x101 

1.2886x101 

raj 

-4.3162x102 

2.5558x102 

1.5490x103 

2.2105x103 

H 

-3.5520x103 

6.9140x103 

6 

5 

Hi 

-2. 559xl0”3 

1.4469xl0'3 

5.4772x10'! 

2.7658x10'! 

Bp  a 

2.0469x102 

0.1525x102 

|| 

-2.3431x103 

0.3582x103 

1 . 951 6x1 0*4 

0. 3671 xlO4 

B 

-6. 7472xl04 

1 . 3609xl04 

7 

3 

i 

-1 . 1862xl0"2 

0.1 723xl0"2 

3.3434x10° 

0.1541x100 

1.0242x102 

0.0371x102 

B 

-2.7622x102 

0.2531x102 
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NO  2 3.5  IN  DIAM  20.1  PSI 

PFAK  OVFRPRFSSURF  20.1  PSI 


FRIFOLANDFR  OFCAY  CONSTANT  * 2.23 


10.5  IN  DIAM  ENO  PLATES 
POSITIVE  DURATION  190.0  MSEC 


timf  after 

DRAG 

arsolute 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

ELUW 

REYNOLDS 

shock  front 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFE I C I ENT 

FRRUK** 

MACH 

No. 

arrival 

PD 

0 

QI 

CD 

NO. 

R 

(MSEC) 

(PSI  ) 

(PSI  ) 

(PERCENT) 

(PSI  ) 

(PSI  ) 

(OI/O) 

(PD/QI ) 

M 

( X 10-5  ) 

2.0 

6.364 

1.783 

28.0 

8.81 

9.72 

1.103 

0.654 

0.183 

0.628 

13.87 

4.0 

6.166 

1.707 

27.6 

8.23 

9,08 

1.103 

0.678 

0.187 

0.614 

13.47 

6.0 

3.967 

1.631 

27*3 

7.66 

a.4e 

1.103 

0.703 

0.192 

0.599 

13. C9 

3.0 

5.768 

1.555 

26.9 

7.17 

7.91 

1.103 

0.72b 

0.196 

C.505 

12.71 

10.0 

3.570 

1.480 

26.5 

6.69 

7,39 

1.103 

0.753 

0.20C 

0.570 

12.34 

12.0 

3.371 

1.404 

26.1 

6.24 

6.89 

1.103 

0.778 

C.  203 

0.556 

11.96 

14.0 

5.172 

1.330 

25.7 

5.82 

6.43 

1.103 

0.8C3 

0.206 

0.543 

11.63 

16. 0 

4.973 

1.255 

25.2 

5.43 

5.99 

1.103 

0.828 

0.209 

0.529 

11.29 

1 A.  0 

4.775 

1.181 

24.7 

5.06 

5.59 

1.103 

0.853 

0.211 

C • 5 1 6 

10.95 

20.0 

4.576 

1.108 

74.2 

4.72 

5.21 

1.103 

0.877 

0.212 

0.503 

10.63 

22.0 

4.377 

1.036 

23.6 

4.40 

4.85 

1.103 

0.901 

0.213 

C.490 

10.31 

24.0 

4.179 

0.964 

73.0 

4.10 

4.52 

1.103 

0.923 

0.213 

0.477 

9.99 

26.0 

3.9P0 

0.894 

22.4 

3.81 

4.21 

1.103 

0.944 

0.212 

0.465 

9.69 

28.0 

3.781 

0.826 

21.8 

3.55 

3.92 

1.103 

0.963 

C.210 

0.452 

9.39 

30.0 

3.582 

0.760 

21.2 

3.30 

3.65 

1.103 

0.980 

0.208 

0.440 

9.10 

32.0 

3.384 

0.696 

20.3 

3.07 

3.39 

1.103 

0.995 

C.2C4 

0.428 

6.82 

34.0 

3.185 

0.636 

19.9 

2.86 

3.16 

1.103 

1.0C7 

0.201 

0.417 

8.54 

36.0 

2.986 

0.582 

19.4 

2.66 

2.93 

1.103 

1.016 

0.198 

0.405 

B. 27 

38.0 

2.787 

0.534 

19.1 

2.47 

2.73 

1.103 

1.020 

0.195 

0.394 

8.00 

40.0 

7.589 

0.495 

19.1 

2.29 

2.53 

1.103 

1.020 

0.195 

C.383 

7.75 

42.0 

7.390 

0.467 

19.5 

2.13 

2.35 

1.103 

1.013 

0.19e 

0.372 

7.50 

44.0 

7.191 

0.452 

20  • 6 

1.98 

2.18 

1.103 

1.000 

0.206 

0.361 

7.25 

46.0 

1.993 

0.431 

22.6 

1.84 

2.03 

1.103 

0.980 

0.221 

0.351 

7.01 

46.0 

1.794 

0.467 

75.7 

1.70 

1.88 

1.103 

0.951 

0.245 

C.341 

6.76 

30.0 

1.595 

0.483 

30.2 

1.58 

1.74 

1.103 

0.912 

0.276 

0.331 

6.55 

32.0 

1.396 

0.513 

36.7 

1.46 

1.62 

1.103 

0.861 

0.316 

0.321 

6.33 

34.0 

1.198 

0.550 

45.9 

1.36 

1.50 

1.103 

C.797 

0.366 

0.311 

6.12 

36.0 

0.999 

0.59C 

59.0 

1.26 

1.39 

1.103 

0.718 

0.424 

0.301 

5.91 

36.0 

0.800 

0.635 

79.3 

1.16 

1.28 

1.103 

0.621 

0.492 

C.  292 

5.70 

60.0 

0.602 

0.682 

113.2 

1 .08 

1.19 

1.103 

0.504 

0.572 

0.283 

5.50 

62.0 

0.403 

0.731 

181.3 

0.99 

1.10 

1.103 

0.365 

0.662 

0.274 

5.31 

64.0 

0.204 

0.782 

383.3 

0.92 

1.01 

1.103 

0.200 

0.767 

0.265 

5.12 

66.0 

0.005 

0.833 

16660.0 

0.85 

0.94 

1.103 

0.0C5 

0.884 

C.  257 

4.94 

68.0 

-0.192 

0.886 

-461.4 

0.78 

0.87 

1.103 

-0.220 

i.oie 

C.248 

4.76 

70.0 

-0.391 

0. 940 

-240.4 

0.72 

0.80 

1.103 

-0.486 

1.170 

0.240 

4.59 

72.0 

-0.590 

0.994 

-168.4 

0.67 

0.74 

1.103 

-0.796 

1.341 

0.23  2 

4.*2 

74.0 

-0.789 

1.048 

-132.8 

0.61 

0.68 

1.103 

-1.154 

1.533 

0.224 

4.25 

76.0 

-0.987 

1.103 

-111.7 

0.57 

0.62 

1.103 

-1.567 

1.751 

0.216 

4.09 

78.0 

-1.186 

1.139 

-97.7 

0.52 

0.57 

1.103 

-2.045 

1.998 

0.209 

3.94 

80.0 

-1.385 

1.213 

-87,7 

0.48 

0.53 

1.103 

-2.594 

2.276 

C.2C1 

3.79 

82.0 

-1.583 

1.271 

-80.2 

0.44 

0.49 

1.103 

-3.224 

2.589 

0.194 

3.64 

• • 3 STANDARD  DEVIATIONS 


TABLE  NO.  16 


UNCLASSIFIED 


UNCLASSIFIED 


STP  453 


NO  3 1(1.0  IN  DI AM  9.7  PSI  1A.0  IN  DIAM  END  PLATES 

PEAK  OVFRPRESSURF  9.7  PSI  POSITIVE  DURATION  233.0  MSEC 

FR I FOLANDFR  DFCAY  CONSTANT  • 0.83 


TI*F  after 

DRAG 

ABSOLUTF 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

FLUW 

REYNOLDS 

SHOCK  PROMT 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PD 

Q 

01 

CD 

NO. 

R 

( MSEC ! 

(PSI  ) 

(PSI  J 

(PERCENT! 

(PSD 

(PSI  ' 

(01/0) 

(PD/OI ) 

M 

( X 10-5 ) 

2.0 

0.847 

1.281 

151.1 

2.35 

2.45 

1.039 

0.345 

0.522 

0.391 

42.66 

4.0 

0.79P 

1.091 

136.6 

2.28 

2.37 

1.039 

0.336 

0.459 

0.386 

42.06 

6*0 

C.734 

0.920 

122.0 

2.21 

2.30 

1.039 

0.327 

0.400 

0.382 

41.46 

8.0 

0.714 

0.767 

107,4 

2.14 

2.22 

1.039 

0.320 

0.344 

C.  377 

40.87 

10.0 

0.678 

0.630 

93.0 

2.07 

2.15 

1.039 

C.  314 

0.292 

0.37  2 

40.29 

12.0 

0.645 

0.512 

79.4 

2.01 

2.08 

1.039 

0.308 

0.245 

0.367 

39.70 

14. C 

C.616 

0.412 

66.9 

1.94 

2.02 

1.039 

C.3C4 

0.204 

C.363 

39.13 

16.0 

0.590 

0.332 

56.2 

1.88 

1.95 

1.039 

0.301 

0.169 

0.358 

38.56 

18.0 

0.568 

0.272 

47.9 

1.82 

1.89 

1.039 

0.299 

0.143 

0.353 

37.99 

20.0 

0.547 

0.235 

43.0 

1.76 

1.83 

1.039 

0.298 

0.126 

C.  349 

37.43 

22.0 

0.530 

0.220 

41.5 

1.70 

1.77 

1.039 

C.299 

0.124 

0.344 

36.87 

24.0 

0.515 

0.220 

42.7 

1.65 

1.71 

1.039 

0 . 3 0 C 

0.128 

C.  339 

36.32 

26.0 

0.503 

0.229 

45.5 

1.59 

1.65 

1.039 

0.303 

0.13e 

0.335 

35.77 

26.0 

0.492 

0.239 

48.5 

1.54 

1.60 

1.C39 

0.3C7 

0.149 

0.330 

35.22 

*0.0 

0.484 

0.248 

51.3 

1.49 

1.55 

1.C39 

0.312 

0.160 

0.326 

34.69 

32.0 

0.476 

0.253 

53.2 

1.44 

1.49 

1.C39 

0.317 

0.169 

0.321 

34.15 

34.0 

0.469 

0.255 

54.2 

1.39 

1.44 

1.039 

C.324 

0.176 

0.317 

33.62 

36.0 

0.465 

0.252 

54.1 

1.34 

1.39 

1.C39 

0.332 

0.180 

0.312 

33.10 

38.0 

0.461 

0.245 

53.1 

1.30 

1.35 

1.039 

C.341 

0.181 

0.308 

32.57 

40.0 

0.457 

0.234 

51.2 

1.25 

1.3C 

1.039 

0.350 

0.179 

0.304 

32.C6 

42.0 

0.454 

0.221 

48.7 

1.21 

1.26 

1.039 

0.360 

0.175 

0.299 

31.55 

44.0 

0.451 

0.207 

45.9 

1.17 

1.21 

1.039 

0.371 

0.170 

0.295 

31.04 

46.0 

0.448 

0.192 

42.8 

1.12 

1.17 

1.039 

0.381 

0.163 

0.291 

30.54 

48.0 

0.445 

0.178 

39.9 

1.09 

1.13 

1.039 

0.393 

0.157 

0.286 

30. C4 

30.0 

0.441 

0.167 

37.8 

1.05 

1.09 

1.C39 

0.404 

0.153 

0.282 

29.54 

32.0 

0.437 

0.160 

36.6 

1.01 

1.05 

1.C39 

0.415 

0.152 

0.278 

29.05 

34.0 

0.433 

0.159 

36.7 

0.97 

1.01 

1.039 

0.426 

0.156 

C.274 

28.57 

36.0 

0.426 

0.163 

38.4 

0.94 

0.97 

1.039 

0.435 

C.  167 

0.270 

2e.C9 

38.0 

0.419 

0.173 

41.3 

0.90 

0.94 

1.039 

0.443 

0.183 

0.265 

27.61 

60.0 

0.410 

0.185 

45.1 

0.87 

0.90 

1.039 

0.451 

0.203 

0.261 

27.14 

62.0 

0.399 

0.198 

49.7 

0.84 

0.87 

1.C39 

0.455 

0.227 

0.257 

26.67 

64.0 

0.385 

0.212 

55.0 

0.81 

0.84 

1.039 

0.457 

0.251 

0.253 

26.20 

66.0 

0.370 

0.224 

60.6 

0.78 

0.81 

1.039 

0.455 

0.276 

0.249 

25.74 

68.0 

0.353 

0.233 

66.2 

0.75 

0.78 

1.039 

0.451 

0.299 

0.245 

25.29 

70.0 

0.332 

0.239 

72.0 

0.72 

0.75 

1.039 

0.442 

0.318 

0.241 

24.84 

72.0 

0.309 

0.241 

78.2 

0.69 

0.72 

1.039 

0.427 

0.334 

0.237 

24.39 

74.0 

0.282 

0.240 

85.0 

0.66 

0.69 

1.039 

0.406 

0.345 

0.233 

23.95 

76.0 

0.252 

0.234 

93.0 

0.64 

0.66 

1.039 

0.377 

0.351 

0.229 

23.51 

78.0 

0.219 

0.226 

103.2 

0.61 

0.64 

1.039 

0.341 

0.352 

0.225 

23.07 

80.0 

0.181 

0.218 

120.5 

0.59 

0.61 

1.039 

0.293 

0.353 

0.221 

22.64 

82.0 

0.138 

0.213 

153.6 

0.56 

0.59 

1.039 

0.234 

0.360 

0.218 

22.21 

84.0 

0.091 

0.218 

239.6 

0.54 

0.56 

1.039 

0.160 

0.383 

0.214 

21.79 

86.0 

0.038 

0.239 

614.4 

0.52 

0.54 

1.039 

0.071 

0.438 

0.210 

21.37 

88.0 

-0.018 

0.279 

-1473.3 

0.50 

0.52 

1.039 

-0.036 

0.533 

0.206 

20.96 

90.0 

-0.080 

0.340 

-424. 1 

0.48 

0.50 

1.039 

-0.159 

0.677 

0.202 

20.55 

92.0 

-0.146 

0.420 

-286.4 

0.46 

0.48 

1.039 

-0.305 

0.873 

0.199 

20.14 

94.0 

-0.217 

0.519 

-238.6 

0.44 

0.46 

1.039 

-0.472 

1.126 

0.195 

19.73 

96.0 

-0.293 

0.635 

-216.5 

0.42 

0.44 

1.039 

-0.665 

1.940 

C.191 

19.34 

98.0 

-0.376 

0.769 

-204.3 

0.40 

0.42 

1.039 

-0.890 

1.820 

0.188 

18.94 

00.0 

-0.464 

0.920 

-198.1 

0.38 

0.40 

1.039 

-1.148 

2.276 

0.184 

18.55 

••  3 STANDARD  DEVIATIONS 

TABLE  NO.  17 

DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  3 - VELOCITY  TRANSDUCER  DATA 

UNCLASSIFIED 


UNCLASSIFIED 


NO  3 18.0  IN  DIANf  9.7  PSI 

PEAK  OVERPRESSURE  9.7  PSI 
FRIEDLANDER  DECAY  CONSTANT  » 0.83 


18.0  IN  D I AM  END  PLATES 
POSITIVE  DURATION  '33.0  MSEC 


STP  453 


TIME  AFTER 

DRAG 

ABSOLUTE 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

FLOW 

REYNOLDS 

SHOCK  FRONT 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PO 

0 

QI 

CD 

NO. 

R 

(MSEC) 

(PSI  ) 

(PSI  > 

(PERCENT) 

(PSI  ) 

(PSI  ) 

(QI/Q) 

(PD/QI ) 

M 

( X 10-5 ) 

2.0 

0.602 

0.163 

27.0 

2.35 

2.45 

1.039 

0.245 

0.066 

0.391 

42.66 

4.0 

0.592 

0*157 

26.5 

2.28 

2.37 

1.039 

0.249 

0.066 

0.366 

42.06 

6.0 

0.582 

0.151 

25.9 

2.2 1 

2.30 

1.039 

0.252 

0.065 

0.382 

41.46 

8.0 

0.571 

0.146 

25.5 

2.14 

2.22 

1.039 

0.256 

0.065 

0.377 

40.87 

10.0 

0.561 

0.140 

24.9 

2.07 

2.15 

1.039 

0.259 

0.064 

0.372 

40.29 

12.0 

0.550 

0.134 

24.3 

2.01 

2.06 

1.039 

0.263 

0.064 

0.367 

39.70 

14.0 

0.540 

0.129 

23.6 

1.94 

2.02 

1.039 

0.266 

0.063 

0.363 

39.13 

16.0 

0.530 

0.123 

23.2 

1.86 

1.95 

1.039 

0.270 

0.062 

0.358 

38.56 

16.0 

0.519 

0.118 

22.7 

1.82 

1.89 

1.039 

0.273 

0.062 

0.353 

37.99 

20.0 

0.509 

0.112 

22.0 

1.76 

1.83 

1.039 

0.277 

0.061 

0.349 

37.43 

22.0 

0.498 

0.107 

21.4 

1.70 

1.77 

1.039 

0.280 

0.C60 

0.344 

36.87 

24.0 

0.488 

0.101 

20.6 

1.65 

1.71 

1.039 

0.284 

0.058 

0.339 

36.32 

26.0 

0.478 

0.096 

20.0 

1.59 

1.65 

1.039 

0.286 

0.057 

0.335 

35.77 

28.0 

0.467 

0.090 

19.2 

1.54 

1.60 

1.039 

C.291 

0.056 

0.330 

35.22 

80.0 

0.457 

0.085 

18.5 

1.49 

1.55 

1.039 

0*294 

0.054 

0.326 

34.69 

32.0 

0.446 

0.080 

17.9 

1.44 

1.49 

1.039 

0.297 

0.053 

0.321 

34.15 

84.0 

0.436 

0.075 

17.2 

1.39 

1.44 

1.039 

0.301 

0.051 

0.317 

33.62 

86.0 

0.426 

o.o*»i 

16.6 

1.34 

1.39 

1.039 

0.304 

0.050 

0.312 

33.10 

88.0 

0.415 

0.066 

15.9 

1.30 

1.35 

1.039 

0.307 

0*048 

C.  308 

32.57 

40.0 

0.405 

0.062 

15.3 

1.25 

1.30 

1.039 

0.310 

0.047 

0.304 

32.06 

42.0 

0.394 

0.058 

14.7 

1.21 

1.26 

1.039 

0.312 

0.046 

0.299 

31.55 

44.0 

0.384 

0.055 

14.3 

1.17 

1.21 

1.039 

0.315 

0.045 

0.295 

31.04 

46.0 

0.374 

0.053 

14.1 

1.12 

1.17 

1.039 

0.318 

0.045 

0.291 

30.54 

48.0 

0.363 

0.051 

14.0 

1.09 

1.13 

1.039 

0.320 

0.045 

0.266 

30.04 

80.0 

0.353 

0.050 

14.1 

1.05 

1.09 

1.039 

0*323 

0.045 

0.282 

29.54 

82.0 

0.342 

0.050 

14.6 

1.01 

1.05 

1.039 

0.324 

0.047 

C.  278 

29.05 

84.0 

0.332 

0.051 

15.3 

0.97 

1.01 

1.039 

0*326 

0.050 

0.274 

28.57 

86.0 

0.322 

0.052 

16.1 

0.94 

0.97 

1.039 

0.328 

0.053 

0.270 

28.09 

88.0 

0.311 

0.054 

17.3 

0.90 

0.94 

1.039 

0.329 

0.057 

0*265 

27.61 

60.0 

0.301 

0.057 

18.9 

0.87 

0.90 

1.039 

0.330 

0.062 

0.261 

27.14 

62.0 

0.290 

0.060 

20.6 

0.84 

0.87 

1.039 

0.331 

0.068 

0*257 

26.67 

64*0 

0.280 

0.063 

22.5 

0.81 

0.84 

1.039 

0.331 

0.074 

0.253 

26.20 

66.0 

0.270 

0.067 

24.8 

0.78 

0.81 

1.039 

0.332 

0.C82 

0.249 

25.74 

68.0 

0.259 

0.071 

27.4 

0.75 

0.78 

1.039 

0.331 

0.090 

0.245 

25.29 

70.0 

0.249 

0.075 

30.1 

0.72 

0.75 

1.039 

0.331 

0.099 

0.241 

24.84 

72.0 

0.238 

0.079 

33.1 

0.69 

0.72 

1.039 

0.329 

0.109 

0.237 

24.39 

74.0 

0.228 

0.083 

36.4 

0.66 

0.69 

1.039 

0.327 

0.119 

0.233 

23.95 

76.0 

0.218 

0.088 

40.3 

0.64 

0.66 

1.039 

0.326 

0.131 

0.229 

23.51 

78.0 

0.207 

0.092 

44.4 

0*61 

0.64 

1.039 

0.322 

0.143 

0.225 

23.07 

80.0 

0.197 

0.096 

48.7 

0.59 

0.61 

1.039 

0.319 

0.155 

0.221 

22.64 

82.0 

0.186 

0.101 

54.3 

0.56 

0.59 

1.039 

0.313 

0.170 

0.218 

22.21 

84.0 

0.176 

0.105 

59.6 

0.54 

0.56 

1.039 

0.309 

0.164 

0.214 

21.79 

86.0 

0.166 

0.110 

66.2 

0*52 

0.54 

1.039 

0.304 

0.201 

0.210 

21.37 

88.0 

0.155 

0.115 

74.1 

0*50 

0.52 

1.039 

0.296 

0.219 

0.206 

20.96 

90.0 

0.145 

0.119 

82.0 

0.48 

0.50 

1.039 

0.288 

0.237 

0.202 

20.55 

92.0 

0.134 

0.124 

92.5 

0.46 

0.48 

1.039 

0.278 

0.257 

0.199 

20.14 

94.0 

0.124 

0.128 

103.2 

0*44 

0.46 

1.039 

0.269 

0.277 

o.m 

19.73 

96.0 

0.114 

0.133 

116.6 

0*42 

0.44 

1.039 

0.258 

0.301 

0.191 

19.34 

98.0 

0.057 

0.081 

142.1 

0.40 

0.42 

1.039 

0.134 

0.191 

0.186 

18.94 

100.0 

0.049 

0.084 

171.4 

0.38 

0.40 

1.039 

0.121 

0.207 

0.184 

18.55 

102.0 

0.041 

0.087 

212.1 

0.37 

0.38 

1.039 

0.106 

0.225 

0.181 

16.16 

104.0 

0.033 

0.090 

272.7 

0.35 

0.36 

1.039 

0.089 

0.243 

0.177 

17.78 

106.0 

0.025 

0.093 

372.0 

0.33 

0.35 

1.039 

0.070 

0.263 

C.  173 

17.40 

104.0 

0.017 

0.096 

564.7 

0.32 

0.33 

1.039 

0.050 

0.284 

0.170 

17.02 

110.0 

0.009 

0.099 

1100.0 

0.30 

0.32 

1.039 

0.027 

0.307 

0.  166 

16.65 

112.0 

0.001 

0.103 

10300.0 

0.29 

0.30 

1.039 

0.003 

0.335 

0.163 

16.28 

114.0 

-0.667 

0.220 

-32.9 

0.28 

0.29 

1.039 

-2.277 

0.751 

0.  160 

15.91 

116.0 

-0.644 

0.199 

-30.9 

0*26 

0.27 

1.039 

-2.307 

0.712 

0.156 

15.55 

118.0 

-0.636 

0.171 

-26.8 

0.25 

0.26 

1.039 

-2.392 

0.643 

0.153 

15.19 

120.0 

-0.645 

0.164 

-25.4 

0*24 

0.25 

1.039 

-2.548 

0.648 

0.149 

14.84 

••  3 STANDARD  DEVIATIONS 

TABLE  NO.  18 

DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  3 - CAMERA  DATA 


UNCLASSIFIED 


mM  copy  FUKHISHSD  VO 


UNClASSIFIED 


STP  453 


*10  4 9.3  IN  01  AM  6,7  PSI  9*3  IN  DIAM  ENO  PLATES 

PFAK  OVFRPRFSSURP  6.7  PS!  POSITIVE  DURATION  249.0  MSEC 

FRIEDLANOFR  OFF  AY  CONSTANT  • 0.31 


TIME  AFTFR 

DRAG 

ARSOLUTE 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

FLOW 

REYNOLDS 

SHOCK  FRONT 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PD 

0 

01 

CO 

NO. 

R 

(MSEC) 

(PSI  ) 

(PSI  ) 

(PERCENT) 

(PSI  ) 

(PSI  ) 

(Ol/Q) 

(PD/QI ) 

M 

U10-5I 

2.0 

0.390 

0.192 

49.3 

1.16 

1.19 

1.022 

0.326 

0.161 

0.296 

16.53 

4.0 

0.376 

0.175 

46.6 

1.14 

1.16 

1.C22 

0.322 

0.150 

0.293 

16.34 

6.0 

0.361 

0.159 

44.0 

1.11 

1.13 

1.022 

0.317 

0.139 

0.290 

16.16 

0.0 

0.347 

0.143 

41.1 

1.06 

1.11 

1.022 

0.313 

0.128 

0.287 

15.97 

10.0 

0.334 

0.129 

38.4 

1.05 

1.08 

1.022 

0.309 

0.118 

C.284 

15.79 

12.0 

0.321 

0.115 

35.7 

1.03 

1.05 

1.022 

0.304 

0.106 

0.281 

15.61 

14.0 

0.309 

0.102 

33.1 

1.00 

1.03 

1.022 

0.300 

0.099 

0.278 

15.43 

16.0 

0.297 

0.091 

30.7 

0.98 

1.00 

1.022 

0.295 

0.090 

0.275 

15.25 

10.0 

0.295 

o.oeo 

29.0 

0.95 

0.97 

1.022 

0.291 

0.081 

0.272 

15.07 

20.0 

0.274 

0.071 

25.9 

0.93 

0.95 

1.022 

0.287 

0.074 

0.269 

14.89 

22.0 

0.263 

0.062 

23.5 

0.90 

0.92 

1.022 

0.263 

0.066 

0.266 

14.71 

24.0 

0.253 

0.054 

21.6 

0.88 

0.90 

1.022 

0.279 

0.060 

0.263 

14.54 

26.0 

0.243 

0.04B 

20.0 

0.86 

0.88 

1.0  22 

0.275 

C.055 

0.261 

14.36 

2(1.0 

0.233 

0.043 

18.7 

0.84 

0.85 

1.022 

0.271 

C.050 

0.258 

14.19 

30.0 

0.223 

0.040 

17.9 

0.81 

0.83 

1.022 

0.267 

0.047 

0.255 

14.01 

32.0 

0.215 

0.037 

17.4 

0.79 

0.81 

1.022 

0.264 

0.046 

0.252 

13.84 

34.0 

0.206 

0.036 

17.6 

0.77 

0.79 

1.022 

0.260 

0.045 

0.249 

13.67 

36.0 

0.199 

0.036 

18.2 

0.73 

0.77 

1.022 

0.257 

0.046 

0.246 

13.50 

30.0 

0.191 

0.036 

18.8 

0.73 

0.75 

1.022 

0.253 

0.047 

0.243 

13.33 

40.0 

0.193 

0.036 

19.9 

0.71 

0.73 

1.022 

0.251 

0.050 

0.241 

13.16 

42.0 

0.177 

0.037 

21.1 

0.69 

0.71 

1.022 

0.246 

0.052 

0.236 

12.99 

44.0 

0.170 

0.037 

22.1 

0.67 

0.69 

1.022 

0.246 

0.054 

0.235 

12.83 

46.0 

0.164 

0.039 

23.5 

0.65 

0.67 

1*022 

0.243 

0.057 

0.232 

12.66 

4(1.0 

0.139 

0.039 

24.0 

0.64 

0.65 

1.022 

0.241 

0.060 

0.229 

12.49 

90.0 

0.152 

0.039 

26.0 

0.62 

0.63 

1.022 

0.238 

0.062 

0.227 

12.33 

92.0 

0.146 

0.039 

27.0 

0.60 

0.61 

1.022 

0.237 

0.064 

0.224 

12.17 

94.0 

0.141 

0.039 

27.7 

0.58 

0.60 

1.022 

0.235 

0.065 

0.221 

12.00 

96.0 

0.137 

0.039 

28.4 

0.57 

0.58 

1.022 

0.234 

0.066 

0.218 

11.84 

99.0 

0.131 

0.03* 

?9.2 

0.55 

0.56 

1.022 

0.232 

0.067 

0.216 

11.68 

AO.  0 

0.126 

0.037 

29.7 

0.53 

0.55 

1.022 

0.230 

0.068 

0.213 

11.52 

62.0 

0.123 

0.036 

29.6 

0.52 

0.53 

1.022 

0.230 

0.068 

0.210 

11.36 

A4.0 

0.119 

0.033 

29.7 

0.50 

0.51 

1.022 

0.229 

0.068 

0.207 

11.20 

AA.O 

0.114 

0.034 

29.6 

0.49 

0.50 

1.022 

0.227 

0.067 

0.205 

11.05 

A0.O 

0.111 

0.033 

29.6 

0.47 

0.48 

1.022 

0.228 

0.067 

0.202 

10.89 

70.0 

0.1C9 

0.031 

29.2 

0.46 

0.47 

1.022 

0.227 

0.066 

0.  199 

10.73 

72.0 

0.105 

0.030 

28.6 

0.44 

0.43 

1.022 

0.228 

0.065 

0.197 

10.58 

74.0 

0.101 

0.029 

20.6 

0.43 

0.44 

1.022 

0.227 

0.065 

0.194 

10.42 

76.0 

0.090 

0.027 

27.9 

0.42 

0.43 

1.022 

0.228 

0.063 

0.191 

10.27 

78.0 

0.093 

0.026 

27.8 

0.40 

0.41 

1*022 

0.228 

0.063 

0.  109 

10.12 

*0.0 

0.092 

0.023 

27.7 

0.39 

0.40 

1.022 

0.228 

0.063 

0.166 

9.97 

82.0 

C.  009 

0.024 

27.7 

0.38 

0.39 

1.022 

0.227 

0.063 

0.184 

9.82 

84.0 

0.0P6 

0.024 

28.3 

0.37 

0.38 

1.022 

0.228 

0.064 

0.181 

9,67 

86.0 

0.083 

0.024 

29.1 

0.36 

0.36 

1.022 

0.227 

0.066 

C.178 

9.52 

80.0 

0.080 

0.024 

31.0 

0.34 

0.35 

1.022 

0.225 

0.070 

0.176 

9.37 

90.0 

0.077 

0.025 

32.8 

0.33 

0.34 

1.022 

0.225 

0.073 

0.173 

9.22 

92.0 

0.074 

0.026 

35.0 

0.32 

0.33 

1.022 

0.224 

0.078 

0.171 

9.08 

94.0 

0.072 

0.026 

37.2 

0.31 

0.32 

1.022 

0.224 

0.083 

0.168 

8.93 

96.0 

0.060 

0.028 

40.9 

0.30 

0.31 

1.022 

0.220 

0.090 

0.165 

8.79 

98.0 

0.066 

C.029 

43.6 

0.29 

0.30 

1.022 

0.'221 

0.096 

0.163 

8.64 

100.0 

0.063 

0.029 

46.7 

0.28 

0.29 

1.022 

0.219 

0.102 

0.160 

8.50 

102.0 

0.061 

0.030 

49.8 

0.27 

0.28 

1.022 

0.220 

0.109 

0.158 

8.36 

104.0 

0.039 

0.032 

54.1 

0.26 

0.27 

1.022 

0.219 

0.119 

0.155 

8.22 

106.0 

0.036 

0.032 

56.9 

0.23 

0.26 

1.022 

0.218 

0.124 

0.153 

8.07 

100.0 

0.033 

0.033 

60.4 

0.24 

0.25 

1 .022 

0.219 

0.132 

0.150 

7.94 

110.0 

0.032 

0.034 

65.0 

0.23 

0.24 

1.022 

0.216 

0.140 

0.148 

7.80 

112.0 

0.049 

0.034 

69.2 

0.22 

0.23 

1.022 

0.212 

0.146 

0.145 

7.66 

114.0 

0.047 

0.034 

72.4 

0.22 

0.22 

1.022 

0.211 

0.153 

0.143 

7.52 

11A.0 

0.044 

0.034 

75.8 

0.21 

0.21 

1.02? 

0.207 

0.157 

0.140 

7.38 

118.0 

0.041 

0.034 

81.1 

0.20 

0.20 

1.022 

0.201 

0.163 

0.138 

7.25 

120.0 

".039 

0.033 

84.3 

0.19 

0.20 

1.022 

0.197 

0.167 

0.136 

7.11 

122.0 

0.036 

0.032 

90.5 

0.18 

0.19 

1.022 

0.188 

0.170 

0.133 

6.98 

124.0 

0.033 

0.032 

97.3 

0.10 

0.18 

1.022 

0.178 

0.173 

0.131 

6.85 

126.0 

0.029 

0.032 

107.3 

0.17 

0.17 

1.02? 

0.167 

0.180 

0.128 

6.71 

128.0 

0.023 

0.C31 

121.1 

0.16 

0.17 

1.02? 

0.152 

0.184 

0.126 

6.56 

130.0 

0.021 

0.031 

144.6 

0.16 

0.16 

1.022 

0.134 

0.194 

0.124 

6.45 

132.0 

0.C17 

0.032 

183.2 

0.15 

0.15 

1.022 

0.112 

0.206 

0.121 

6.32 

134.0 

0.013 

0.033 

239.3 

0.14 

0.15 

1.022 

0.092 

0.222 

0.119 

6.19 

136.0 

0.010 

0.033 

346.7 

0.14 

0.14 

1.022 

0.071 

0.246 

0.116 

6.06 

130.0 

0.007 

0.038 

511.2 

0.13 

0.13 

1.022 

0.054 

0.280 

0.114 

5.94 

140.0 

0.CD4 

0.042 

672.3 

0.12 

0.13 

1.022 

0.037 

0.323 

0.112 

5.81 

142.0 

0.001 

0.047 

2390.7 

0.12 

0.12 

1.022 

0.015 

0.375 

0.109 

5.66 

144.0 

-c.coo 

0.C33 

-13129.2 

0.11 

0.12 

1.022 

-0.002 

0.446 

0.1C7 

5.56 

146.0 

•0.003 

C.060 

-1743.3 

0.11 

0.11 

1.022 

-0.030 

0.526 

0.105 

5.43 

148.0 

-0.007 

C.C68 

-916.9 

0.10 

0.10 

1.022 

—0.060- 

0.626 

0.102 

5.31 

130.0 

-0.011 

0.077 

-705.8 

0.10 

0.10 

1.022 

—0.103 

0.745 

0.100 

5.19 

132.0 

-0.C13 

0.088 

-304.6 

0.09 

0.09 

1.022 

-0.151 

0.886 

0.098 

5.06 

134.0 

-0.C19 

0.098 

-304.1 

0.09 

0.09 

1.022 

-0.207 

1.047 

0.096 

4.94 

136.0 

-0.024 

0.110 

-430.3 

o.oe 

0.08 

1.022 

-0.269 

1.236 

0.093 

4.82 

138.0 

-0.030 

0.123 

-400.3 

0.08 

0.08 

1.022 

-0.355 

1.454 

0.091 

4.70 

3 STANDARD  DEVIATIONS 

TABLE  WO.  19 


DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  4 - VELOCITY  TRANSDUCER  DATA 


UNCLASSIFIED 


; yr~ 


NO  4 9.5  IN  0 1 AM  6.7  PSI 

PEAK  OVFRPRESSURE  6.7  PSI 
FRIEDLANDFR  DECAY  CONSTANT  ■ 0.51 


9.5  IN  DIAM  END  PLATES 
POSITIVE  DURATION  2<»9.0  MSEC 


STP  453 


TIME  AFTER 

DRAG 

ABSOLUTE 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

FLOW 

REYNOLDS 

SHOCK  FRONT 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PD 

0 

01 

CD 

NO. 

R 

(MSEC) 

(PSI  ) 

(PSI  ) 

(PERCENT) 

(PSD 

(PSI  ) 

(QI/O) 

(PD/QI ) 

M 

(X10-5) 

2.0 

0.49  8 

0.067 

13.4 

1.16 

1.19 

1.022 

0.416 

0.056 

0.296 

16.53 

4.0 

0.482 

0.060 

12*4 

1.14 

1.16 

1.022 

0.413 

0.051 

0.293 

16.34 

6.0 

0.465 

0.052 

11.1 

1.11 

1.13 

1.022 

0.408 

0.045 

0.290 

16.16 

8.0 

0.449 

0.046 

10.2 

1.08 

1.11 

1.022 

0.404 

0.041 

0.287 

15.97 

10.0 

0.434 

0.040 

9.2 

1.05 

1.08 

1.022 

0.400 

0.036 

0.284 

15.79 

12.0 

0.418 

0.034 

8.1 

1.03 

1.05 

1.022 

0.395 

0.032 

0.261 

15.61 

14.0 

0.403 

0.029 

7.1 

1.00 

1.03 

1.022 

0.391 

0.026 

0.278 

15.43 

16.0 

0.388 

0.024 

6.1 

0.98 

1.00 

1.022 

0.386 

0.023 

0.275 

15.25 

18.0 

0.374 

0.020 

5.3 

0.95 

0.97 

1.022 

0.382 

0.020 

0.272 

15.07 

20.0 

0.360 

0.017 

4.7 

0.93 

0.95 

1.022 

0.377 

0.017 

0.269 

14.89 

22.0 

0.346 

0.013 

3.7 

0.90 

0.92 

1.022 

0.372 

0.013 

0.266 

14.71 

24.0 

0.332 

0.011 

3.3 

0.68 

0.90 

1.022 

0.366 

0.C12 

0.263 

14.54 

26.0 

0.319 

0.009 

2.8 

0.86 

0.88 

1.022 

0.361 

0.010 

0.261 

14.36 

28.0 

0.306 

0.007 

2.2 

0.84 

0.85 

1.022 

0.355 

0.008 

0.256 

14.19 

30.0 

0.293 

0.006 

2.0 

0.81 

0.83 

1.022 

0.349 

0.007 

0.255 

14.01 

32.0 

0.281 

0.006 

2.1 

0.79 

0.81 

1.022 

0.344 

0.007 

0.252 

13.84 

34.0 

0.268 

0.006 

2.2 

0.77 

0.79 

1*022 

0.337 

0.007 

0.249 

13.67 

36.0 

0.25* 

0.007 

2.7 

0.75 

0.77 

1*022 

0.332 

0.009 

0.246 

13.50 

38.0 

0.245 

0.007 

2.8 

0.73 

0.75 

1.022 

0.325 

0.009 

0.243 

13.33 

40.0 

0.234 

0.007 

2.9 

0.71 

0.73 

1.022 

0.319 

0.009 

0.241 

13.16 

42.0 

0.223 

0.008 

3.5 

0.69 

0.71 

1.022 

0.313 

0.011 

0.238 

12.99 

44.0 

0.212 

0.008 

3.7 

0.67 

0.69 

1.025 

0.306 

0.011 

0.235 

12.63 

46.0 

0.202 

0.008 

3.9 

0.65 

0.67 

1.022 

0.299 

0.011 

0.232 

12.66 

46.0 

0.192 

0.008 

4.1 

0.64 

0.65 

1.022 

0.293 

0.C12 

0.229 

12.49 

50.0 

0.182 

0.007 

3.8 

0.62 

0.63 

1.022 

0.285 

0.010 

0.227 

12.33 

52.0 

0.173 

0.007 

4.0 

0.60 

0.61 

1.022 

0.279 

0.011 

0.224 

12.17 

54.0 

0.164 

0.006 

3.6 

0.58 

0.60 

1.022 

0.272 

0.009 

0.221 

12.00 

56.0 

0.155 

0.006 

3.8 

0.57 

0.58 

1.022 

0.265 

0.010 

0.218 

11.84 

58.0 

0.147 

0.005 

3.4 

0.55 

0.56 

1.022 

0.259 

0*008 

0.216 

11.66 

60.0 

0.139 

0.005 

3.5 

0.53 

0.55 

1.022 

0.252 

0.009 

0.213 

11.52 

62.0 

0.131 

0.005 

3.8 

0.52 

0.53 

1.022 

0.244 

0.009 

0.210 

11.36 

64.0 

0.123 

0.004 

3.2 

0.50 

0.51 

1.022 

0.236 

0.007 

0.207 

11.20 

66.0 

0.116 

0.004 

3.4 

0.49 

0.50 

1.022 

0.230 

0.007 

0.205 

11.05 

68.0 

0.109 

0.004 

3.6 

0.47 

0.48 

1.022 

0.222 

0.008 

£ .202 

10.69 

70.0 

0.102 

0.C04 

3.9 

0.46 

0.47 

1.022 

0.215 

0.008 

0.199 

10.73 

72.0 

0.096 

0.005 

5.2 

0.44 

0.45 

1*022 

0.208 

0.C10 

0.197 

10.58 

74.0 

0.090 

0.005 

5.5 

0.43 

0.44 

1*022 

0.201 

0.011 

0.194 

10.42 

76.0 

0.084 

0.006 

7.1 

0.42 

0.43 

1.022 

0.194 

0.013 

0.191 

10.27 

78.0 

C.079 

0.006 

7.5 

0.40 

0.41 

1.022 

0.168 

0.014 

0.189 

10.12 

80.0 

0.073 

0.007 

9.5 

0.39 

0.40 

1.022 

0.179 

0.017 

0.186 

9.97 

82.0 

0.069 

0.007 

10.1 

0.38 

0.39 

1.022 

0.175 

0.017 

0.184 

9.82 

84.0 

0.064 

0.007 

10.9 

0.37 

0.38 

1.022 

0.166 

0.018 

0.181 

9.67 

86.0 

0.060 

0.008 

13.3 

0.36 

0.36 

1.022 

0.162 

0.021 

0.178 

9.52 

88.0 

0.056 

0.008 

14.2 

0.34 

0.35 

1.022 

0.157 

0.022 

0.176 

9.37 

90. C 

0.052 

0.008 

15.3 

0.33 

0.34 

1.022 

0.150 

0.023 

0.173 

9.22 

92.0 

0.049 

0.008 

16.3 

0.32 

0.33 

1.022 

0.147 

0.024 

0.171 

9.06 

94.0 

0.046 

0.007 

15.2 

0.31 

0.32 

1.022 

0.142 

0.021 

0.168 

8.93 

96.0 

0.043 

0.007 

16.2 

0.30 

0.31 

1.022 

0.136 

0.022 

0.165 

8.79 

98.0 

0.041 

0.007 

17.0 

0.29 

0.30 

1.022 

0.136 

0.023 

0.163 

6.64 

100.0 

0.038 

0.006 

15.7 

0.26 

0.29 

1.022 

0.1  SO 

0.020 

0.160 

8.50 

102.0 

0.037 

0.006 

16.2 

0.27 

0.28 

1.022 

0.132 

0.021 

0.158 

8.36 

104.0 

0.035 

0.007 

20.0 

0.26 

0.27 

1.022 

0.129 

0.025 

0.155 

8.22 

106.0 

0.034 

0.008 

23.5 

0.25 

0.26 

1.022 

0.130 

0.030 

0.153 

8.07 

108.0 

0. 033 

0.010 

30.3 

0.24 

0.25 

1.022 

0.131 

0.039 

0.150 

7,94 

110.0 

0.032 

0.012 

37.5 

0.23 

0.24 

1.022 

0.131 

0.049 

0.148 

7.80 

112.0 

0.032 

0.015 

46.8 

0.22 

0.23 

1.022 

0.136 

0.064 

0.145 

7.66 

114.0 

0.032 

0.018 

56.2 

0.22 

0.22 

1.022 

0.142 

0.080 

0.143 

7.52 

116.0 

0.032 

0.022 

68.7 

0.21 

0.21 

1.022 

0.147 

0.101 

0.140 

7.38 

118.0 

0.033 

0.027 

81.8 

0.20 

0.20 

1.022 

0.156 

0.129 

0.138 

7.25 

120.0 

0.033 

0.031 

93.9 

0.19 

0.20 

1.022 

0.164 

0.154 

0.236 

7.11 

122.0 

0.035 

0.037 

105.7 

0.16 

0.19 

1.022 

0.161 

0.192 

0.133 

6.98 

124.0 

0.036 

0.043 

119.4 

0.16 

0.18 

1.022 

0.194 

0.232 

0.131 

6,65 

126.0 

0.038 

0.049 

128.9 

0.17 

0.17 

1.022 

0.213 

0.275 

0.128 

6.71 

128.0 

0.040 

0.056 

140.0 

0.16 

0.17 

1.022 

0.234 

0.328 

0.126 

6.56 

130.0 

0.042 

0.063 

150.0 

0.16 

0.16 

1.022 

0.256 

0.365 

0.124 

6.45 

132.0 

0.045 

0.071 

157.7 

0.15 

0.15 

1.022 

0.266 

0.452 

0.121 

6.32 

•*  i standard  deviations 

TABLE  NO.  20 

DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  4 - CAMERA  DATA 
UNCLASSIFIED 


SET : 

— - - — 

■fci  Ar " - 

THIS  PACE  IS  BEST  QUALITY  PBtWSHCAMi 

7B0«  COPY  FURBISHED  TODDC  


UNCLASSIFIED 


STP  453 


NO  A 9.*  IN  DIAM  6,7  PSI 
PEAK  OVERPRESSURE  6.7  PSI 
FRIFOLANOPR  OFCAY  CONSTANT  • 0.51 


28.5  IN  0 1 AM  END  PLATES 
POSITIVE  DURATION  249.0  MSEC 


TIMF  AcTFp 

DRAG 

ARSOLUTF 

RELATIVE 

DYNAMIC 

impact 

RATIO 

DRAG 

ABSOLUTE 

FLOW 

REYNOLDS 

SHOCK  FRONT 

PRESSURE 

ERROR** 

error** 

PRESSURE 

PRESSURE 

coefficient 

ERROR** 

MACH 

NO. 

ARRIVAL 

PD 

0 

01 

CD 

NO. 

R 

(MSEC) 

(PSI  ) 

(PSI  ) 

(PERCENT) 

(PSI  ) 

(PSI  ) 

(01/0) 

(PD/OI ) 

M 

(X10-5) 

2.0 

0.527 

C.  139 

26.3 

1.16 

1.19 

1.022 

0.441 

0.116 

0.296 

16.53 

4.0 

o.519 

0.134 

25.8 

1.14 

1.16 

1.022 

0.444 

0.115 

0.293 

16.34 

6.0 

0.5C9 

0.129 

25.3 

1.11 

1.13 

1.022 

0.447 

0.113 

0. 290 

16.16 

P.0 

0.500 

C.  124 

24.9 

1.08 

1.11 

1.022 

0.451 

0.112 

0.287 

15.97 

10.0 

0.491 

0.120 

24.4 

1.05 

1.08 

1.022 

0.453 

0.110 

0.284 

15.79 

12.0 

0.482 

0.115 

23.9 

1.03 

1.05 

1.022 

0.456 

0.109 

0.281 

15.61 

14.0 

0.473 

o.no 

23.3 

1.00 

1.03 

1.02? 

0.460 

0.107 

0.278 

15.43 

1 A.  0 

0.465 

0.106 

22.8 

C .98 

1.00 

1.022 

0.463 

0.105 

0.275 

15.25 

1*».0 

0.455 

C.101 

22.2 

0.95 

0.97 

1.022 

0.465 

0.103 

0.272 

15.07 

20.0 

0.447 

0.097 

21.7 

0.93 

0.95 

1.02  2 

0.466 

0.101 

0.269 

14.89 

22.0 

0.438 

0.092 

21.0 

0.90 

0.92 

1.022 

0.471 

0.099 

0.266 

14.71 

24.0 

0.429 

0.088 

20.5 

0.88 

0.90 

1.02? 

0.474 

0.097 

0.263 

14.54 

26.0 

0.420 

0.083 

19.7 

0.86 

0.88 

1.022 

0.476 

0.094 

0.261 

14.36 

?®.o 

0.411 

0.078 

19.1 

0.84 

0.85 

1.022 

0.478 

0.091 

0.258 

14.19 

70.0 

C.  403 

0.074 

18.5 

0.81 

0.C3 

1.02? 

0.481 

0.089 

0.255 

14.01 

72.0 

0.393 

0.070 

17.8 

0.79 

0.81 

1.022 

0.483 

0.086 

0.252 

13.84 

74.0 

0.385 

0.066 

17.1 

0.77 

0.79 

1.022 

0.485 

0.083 

0.249 

13.67 

•»A.O 

0.376 

0.062 

16.5 

0.75 

0.77 

1.022 

0.486 

0.080 

0.246 

13.50 

7«.0 

0.367 

0.058 

15.9 

0.73 

0.75 

1.02? 

0.488 

0.077 

0.243 

13.33 

40.0 

0.357 

0.054 

15.2 

0.71 

0.73 

1.022 

0.489 

0.074 

0.241 

13.16 

42.0 

0.348 

0.051 

14.6 

0.69 

0.71 

1.02? 

0.490 

0.071 

0.238 

12.99 

44.0 

0.338 

0.048 

14.3 

0.67 

0.69 

1.022 

0.489 

C.070 

0.235 

12.83 

4 A , 0 

0.328 

0.045 

13.9 

0.65 

0.67 

1.022 

0.487 

0.066 

0.232 

12.66 

A 8 . 0 

0.317 

0.043 

13.8 

0.64 

0.65 

1.022 

0.484 

0.066 

0.229 

12.49 

50.0 

0.306 

0.042 

13.8 

0.62 

0.63 

1.022 

0.480 

0.066 

0.227 

12.33 

52.0 

0.294 

0.041 

14.0 

0.60 

0.61 

1.022 

0.475 

0.066 

0.224 

12.17 

54.0 

0.282 

0.041 

14.5 

0.58 

0.60 

1.022 

0.469 

0.068 

0.221 

12.00 

56.0 

0.270 

0.040 

15.1 

0.57 

0.58 

1 .022 

0.462 

0.070 

0.218 

11.84 

58.0 

0.258 

0.041 

15.9 

0.55 

0.56 

1.022 

0.455 

0.072 

0.216 

11.66 

50.0 

0.247 

0.04? 

17.2 

0.53 

0.55 

1.02? 

0.449 

0.077 

0.213 

11.52 

A2.0 

0.237 

0.043 

18.5 

0.52 

0.53 

1.022 

0.443 

0.062 

0.210 

11.36 

64,0 

0.226 

0.045 

20.1 

0.50 

0.51 

1.022 

0.435 

0.087 

0.207 

11.20 

66.0 

C.217 

0.048 

22.1 

0.49 

0.50 

1.022 

0.431 

0.095 

0.205 

11.05 

68.0 

0.207 

0.050 

24.4 

0.47 

0.48 

1.022 

0.424 

0.103 

0.202 

10.89 

70.0 

0.198 

0.053 

27.0 

0.46 

0.47 

1.02? 

0.419 

0.113 

0.199 

10.73 

72.0 

0.189 

0.056 

29.6 

0.44 

0.45 

1.022 

0.411 

0.122 

0.197 

10.58 

74.0 

0.181 

0.060 

33.1 

0.43 

0.44 

1.022 

0.406 

0.134 

0.194 

10.42 

76.0 

0.172 

0.063 

36.6 

0.42 

0.43 

1.022 

0.399 

0.146 

0.191 

10.27 

78.0 

0.163 

0.067 

41.0 

0.40 

0.41 

1.02? 

0.390 

0.160 

0.189 

10.12 

80.0 

0.155 

0.070 

45.5 

0.39 

0.40 

1.C2? 

0.382 

0.174 

0.186 

9. 97 

82.0 

0.146 

0.074 

50.7 

0.38 

0.39 

1.022 

0.372 

0.186 

0.184 

9.8? 

84.0 

0.137 

0.078 

56.7 

P • 37 

0.38 

1.0  22 

0.361 

0.205 

C.  181 

9.67 

8 6.0 

0.128 

0.081 

63.5 

0.36 

0.36 

1.022 

0.350 

0.222 

0.170 

9.52 

88.0 

0.120 

0.085 

71.5 

0.34 

0.35 

1.02? 

0.337 

0.241 

0.176 

9.37 

90. 0 

0.111 

0.089 

80.6 

0.33 

0.34 

1.022 

0.323 

0.260 

0.173 

9.22 

92.0 

0.102 

0.093 

91.6 

0.32 

0.33 

1.022 

0.306 

0.281 

0.171 

9.08 

94.0 

0.093 

0.097 

104.1 

0.31 

0.32 

1.02? 

0.291 

0.303 

0.168 

8.93 

«6.0 

0.084 

c.ioi 

119.4 

0.30 

0.31 

1.022 

0.273 

0.326 

0.165 

6.79 

**  * standard  deviations 


TABLE  NO.  21 


DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  5 - VELOCITY  TRANSDUCER  DATA 


UNCLASSIFIED 


*1 


1S  0BS1  QlttWff 


UNCLASSIFIED 


STP  453 


NO  5 9,5 

IN  01 AM 

6.7 

PSI 

28.5 

IN  DIAM 

END  PLATES 

PEAK  OVERPRESSURE 

6.7 

PSI 

POSITIVE 

DURATION 

249.0  MSEC 

FRT  EDLANDER 

DECAY  CONSTANT  • 0.51 

TIME  AFTER 

drag 

ABSOLUTE 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

FLOW 

Rl-NOLDS 

SHOCK  FRONT 

PRESSURE 

ERROR 

ERROR  1 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR 

MACH 

NO. 

ARRIVAL 

PD 

0 

01 

CD 

NO. 

R 

(MSEC) 

(PSD 

(PSI  ) 

(PERCENT) 

(PSI  ) 

(PSI  ) 

(01/0) 

(PD/QI ) 

M 

(X10-5) 

2.0 

0.677 

0.280 

41.3 

1.16 

1.19 

1.022 

0.566 

0.234 

0.296 

16.53 

4.0 

0.657 

0.247 

37.5 

1*14 

1.16 

1.022 

0.563 

0.211 

0.293 

16.34 

6*0 

0.637 

0.216 

33.9 

1.11 

1*13 

1.022 

0.559 

0.189 

0.290 

16.16 

8.0 

0.618 

0.188 

30.4 

1.08 

1.11 

1.022 

0.556 

0.169 

0.287 

15.97 

10.0 

0.600 

0.162 

26.9 

1.05 

1.08 

1.022 

0.553 

0.149 

0.284 

15.79 

12.0 

0.581 

0.138 

23.7 

1.03 

1.05 

1.022 

0.549 

0.130 

0.261 

15.61 

14.0 

0.564 

0.116 

20.5 

1.00 

1.03 

1.022 

0.547 

0.112 

0.278 

15.43 

16.0 

0.547 

0.097 

17.7 

0.98 

1.00 

1.022 

0.544 

0.096 

0.275 

15.25 

18.0 

0.530 

0.080 

15.0 

0.95 

0.97 

1.022 

0.541 

0.081 

0.272 

15.07 

20.0 

0.514 

0.064 

12.4 

0.93 

0.95 

1.02  2 

0.538 

0.067 

0.269 

14.69 

22.0 

0.499 

0.052 

10.4 

0.90 

0.92 

1.022 

0.536 

0.055 

0.266 

14.71 

24.0 

0.484 

0.041 

8.4 

0.68 

0.90 

1.022 

0.534 

0.045 

0.263 

14.54 

26.0 

0.469 

0.034 

7.2 

0.86 

0.88 

1.022 

0.531 

0.038 

0.261 

14.36 

24.0 

0.455 

0.029 

6.3 

0.84 

0.85 

1.022 

0.529 

0.033 

0.258 

14.19 

40.0 

0.441 

0.028 

6.3 

0.81 

0.83 

1.022 

0.526 

0.033 

0.255 

14.01 

32.0 

0.427 

0.028 

6.5 

0.79 

0.81 

1.022 

0.523 

0.034 

0.252 

13.84 

34.C 

0.414 

0.029 

7.0 

0.77 

0.79 

1.022 

0.521 

0.036 

0.249 

13.67 

36.0 

0.401 

0.031 

7.7 

0.75 

0.77 

1.022 

0.518 

0.040 

0.246 

13.50 

38.0 

0.389 

0.032 

8.2 

0.73 

0.75 

1.022 

0.517 

0.042 

0.243 

13.33 

40.0 

0.377 

0.033 

8.7 

0.71 

0.73 

1.022 

0.515 

0.045 

0.241 

13.16 

4 2.0 

0.365 

0.034 

9.3 

0.69 

0.71 

1.022 

0.512 

0.047 

0.238 

12.99 

44.0 

0.354 

0.034 

9.6 

0.67 

0.69 

1.022 

0.511 

0.049 

0.235 

12.83 

46.0 

0.343 

0.033 

9.6 

0.65 

0.67 

1.022 

0.509 

0.046 

0.232 

12.66 

48.0 

0.332 

0.032 

9.6 

0.64 

0.65 

1.022 

0.506 

0.046 

0.229 

12.49 

50.C 

0.322 

0.030 

9.3 

0.62 

0.63 

1.022 

0.5C5 

0.047 

0.227 

12.33 

52.0 

0.312 

0.028 

8.9 

0.60 

0.61 

1.022 

0.504 

0.045 

0.224 

12.17 

54.0 

0.302 

0.026 

8.6 

0.58 

0.60 

1.022 

0.502 

0.043 

0.221 

12.00 

56.0 

C.292 

0.024 

8.2 

0.57 

0.58 

1.022 

0.499 

0.041 

0.218 

11.84 

58.0 

0.282 

0.022 

7.e 

0.55 

0.56 

1.022 

0.496 

0.038 

0.216 

11.68 

60.0 

0.273 

0.020 

7.3 

0.53 

0.55 

1.022 

0.495 

0.036 

0.213 

11.52 

62.0 

0.264 

0.019 

7.1 

0.52 

0.53 

1.022 

0.493 

0.035 

0.210 

11.36 

64.0 

0.255 

0.019 

7.4 

0.50 

0.51 

1.022 

0.491 

0.036 

0.207 

11.20 

66.0 

0.246 

0.019 

7.7 

0.49 

0.50 

1.022 

0.486 

0.037 

0.205 

11.05 

68.0 

0.237 

0.021 

8.8 

0.47 

0.48 

1.022 

0.484 

0.042 

0.202 

10.89 

70.0 

0.229 

0.02? 

9.6 

0.46 

0.47 

1.022 

0.482 

0.046 

0.199 

10.73 

72. C 

0.220 

0.025 

11.3 

0.44 

0.45 

1.022 

0.478 

0.054 

0.197 

10.58 

74.0 

0.212 

0.027 

12.7 

0.43 

0.44 

1.02? 

0.475 

0.060 

0.194 

10.42 

76.0 

0.204 

0.029 

14.2 

0.42 

0.43 

1.022 

0.472 

0.067 

0.191 

10.27 

78.0 

0.196 

0.031 

15.8 

0.40 

C .4 1 

1.022 

0.468 

0.074 

0.1H9 

10.12 

80.0 

0.188 

0.032 

17.0 

0.39 

0.40 

1.02? 

0.463 

0.076 

0.166 

9.97 

82.0 

0.180 

0.033 

18.3 

0.36 

0.39 

1.022 

0.458 

0.083 

0.184 

9.82 

84.0 

0.172 

0.034 

19.7 

0.37 

0.38 

1.022 

0.452 

0.089 

0.181 

9.67 

86.0 

0.164 

C.034 

20.7 

0.36 

0.36 

1.022 

0.445 

0.092 

0.178 

9.52 

88.0 

0.156 

0.033 

21.1 

0.34 

0.35 

1.022 

0.437 

0.092 

0.176 

9.37 

90.0 

0.148 

0.032 

21.6 

0.33 

0.34 

1.022 

0.429 

0.092 

0.173 

9 .22 

92.0 

0.140 

0.030 

21.4 

0.32 

0.33 

1.022 

0.420 

0.090 

0.171 

9.08 

94.0 

0.132 

0.029 

21.9 

0.31 

0.32 

1.022 

C.409 

0.090 

0.168 

8.93 

96.0 

0.124 

0.C28 

22.5 

0.30 

0.31 

1.022 

0.398 

0.089 

0.165 

8.79 

98.0 

0.116 

0.028 

24.1 

0.29 

0.30 

1.022 

0.385 

0.093 

0.163 

8.64 

100.0 

0.108 

0.031 

28.7 

0.28 

0.29 

1.022 

0.371 

0.106 

0.160 

8.50 

102.0 

0.09'' 

0.036 

36.3 

0.27 

0.28 

1.02? 

0.353 

0.128 

0.158 

8.36 

104.0 

0.091 

0.045 

49.4 

0.26 

0.27 

1.022 

0.336 

0.166 

0.155 

8 .22 

106.0 

C.083 

0.056 

67.4 

0.25 

0.26 

1.02? 

0.318 

0.214 

0.153 

8.07 

108.0 

0.074 

0.070 

94.5 

0.24 

0.25 

1.022 

0.294 

0.276 

0.150 

7.94 

110.0 

C.  065 

0.086 

132.3 

0.23 

0.24 

1.022 

0.268 

0.354 

0.148 

7.80 

112.0 

C.056 

0.104 

185.7 

0.2  2 

0.23 

1.022 

0.239 

0.445 

0.145 

7.66 

114.0 

0.047 

0.124 

263.8 

0.22 

0.22 

1.022 

0.208 

0.551 

0.143 

7.52 

116.C 

C.C38 

0.146 

384.2 

0.21 

0.21 

1.022 

0.175 

0.674 

0.140 

7.38 

119.0 

0.028 

0.171 

610.7 

0.20 

0.20 

1.022 

0.134 

0.820 

0.138 

7.25 

120.0 

0.C19 

0.198 

1042.1 

0.19 

0.20 

1.022 

0.094 

0.988 

0.136 

7.11 

122.0 

0.CC9 

0.227 

2522.2 

0.18 

0.19 

1.022 

0.046 

1.178 

0.133 

6.98 

1 ?4. 0 

C.000 

0.259 

25.9 

C.  18 

0.18 

1.022 

0.000 

1.399 

0.131 

6.85 

TABLE  NO.  22 


DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  5 - CAMERA  DATA 
UNCLASSIFIED 


t 


UNCLASSIFIED 


STP  453 


NO  6 3.5 

IN  DIAM 

9.7  PSI 

3.5 

IN  0 I AM 

END  PLATES 

raoa 

0OPY  yuRttlSHED  Iv  uwv 

PEAK  OVERPRESSURE 

9.7  PSI 

POSITIVE 

DURATION 

233.0  MSEC 

FRI FOLANPPP 

OFrAY  CONSTANT  ■ 0.83 

time  aft f° 

ORAG 

arsolutf 

RELATIVE 

dynamic 

impact 

RATIO 

DRAG 

ABSOLUTE 

FLO* 

REYNOLDS 

Shock  front 

PRFSSURF  ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PO 

0 

01 

CD 

NO. 

R 

(MSEC) 

(PSD 

(PSI  ) 

(PERCENT) 

(PSI  ) 

(PSI  ) 

(01/0) 

(PD/QI ) 

M 

(X1C-5 ) 

?.o 

0.777 

0.184 

23.7 

2.35 

2.45 

1.039 

0.317 

0.075 

0.391 

8.29 

4*0 

0.746 

0.167 

22.3 

2*28 

2.37 

1.039 

0.314 

0.070 

U.  366 

8.17 

6*0 

0.717 

0.149 

20.9 

2. 21 

2.30 

1.039 

0.311 

0.065 

0.382 

8.06 

0.0 

0.660 

0. 1 34 

19.5 

2.14 

2.22 

1.039 

0. 3f  9 

0.06C 

0.377 

7.94 

10*0 

0.661 

0.120 

10.1 

2.07 

2.15 

1.039 

0. 306 

0.055 

0.372 

7.83 

12.0 

0.636 

0.106 

16.7 

2.01 

2.08 

1.039 

0.304 

0.051 

C.367 

7.72 

14.0 

0.61? 

0.094 

15.4 

1.94 

2.0  2 

1.C39 

C.302 

0.046 

0.363 

7.60 

16.0 

0.5«9 

0.063 

14.1 

1.88 

1.95 

1.039 

C.3C1 

0.042 

0.358 

7.49 

16.0 

0.567 

0.073 

12.9 

1.82 

1.89 

1.039 

0.299 

0.03b 

0.353 

7.30 

20.0 

0.346 

C.  064 

11.8 

1.76 

1.83 

1.039 

0.298 

0.035 

C.  349 

7-27 

22.0 

0.526 

0.057 

10.8 

1.70 

1.77 

1.039 

0.296 

C.032 

0.344 

7.17 

24.0 

0.506 

0.051 

10.1 

1.65 

1.71 

1.039 

0.295 

0.029 

0.339 

7.06 

26.0 

0.465 

0.046 

9.5 

1.59 

1.65 

1.039 

0.292 

0.026 

0.335 

6.95 

20.0 

0.466 

0.043 

9.3 

1.54 

1.60 

1.039 

0.290 

0.027 

0.330 

6.85 

30.0 

0.440 

0.041 

9.2 

1.49 

1.55 

1.039 

0.209 

0.026 

0.326 

6.74 

32.0 

0.430 

0.040 

9.4 

1.44 

1.49 

1.039 

0.207 

0.02  7 

0.321 

6.64 

34.0 

0.413 

0.040 

9.7 

1.39 

1.44 

1.039 

0.285 

0.027 

C.317 

6.53 

36.0 

0.309 

0.040 

10.1 

1.34 

1.39 

1 .039 

0.285 

0.026 

0.312 

6.43 

30.0 

0.3«6 

0.040 

10.5 

1.30 

1.35 

1.039 

0.205 

0.03C 

0.3C6 

6.33 

40.0 

0.373 

0.040 

10.0 

1.25 

1.30 

1.039 

0.206 

0.031 

0.304 

6.23 

42.0 

0.362 

0.041 

11.3 

1*21 

1.26 

1.039 

0.287 

0.032 

C.  299 

6.13 

44.0 

0.331 

0.041 

11.7 

1.17 

1.21 

1.039 

0.289 

0.033 

0.295 

6 • C 3 

46.0 

o.342 

0.041 

11.9 

1.12 

1.17 

1.039 

0.291 

0.034 

C.291 

5.93 

46. 0 

C.333 

0.040 

12.2 

1.09 

1.13 

1.039 

0.294 

0.036 

0.286 

5.64 

30.0 

0.325 

0.040 

12.3 

1.05 

1.09 

1.039 

0.298 

0.036 

0.282 

5.74 

32.0 

0.310 

0.039 

12.3 

1.01 

1.05 

1.039 

0.301 

C.037 

0.2/8 

5.65 

94.0 

0.310 

0.030 

12.3 

0.97 

1.01 

1.039 

0.305 

0.037 

0.274 

5.55 

36.0 

0.304 

0.037 

12.2 

0.94 

0.97 

1.039 

0.310 

0.038 

0.270 

5.46 

36.0 

0.297 

0.036 

12.1 

0.90 

0.94 

1.039 

0.315 

0.038 

0.265 

5.36 

60.0 

0.292 

0.035 

11.9 

0.87 

0.90 

1.039 

0.321 

C.038 

0.261 

5.27 

62.0 

0.207 

0.033 

11.6 

0.84 

0.87 

1.039 

0.327 

0.038 

0.257 

5.18 

64.0 

0.262 

0.032 

11.5 

0.81 

0.84 

1.039 

0.334 

0.038 

0.253 

5.09 

66.0 

0.276 

0.031 

11.2 

0.78 

0.81 

1.039 

0.340 

0.036 

0.249 

5.00 

60.0 

0.270 

0.030 

11.0 

0.75 

0.78 

1.039 

0.346 

0.036 

0.245 

4.91 

70. P 

0.266 

0.029 

11.0 

0.72 

0.75 

1.039 

0.353 

0.038 

0.241 

4.83 

72.0 

0.260 

0.029 

11.1 

0.69 

0.72 

1.039 

0.360 

0.040 

C.  237 

4.74 

74.0 

0.256 

0.020 

11.1 

0.66 

0.69 

1.039 

0.368 

0.041 

0.233 

4.65 

76.0 

0.251 

0.020 

11.3 

0.64 

0.66 

1.039 

0.375 

0.042 

0.229 

4.57 

70.0 

0.247 

0.028 

11.5 

0.61 

0.64 

1.039 

0.385 

0.044 

0.225 

4.48 

60.0 

0.243 

0.029 

11.9 

0.59 

0.61 

1.039 

0.396 

0.047 

0.221 

4.40 

62.0 

0.240 

0.029 

12.2 

0.56 

0.59 

1.039 

0.405 

0.049 

0.218 

4.32 

04.0 

0.237 

0.030 

12.6 

0.54 

0.56 

1.039 

0.417 

0.052 

0.214 

4.23 

• 6.0 

0.234 

0.031 

13.2 

0.52 

0.54 

1.039 

0.429 

0.056 

0.210 

4.15 

60.0 

0.231 

0.031 

13.7 

0.50 

0.52 

1.C39 

0.442 

0.060 

0.206 

4.07 

90.0 

0.229 

0.03? 

14.1 

0.48 

0.50 

1.039 

0.456 

0.064 

0.202 

3.99 

02.0 

0.226 

0.033 

14.8 

0.46 

0.48 

1.039 

0.471 

0.069 

0.199 

3.91 

04.0 

0.224 

0.034 

15.2 

0.44 

0.46 

1.039 

0.486 

0.074 

0.195 

3.83 

06.0 

0.221 

0.034 

15.6 

0.4  2 

0.44 

1.039 

0.502 

0.078 

0.191 

3.76 

90.0 

0.216 

0.033 

16.0 

0.40 

0.42 

1.039 

0.516 

0.062 

0.  188 

3.68 

100.0 

0.214 

0.034 

16.2 

0.30 

0.40 

1.039 

0.530 

0.066 

0.184 

3.60 

102.0 

0.211 

0.035 

16.6 

0.37 

0.38 

1.039 

0.545 

0.090 

0.181 

3.53 

104.0 

0.207 

0.034 

16.7 

0.35 

0.36 

1.039 

0.561 

0.094 

0.177 

3.45 

106.0 

0.203 

0.034 

16.9 

0.33 

0.35 

1.039 

0.576 

0.097 

0.173 

3.38 

100.0 

0.199 

0.034 

17.4 

0.32 

0.33 

1.039 

0.590 

0.103 

0.170 

3.31 

110.0 

0.193 

0.034 

17.9 

0.30 

0.32 

1.039 

0.601 

0.106 

0.  166 

3.23 

112.0 

0.1*0 

0.033 

18.6 

0.29 

0.30 

1.039 

0.614 

0.114 

0.163 

3.16 

114.0 

C.102 

0.035 

19.5 

0.28 

0.29 

1.039 

0.622 

0.121 

0.160 

3.09 

116.0 

0.176 

0.036 

20.3 

0.26 

0.27 

1.039 

0.633 

0.129 

0.156 

3.02 

110.0 

0.169 

0.037 

21.0 

0.25 

0.26 

1.039 

0.639 

0.139 

0.153 

2.95 

120.0 

0.163 

0.037 

22.7 

0.24 

0.25 

1.039 

0.646 

0.146 

0.149 

2.88 

122.0 

0.135 

0.030 

24.0 

0.23 

0.24 

1.039 

0.645 

0.160 

0.146 

2.81 

124.0 

0.147 

0.040 

27.1 

0.22 

0.22 

1.039 

0.644 

0.175 

0.143 

2.74 

126.0 

0.136 

0.042 

30.5 

0.20 

0.21 

1.039 

0.636 

0.194 

0.139 

2.68 

120.0 

0.129 

0.046 

35.7 

0.19 

0.20 

1.039 

0.627 

0.224 

0.136 

2.61 

130.0 

0.110 

0.091 

43.2 

0.18 

0.19 

1.039 

0.605 

0.261 

C.  133 

2.55 

132.0 

0.107 

0.038 

53.8 

0.17 

0.18 

1.039 

0.580 

0.312 

0.130 

2.48 

• OFV!  AT  IONS 


TABLE  NO.  23 

DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  6 - VELOCITY  TRANSDUCER  DATA 

UNCLASSIFIED 


BIS  PA®  IS  BEST  QUALITY  PRACTICABLE 
ISMH  COPY  PUBALSHJCO  10  DDQ  ■ - 

UNCLASSIFIED  STP  453 


NO  6 3.5  IN  DIAM  9.7  PSI  3.3  IN  DIAM  EnO  PLATES 

*>FA*  OVERPRESSURE  9.7  PSI  POSITIVE  OUfiATION  233.0  MSEC 

FRIFDLANDFR  DECAY  CONSTANT  ■ 0.83 


f TIME  AFTER 

DRAG 

ABSOLUTE 

RELATIVE 

DYNAMIC 

IMPACT 

RATIO 

DRAG 

ABSOLUTE 

flow 

REYNOLDS 

f SMOCK  FRONT 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PD 

0 

01 

CO 

NO. 

R 

[ («SEC) 

(PSI  > 

(PSI  ) 

< PERCENT  1 

(PSI  ) 

IPSI ) 

(QI/OI 

(PD/QI ) 

M 

C X 10-5 1 

2.0 

1.161 

0.241 

20.7 

2.35 

2.45 

1.039 

0.473 

0.096 

0.391 

8.29 

4.0 

1.084 

0.207 

19.0 

2.28 

2.37 

1.039 

0.456 

0.087 

0.386 

8.17 

6.0 

1.012 

0.176 

17.3 

2.21 

2.30 

1.039 

0.439 

0.076 

0.362 

8.06 

8.0 

0.945 

0.148 

15.6 

2.14 

2.22 

1.039 

0.424 

0.066 

0.377 

7.94 

10.0 

0.883 

0.123 

13.9 

2.07 

2.15 

1.039 

0.409 

0.057 

0.37  2 

7.83 

12.0 

0.826 

0.100 

12.1 

2.01 

2.08 

1.039 

0.395 

0.047 

0.367 

7.7  2 

14.0 

0.774 

0.080 

10.3 

1.94 

2.02 

1.039 

0.382 

0.039 

0.363 

7.60 

16. C 

0.725 

0.063 

8.6 

1.88 

1.95 

1.039 

0.370 

0.032 

0.358 

7.49 

18.0 

0.681 

0.049 

7.1 

1*82 

1.89 

1.039 

0.359 

0.025 

0.353 

7.38 

20.0 

0.641 

0.038 

5.9 

1.76 

1.83 

1.039 

0.349 

0.020 

0.349 

7.27 

22.0 

0.604 

0.029 

4.8 

1.70 

1.77 

1.039 

0.340 

0.016 

0.344 

7.17 

2 4.0 

0.571 

0.025 

4.3 

1.65 

1.71 

1.039 

0.332 

0.014 

0.339 

7.06 

26.0 

0.54  2 

0.024 

4.4 

1.59 

1.65 

1.039 

0.326 

0.014 

0.335 

6.95 

29.0 

0.515 

0.025 

4.8 

1.54 

1.60 

1.039 

0.321 

0.015 

0.330 

6.85 

*0.0 

0.492 

0.026 

5.2 

1.49 

1.55 

1.039 

0.317 

0.016 

0.326 

6.74 

32.0 

0.471 

0.028 

5.9 

1.44 

1.49 

1.039 

0.314 

0.018 

0.321 

6.64 

34.0 

0.452 

0.029 

6.4 

1.39 

1.44 

1.039 

0.312 

0.020 

0.317 

6.53 

36.0 

C.436 

0.029 

6.6 

1.34 

1.39 

1.039 

0.311 

0.020 

0.312 

6.43 

38.0 

0.423 

0.029 

6.8 

1.30 

1.35 

1.039 

0.312 

0.021 

0.308 

6.33 

40.0 

0.411 

0.028 

6.8 

1.25 

1.30 

1.039 

0.314 

0*021 

0.304 

6.23 

42.0 

C.401 

0.026 

6.4 

1.21 

1.26 

1.039 

0.318 

0.020 

0.299 

6.13 

44.0 

C.392 

0.024 

6.1 

i»l? 

1.21 

1.039 

0.322 

0.019 

0.295 

6.03 

46.0 

0.385 

0.022 

5.7 

1.12 

1.17 

1.039 

0.327 

0.018 

0.291 

5.93 

48.0 

0.379 

0.020 

5.2 

1.09 

1.13 

1.039 

0.334 

0.017 

0.266 

5.84 

50.0 

0.374 

0.C18 

4.8 

2.05 

2.09 

1.039 

0.342 

0.-016 

0.262 

5.74 

32.0 

0.370 

0.017 

4.5 

1.01 

1.05 

1.039 

0.351 

0.016 

0.278 

5.65 

34.0 

0.366 

0.016 

4.3 

0.97 

1.01 

1.039 

0.360 

0.015 

0.2  74 

5.55 

36.0 

0.363 

0.017 

4.6 

0.94 

0.97 

1.039 

0.370 

0.017 

0.270 

5.46 

33.0 

0.359 

0.018 

5.0 

0.90 

0.94 

1.039 

0.38C 

0.019 

0.265 

5.36 

AO.O 

0.356 

0.021 

5.8 

0.87 

0.90 

1.039 

0.391 

0.023 

0.261 

5.27 

62.0 

C.353 

0.C23 

6.5 

0.84 

0.87 

1.C39 

0.4U2 

0.026 

0.257 

5.18 

6 4.0 

0.349 

0.025 

7.1 

0.81 

0.84 

1.039 

0.413 

0.029 

0.253 

5.09 

66. C 

0.345 

0.027 

7.8 

0.78 

0.81 

1.039 

0.424 

0.033 

0.249 

5.00 

68.0 

0.340 

0.028 

8.2 

0.75 

0.78 

1.039 

0.435 

0.035 

0.245 

4.91 

70.0 

0.333 

0.029 

8.7 

0.72 

0.75 

1.039 

0.442 

0.036 

0.241 

4.83 

T2.0 

0.326 

0.029 

8.8 

0.69 

0.72 

1.039 

0.45C 

0.040 

0.237 

4.74 

?4.0 

0.31? 

0.029 

9.1 

0.66 

0.69 

1.039 

0.455 

0.041 

0.233 

4.65 

76.  C 

0.307 

0.027 

8.7 

0.64 

0.66 

1.039 

0.459 

0.040 

0.229 

4.57 

*8.0 

0.295 

0.026 

B.6 

0.61 

0.64 

1.039 

0.459 

0.040 

0.225 

4.48 

• O.C 

0.281 

0.024 

0.5 

0.59 

0.61 

1.039 

0.455 

0.036 

0.221 

4.40 

*2.0 

0.264 

0.024 

9.0 

0.56 

0.59 

1.039 

0.445 

0.040 

0.218 

4.32 

*4.r 

C.246 

0.026 

10.5 

0.54 

0.56 

1.039 

0.432 

0.045 

0.214 

4.23 

*6.0 

0.22* 

0.03? 

14.2 

0.52 

0.54 

1.039 

0.41C 

0.058 

0.210 

4.15 

88.9 

0.200 

0.041 

20.5 

0.50 

0.52 

1.039 

0.362 

0.078 

0.206 

4.07 

O0.0 

0.173 

0.053 

30.6 

0.40 

0.50 

1.039 

0.344 

0.105 

0.202 

3.99 

02. 0 

0.143 

0.060 

47.5 

0.46 

0.40 

1.039 

0.297 

0.141 

0.199 

3.91 

04.0 

0.110 

0.086 

78.1 

0.44 

0.46 

1.039 

0.238 

0.186 

0.195 

3.83 

06.0 

0.073 

0.107 

146.5 

0.42 

0.44 

1.039 

0.165 

0.242 

0.191 

3.76 

08.0 

0.032 

0.130 

406.2 

0.40 

0.42 

1.039 

0.075 

0.307 

0.188 

3.68 

10C.0 

-0.012 

0.156 

-1300.0 

0.38 

0.40 

1.039 

-0.029 

0.385 

0.184 

3.60 

•*  3 STANDARD  DEVIATIONS 


nnsw®!?. 


FRO»  COP* 


TUR»1SHED 


xo  ddc 


NO  7 3.8  IN  ''JAM  9.7  PSI 

PEAK  OVERPRESSURE  9.7  PSI 


FRIFOLANDFR  RFC  AY  CONSTANT  ■ 0.83 


UNCLASSIFIED 


17.5  IN  0 1 AM  ENO  PLATES 
POSITIVE  DURATION  233.0  MSEC 


T IMF  AFTFR 

DRAG 

ARSOLUTE 

RELATIVE 

DYNAMIC 

impact 

RATIO 

DRAG 

ABSOLUTE 

FLOW 

REYNOLD 

SMOCK  FRONT 

PRESSURE 

ERROR** 

error** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

(MSFC) 

pn 
(PSI  ) 

(PSI  ) 

(PERCENT) 

Q 

(PSD 

01 

(PSI  1 

(01/0) 

CD 

(PD/OI ) 

NO. 

M 

R 

(X10-5) 

2.0 

0.800 

0.207 

29.9 

2*39 

2.49 

1.039 

0.326 

0.084 

0.391 

8.29 

4.0 

0.768 

0.187 

24.3 

2*26 

2.37 

1.039 

0.323 

0.078 

0.366 

8.17 

6.0 

0.739 

0*168 

22*8 

2*21 

2.30 

1.039 

0.321 

0.073 

0.382 

8.06 

8.0 

0.710 

0.190 

21.2 

2*14 

2.22 

1.039 

0.318 

0.067 

0.377 

7.94 

10.0 

0.682 

0.134 

19*7 

2*07 

2.19 

1.039 

0.316 

0.062 

0.372 

7.83 

12*0 

0.699 

0.119 

18.2 

2*01 

2.08 

1.C39 

C.313 

0.057 

0.367 

7.72 

14.0 

0.629 

0.109 

16.8 

1*94 

2.02 

1.039 

0.311 

0.052 

0.363 

7.60 

16.0 

0.604 

0.093 

19.4 

1*88 

1.99 

1.039 

0.3C6 

0.047 

0.358 

7,49 

1B.0 

0.981 

0.082 

14.1 

1*82 

1.89 

1.039 

0.306 

0.043 

0.353 

7,38 

20.0 

0.998 

0.072 

12*9 

1*76 

1.83 

1.039 

0.304 

0.039 

0.349 

7.27 

22.  0 

0.936 

0.063 

11*9 

1*70 

1.77 

1.039 

0.3C2 

0.036 

0.344 

7.17 

24.0 

0.914 

0.096 

11.0 

1*69 

1.71 

1.039 

0.300 

0.033 

0.339 

7.06 

26#  0 

0.493 

0.091 

10.4 

1.99 

1.69 

1.039 

0.297 

0.031 

0.335 

6.95 

28.0 

0.474 

0.047 

10«0 

1*94 

1.60 

1.039 

0.295 

0.029 

0.330 

6.85 

30.0 

0.499 

0.044 

9.8 

1*49 

1.99 

1.039 

0.293 

0.029 

0.326 

6.74 

32.0 

0.436 

0.043 

10.0 

1*44 

1.49 

1.039 

0.291 

0.029 

0.321 

6.64 

34.0 

0.419 

0.044 

10.9 

1.39 

1.44 

1.039 

0.289 

0.030 

0.317 

6.53 

36.0 

0.402 

0.044 

10.9 

1*34 

1.39 

1.039 

0.287 

0.031 

0.312 

6.43 

38.0 

0.386 

0.044 

11*4 

1*30 

1.39 

1.039 

0.285 

0.032 

0. 308 

6.33 

40.0 

0.371 

0.049 

12*1 

1*29 

1.30 

1.039 

0.284 

0.034 

0.304 

6.23 

42.0 

0.397 

0.049 

12.7 

1*21 

1.26 

1.039 

0.283 

0.036 

0.299 

6.13 

44.0 

0.343 

0.049 

13.2 

1*17 

1.21 

1.039 

0.282 

0.037 

0.295 

6.03 

46.0 

0.330 

0.049 

13*8 

1*12 

1.17 

1.039 

0.281 

0.039 

0.291 

5.93 

48.0 

0.317 

0.049 

14.2 

1.09 

1.13 

1.039 

0.260 

0.039 

0.286 

5.84 

30.0 

0.306 

0.044 

14.6 

1.09 

1.09 

1.039 

0.260 

0.041 

0.282 

5.74 

92.0 

0.299 

0.043 

14*8 

1.01 

1.09 

1.039 

0.280 

0.041 

0.278 

5.65 

34.0 

0.284 

0*042 

14.9 

0.97 

1.01 

1.039 

0.279 

0.041 

0.274 

5.55 

36.0 

0.274 

0*041 

19*0 

0.94 

0.97 

1.039 

0.280 

0.042 

C.  270 

5.46 

98.0 

0.264 

0*039 

14.9 

0.90 

0.94 

1.039 

0.280 

0.042 

0.265 

5.36 

60.0 

0.299 

0*038 

14*9 

0.87 

0.90 

1.039 

C.281 

0.042 

0.261 

5.27 

62.0 

0.247 

0.036 

14.8 

0.84 

0.87 

1.039 

0.282 

0.041 

0.257 

5.18 

64.0 

0.230 

0.C34 

14.9 

0.81 

0.84 

1.039 

0.263 

0.041 

0.253 

5.09 

66.0 

0.231 

C.033 

14.4 

0.78 

0.81 

1.039 

0.285 

0.041 

0.249 

5.00 

68.0 

0.224 

0.032 

14.2 

0.79 

0.78 

1.039 

0.287 

0.040 

0.245 

4.91 

70.0 

0.217 

0.030 

14.0 

0.72 

0.79 

1.039 

0.288 

0.040 

0.241 

4.83 

72.0 

0.210 

0.029 

14.0 

0.69 

0.72 

1.039 

0.290 

0.040 

0.237 

4.74 

74.0 

0.204 

0.028 

14.0 

0.66 

0.69 

1.039 

0.293 

0.041 

0.233 

4.65 

76.0 

0.197 

0.028 

14.3 

0.64 

0.66 

1.039 

0.296 

0.042 

0.229 

4.57 

78.0 

0.101 

0.028 

14.7 

0.61 

0.64 

1.039 

0.297 

0.043 

0.225 

4.48 

80.0 

0.189 

0.028 

19.2 

0.99 

0.61 

1.039 

0.301 

0.045 

0.221 

4.40 

•2.0 

0.180 

0.029 

16.2 

0.96 

0.59 

1.039 

0.305 

0.049 

0.218 

4.32 

84.0 

0.179 

0.030 

17.2 

0.94 

0.56 

1.039 

C.309 

0.053 

0.214 

4.23 

•6.0 

0.172 

0.030 

18.0 

0.92 

0.54 

1.039 

0.315 

0.056 

0.210 

4.15 

88.0 

0.167 

0.031 

19.0 

0.90 

0.52 

1.039 

0.320 

0.061 

0.206 

4.07 

90.0 

0.163 

0.033 

20.3 

0.48 

0.50 

1.039 

0.326 

0.066 

0.202 

3.99 

92.0 

0.160 

0.034 

21*2 

0.46 

0.48 

1.039 

0.333 

0.070 

0.199 

3.91 

94.0 

0.197 

0.039 

22. 6 

0.44 

0.46 

1.039 

0.340 

0.077 

0.195 

3.83 

06.0 

0.193 

0.036 

23.4 

0.42 

0.44 

1.039 

0.347 

0.081 

0.191 

3.76 

R8.0 

A.  1 90 

0.036 

24.4 

0.40 

0.4? 

1.039 

0.355 

0.087 

0.188 

3.68 

100.0 

A.147 

0.037 

29.3 

0.38 

0.40 

1.039 

0.364 

0.092 

0.184 

3.60 

102.0 

0.144 

0.037 

29.9 

0.37 

0.38 

1.039 

0.373 

0.096 

0.161 

3.53 

104.0 

0.141 

0.037 

26.4 

0.39 

0.36 

1.039 

0.362 

0.101 

0.177 

3.45 

106.0 

0.138 

0.037 

26.8 

0.33 

0.35 

1.039 

0.391 

0.105 

0.173 

3.38 

108.0 

0.139 

0.036 

26.9 

0.32 

0.33 

1.C39 

0.401 

0.108 

0.170 

3.31 

110.0 

0.132 

0.039 

27.1 

0.30 

0.32 

1.039 

0.410 

0.111 

0.  166 

3.23 

112.0 

0.128 

0.039 

27.6 

0.29 

0.30 

1.039 

0.419 

0.115 

0.163 

3.16 

114.0 

0.129 

0.039 

27.9 

0.28 

0.29 

1.039 

0.430 

0.120 

0.160 

3.09 

116.0 

0.122 

0.039 

28.7 

0.26 

0.27 

1.C39 

0.439 

0.126 

0.156 

3.02 

11».0 

O.lio 

0.039 

29.6 

0.29 

0.26 

1.039 

0.450 

0.133 

0.153 

2.95 

120.0 

0.116 

0.036 

31.9 

0.24 

0.25 

1.039 

0.461 

0.145 

0.  149 

2.88 

122.0 

0.114 

0.C38 

33*8 

0.23 

0.24 

1.039 

0.474 

0.160 

C.146 

2.81 

124.0 

o.ni 

0.041 

37*0 

0 .22 

0.22 

1.039 

0.487 

0.181 

0.  143 

2. 74 

126.0 

0.109 

0.044 

40.9 

0.20 

0.21 

1.039 

0.503 

0.206 

0.139 

2.68 

128.0 

0. 109 

0.049 

46.9 

0.19 

0.20 

1.039 

0.513 

0.240 

0.136 

2.61 

130.0 

0.103 

0.096 

94*4 

0.18 

0.19 

1.039 

0.527 

0.287 

0.133 

2.55 

132.0 

0.0  oo 

0.064 

64.4 

0.17 

0.18 

1.039 

0.537 

0.346 

0.130 

2.48 

«•  3 standard  DEVIATIONS 


TABLE  NO.  25 


DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  7 - VELOCITY  TRANSDUCER  DATA 

UNCLASSIFIED 


^sssssra^ 

UNCUSS  I FI  ED 


STP  453 


HO  7 S.S  IN  DIAm  9.7  PSI  17.5  tN  DIAM  ENO  PLATES 

PEAK  OVFRPRESSURE  9.7  PSI  POSITIVE  DURATION  JSJ.O  MSEC 

ERIEOLANDER  DECAY  CONSTANT  • 0.8J 


T IMF  AFTER 

DRAG 

ABSOLUTE 

RELATIVE 

DYNAMIC 

impact 

RATIO 

DRAG 

ABSOLUTE 

flow 

REYNOLDS 

SMOC<  front 

PRESSURE 

ERROR** 

ERROR** 

PRESSURE 

PRESSURE 

COEFFICIENT 

ERROR** 

MACH 

NO. 

ARRIVAL 

PD 

0 

01 

CD 

NO. 

R 

(MSFC) 

(PSI  ) 

(PSI  I 

(PERCENT! 

(PStI 

(PStl 

(01 /QJ 

(PD/Q1 ) 

M 

U10-5) 

2.0 

0.649 

0.068 

10.4 

2.35 

2.45 

1 

039 

0.264 

0.027 

0.391 

8.29 

4.0 

0.639 

0.065 

10. 1 

2.26 

2.37 

1 

039 

0.269 

0.027 

0.386 

8.17 

6*0 

0.628 

0.063 

10.0 

2.21 

2.30 

1 

039 

0.272 

0.027 

0.382 

8.06 

8.0 

0.617 

0.060 

9.7 

2.14 

2.22 

1 

039 

0.276 

0.026 

0.377 

7.94 

10. 0 

0.607 

0.057 

9.3 

2.07 

2.15 

1 

039 

0.2C1 

0.026 

0.372 

7.83 

12.0 

0.596 

0.054 

9.0 

2.01 

2.08 

1 

039 

0.285 

0.025 

0.367 

7.72 

14.0 

0.585 

0.051 

8.7 

1.94 

2.02 

1 

039 

0.289 

0.025 

0.363 

7.60 

16.0 

0.575 

0.048 

8.3 

1.88 

1.95 

1 

039 

0.293 

0.024 

0.358 

7.49 

10.0 

0.564 

0.045 

7.9 

1.82 

1.89 

1 

039 

0.297 

0.023 

0.353 

7.38 

20.0 

0.553 

0.043 

7.7 

1.76 

1.83 

1 

039 

0.301 

0.023 

0.349 

7.27 

22. 0 

0.542 

0.040 

7.3 

1.70 

1.77 

1 

039 

0.305 

0.022 

0.344s 

7.17 

24.0 

0.532 

0.037 

6.9 

1.65 

1.71 

1 

039 

0.310 

0.021 

0.339 

7.06 

26.0 

0.521 

0.034 

6.5 

1.59 

1.65 

1 

039 

0.314 

0.020 

0.335 

6.95 

28.0 

0.510 

0.031 

6.0 

1.54 

1.60 

1 

039 

0.317 

0.019 

0.330 

6.85 

30.0 

0.500 

0.029 

5*8 

1.49 

1.55 

1 

039 

0.322 

0.C16 

0.326 

6.74 

32.0 

0.489 

0.026 

5.3 

1.44 

1.49 

1 

039 

0.326 

0.017 

0.321 

6.64 

34.0 

0.478 

0.023 

4.8 

1.39 

1.44 

1 

039 

0.330 

0.015 

0.317 

6.53 

36.0 

0.468 

0.021 

4.4 

1.34 

1.39 

1 

039 

0.334 

0.015 

0.312 

6.43 

38.0 

0.457 

0.019 

4.1 

1.30 

1.35 

1 

039 

0.338 

0.014 

0.308 

6.33 

40.0 

0.446 

0.016 

3.5 

1.25 

1.30 

1 

039 

0.341 

0.C12 

0.304 

6.23 

42.0 

0.436 

0.014 

3.2 

1.21 

1.26 

1 

039 

0.345 

0.011 

0.299 

6.13 

44.0 

0.425 

0.013 

3.0 

1.17 

1.21 

1 

039 

0.349 

0.C10 

0.295 

6.03 

46.0 

0.414 

0.012 

2.8 

1.12 

1.17 

1 

039 

0.352 

0.010 

0.291 

5.93 

48.0 

0.404 

0.011 

2.7 

1.09 

1.13 

1 

039 

0.356 

0*009 

0.286 

5.84 

50.0 

0.393 

0.012 

3.0 

1.05 

1.09 

1 

039 

0.359 

0.010 

0.282 

5.74 

52.0 

0.382 

0.013 

3.4 

i.Oi 

1.05 

1 

039 

0.362 

0.012 

0.278 

5.65 

54.0 

0.371 

0.014 

3.7 

0.97 

1.01 

1 

039 

0.365 

0.013 

0.274 

5.55 

56.0 

0.361 

0.016 

4.4 

0.94 

0.97 

1 

039 

0.366 

0.016 

0.270 

5.46 

58.0 

0.350 

0.018 

5.1 

0.90 

0.94 

1 

039 

0.370 

0.019 

0.265 

5.36 

6U.0 

0.339 

0.021 

6.1 

0.87 

0.90 

1 

039 

0.372 

0.023 

0.261 

5.27 

*2.0 

0.329 

0.023 

6.9 

0.84 

0.87 

1 

039 

0.375 

0.026 

0.257 

5.18 

*4.0 

0.318 

0.026 

8.1 

0.91 

0.84 

1 

039 

0.376 

0.030 

0.253 

5.09 

*6.0 

0.307 

0.028 

9.1 

0.78 

0.81 

1 

039 

0.378 

0.034 

0.249 

5.00 

68.0 

0.297 

0.031 

10.4 

0.75 

0.78 

1 

039 

0.380 

0.039 

0.245 

4.91 

70.0 

0.286 

0.034 

11.8 

0.72 

0.75 

1 

039 

0.380 

0.045 

0.241 

4.83 

72.0 

0.275 

0.036 

13.0 

0.69 

0.72 

1 

039 

0.380 

0.049 

0.237 

4.74 

74.0 

0.265 

0.039 

14.7 

0.66 

0.69 

1 

039 

0.381 

0.056 

0.233 

4.65 

76.0 

0.254 

0.042 

16.5 

0.64 

0.66 

1.039 

0.380 

0.062 

0.229 

4.57 

78.0 

0.243 

0.045 

18.5 

0.61 

0.64 

1.039 

0.378 

0.070 

0.225 

4.48 

AO.  0 

0.233 

C.04R 

20.6 

0.59 

0.61 

1.039 

0.377 

0.077 

0.221 

4.40 

82.0 

0.222 

0.051 

22. 9 

0.56 

0.59 

1.039 

0.374 

0.086 

0.218 

4.32 

•4.0 

0.211 

0.053 

25.1 

0.54 

0.56 

1.039 

0.371 

0.093 

0.214 

4.23 

86.0 

0.200 

0.056 

28.0 

0.52 

0.54 

1.039 

0.366 

0.102 

0.210 

4.15 

• 0.0 

0.190 

0.059 

31.0 

0.50 

0.52 

1.039 

0.363 

0.112 

0.206 

4.07 

90.0 

0.179 

0.062 

34.6 

0.48 

0.50 

1.039 

0.356 

0.123 

0.202 

3.99 

92.0 

0.168 

0.065 

38.6 

0.46 

0.48 

1.039 

0.349 

0.135 

0.199 

3.91 

94.0 

0.158 

0.060 

43.0 

0.44 

0.46 

1.039 

0.342 

0.147 

0.195 

3.83 

96.0 

0.147 

0.071 

48.2 

0.42 

0.44 

1.039 

0.333 

0.160 

0.191 

3.76 

90.0 

0.136 

0.074 

54.4 

0.40 

0.42 

1.039 

0.321 

0.175 

0.188 

3.68 

10C.0 

0.126 

0.077 

61.1 

0.38 

0.40 

1.039 

0.311 

0.190 

0.184 

3.60 

102.0 

0.115 

0.079 

68.6 

0.37 

0.38 

1.039 

0.297 

0.204 

0.181 

3.53 

104.0 

0.104 

C.082 

78.8 

0.35 

0.36 

1.039 

0.281 

0.221 

0.177 

3.45 

••  1 STANDARD  DEVIATIONS 

TABLE  NO.  26 

DRAG  COEFFICIENT  VS  TIME  FOR  CYLINDER  7 - CAMERA  DATA 


UNCLASSIFIED 


UNCLASSIFIED 


TABLE  NO.  27 

SUMMARY  OF  RESULTS  FOR  EACH  CYLINDER 


Cylinder  - f n t , Useful  Data 

Number  Type  of  0ata  Obtained  Remarks 

1 V.T.  (East)  No 

V.T.  (West)  No 

Cam.  (Westj  No 

Cylinder  failed  to  undergo  free  flight  due  to  sidewise  blast 
anomaly  at  20  psi  peak  overpressure  location  acting  on 
. large  end  plates. 

^ u’t"  i^asti  ^o  _ ^ast  ^net  broke  shortly  after  shock  arrival. 

r ' Yes  ' Poor  signal/noise  due  to  error  in  circuit  controlling  sensitivity. 

Cam.  (West)  No  - Non-constant  film  speed;  violent  displacement  of  camera  post. 

3 V.T.  (East)  Yes 

V.T.  (West)  Yes 

Cam.  (West)  Yes 

Both  - - 

(Large  amplitude  79  Hz  oscillations  on  signal  produced  large 

L uncertainties  in  derived  drag  pressure. 

Oscillations  absent;  smaller  uncertainties  than  for  V.T.  data. 
East-West  V.T.  results  consistent  over  range  of  measurement.  V.T. 
and  Cam.  results  consistent  over  entire  range  of  measurement. 

4 V.T.  (East)  Yes 

V.T.  (West)  Yes  *| 

Cam.  (West)  Yes 

Both  - - 

(Moderately  large  60  Hz  oscillations  on  signals  produced  increased 

L uncertainties  in  derived  drag  pressure. 

Excellent  data;  only  small  corrections  for  camera  motion. 

East-West  V.T.  results  consistent  over  range  of  measurement.  V.T. 
and  Cam.  data  consistent  over  most  of  range  of  measurement. 

5 V.T.  (East)  Yes 

V.T.  (West)  Yes 

Cam.  (West)  Yes 

Both  - 

Moderately  large  60  Hz  oscillations  on  signals;  signals  terminated 
prematurely  due  to  contact  of  cylinder  shaft  with  photomarker 
- plate. 

Excellent  data;  only  small  corrections  for  camera  motion. 

Linear  drag  pressure  from  V.T.  data  in  agreement  over  most  of 
range  with  (approximately)  linear  drag  pressure  from  higher 
order  fit  to  Cam.  data. 

6 V.T.  (East)  Yes  f 

V.T.  (West)  Yes 

Cam.  (West)  Yes 

Both  - - 

Smaller  59  Hz  oscillations  with  irrigular  fluctuations  superposed. 

Excellent  data;  larger  corrections  for  camera  motion  than  for 

Cyl.  4,  5. 

East-West  V.T.  results  consistent  over  entire  range  of  measurement. 
V.T.  and  Cam.  data  consistent  over  most  of  range  of  measurement. 

7 V.T.  (East)  Yes 

V.T.  (West)  Yes  'L 

Cam.  (West)  Yes 

Both 

Smaller  59  Hz  oscillations  with  irregular  fluctations  superposed. 

Uncertainties  in  drag  pressure  competitive  with  Cam.  data. 

Large  corrections  for  camera  post  motion  forced  linear  fit  to  drag 
pressure. 

East-West  V.T.  results  consistent  over  entire  range  of  measurement. 
V.T.  and  Cam.  data  consistent  over  most  of  range  of  measurement. 

* 


V.T.  - Velocity  Transducer.  Cam.  - High-Speed  Camera. 
(East-West)  - refers  to  end  of  cylinder  where  measurement  recorded. 
Both  - Comments  refer  to  both  V.T.  and  Cam.  data. 


UNCLASSIFIED 


UNCLASSIFIED 


TABLE  NO.  28 

SUMMARY  OF  BEST  VALUES  FOR  DRAG  COEFFICIENT* 


Cyl inder 
Diameter 
(inches) 

Time  After 
Shock  Arrival 
(msec) 

Average 

Reynolds  Number 

Average 

Drag  Coefficient 

3.5 

3-25 

7.66 

0.274  ± 0.054** 

25-50 

6.34 

0.255  ± 0.042 

50-75 

5.14 

0.251  ± 0.049 

75-100 

4.08 

0.282  ± 0.065 

100-125 

3.13 

0.382  ± 0.148 

9.5 

3-25 

15.5 

0.508  ± 0.106 

25-50 

13.3 

0.482  ± 0.085 

50-75 

11.3 

0.455  ± 0.106 

75-100 

9.4 

0.401  ± 0.160 

18.0 

3-25 

39.4 

0.420  ± 0.121 

25-50 

32.6 

0.449  ± 0.099 

50-75 

26.5 

0.625  ± 0.137 

Data  plotted  in  Figure  28. 

3 Standard  deviations  of  uncertainty  (99%  confidence  interval) 
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FIG. I C0(M)  FDR  DIFFERENT  L/D  8 R VALUES 
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DIAMETER  - LENGTH  RATIO , 0/  L 
FIG. 2 Cn  (0/L)  FOR  DIFFERENT  (M , R ) VALUES 


FIG.  3 RANGE  OF  (M,R)  VALUES  MCLUOEO  IN  ALL  PREVIOUS  MEASUREMENTS  & IN  DICE  THROW 


FIG.  4 CYLINDER  DRAG  ASSEMBLY  - LOCATION  1 
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DRAG  CYLINDER  A (9.5  INCH  DIAMETER) 

DRAG  CYLINDER  B (3.5  INCH  DIAMETER) 

DRAG  CYLINDER  C (18.0  INCH  DIAMETER) 

LATTICE  MAST 

UHF  POLEMAST 

WHIP  ANTENNA 

CAMERA  POSITIONS 

PRESSURE  TRANSDUCER 

DUST  COLLECTOR 

CAMERA  LINE  OF  SIGHT 
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FIG.  6 LAYOUT  OF  CYLINDER  DRAG  PROJECT 
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FIG.  12  EAST  SUPPORT  FOR  CYLINDER  1 (POST-SHOT) 
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CYLINDER  5 (POST-SHOT)  AT  7 PSI  PEAK  OVERPRESSURE 
LOCATION  SHOWING  EXTENT  OF  LATERAL  DISPLACEMENT 
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FIG.  19  CYLINDER  4 - HIGH  SPEED  CAMERA  DATA 
(6.7  PSI  PEAK  OVERPRESSURE) 
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FIG.  20  CYLINDER  5 - VELOCITY  TRANSDUCER  DATA 
(6.7  PSI  PEAK  OVERPRESSURE) 
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FIG.  21  CYLINDER  5 • HIGH  SPEED  CAMERA  DATA 
(6.7  PSI  PEAK  OVERPRESSURE) 
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FIG.  25  CYLINDER  7 - HIGH  SPEED  CAMERA  DATA 
(9.7  PSI  PEAK  OVERPRESSURE) 
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FIG.  28  CYLINDER  5 - DRAG  COEFFICIENT  VS  TIME 

Top  - Velocity  Transducer  Data 
Bottom  - Camera  Data 
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FIG.  33  MEASURED  DRAG  COEFFICIENT  VS  REYNOLDS  NUMBER  - 3.5  INCH  DIAMETER  CYLINDER 
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FIG.  34  MEASURED  DRAG  COEFFICIENT  VS  REYNOLDS  NUMBER 
9.5  INCH  DIAMETER  CYLINDER 
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FIG.  35  MEASURED  DRAG  COEFFICIENT  VS  REYNOLDS  NUMBER 
18.0  INCH  DIAMETER  CYLINDER 
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FIG.  36  MEASURED  DRAG  COEFFICIENT  VS  REYNOLDS  NUhBER 
COMPOSITE  PLOT  FOR  ALL  CYLINDER  DIAMETERS 


UNCLASSIFIED 


STP  453 


in 

O 

X 

CL 

CL 

Ul 

00 


o 

z 

>- 

Ui 

CL 


I 


°0  1N3IDIJJ303  9V«a 
UNCLASSIFIED 


FIG.  38  DRAG  COEFFICIENT  VS  REYNOLDS  NUMBER 
BEST  VALUES  (AFTER  CORRECTIONS) 
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